A stope overburden load three-dimensional model construction method

Abstract

The application discloses a stope overburden load three-dimensional model construction method, relates to the model construction technical field, and comprises the following steps: obtaining a sample data scatter diagram based on sample data of an overburden to be detected; obtaining a first stress threshold value and a second stress threshold value based on the sample data scatter diagram; obtaining a stress representation value based on the first stress threshold value and the second stress threshold value; obtaining a stress change function based on the stress representation value; obtaining generated data points based on the overburden to be detected; obtaining a first detection point based on the generated data points; obtaining a real-time stress value based on the first detection point and the stress change function; obtaining a function fitting stress based on the stress change function and the real-time stress value; and generating a stope overburden load three-dimensional model based on the function fitting stress. The application is used to solve the problem of low stress acquisition accuracy of an instrument at an undetected position in the prior art model construction technology, thereby improving the accuracy of the stope overburden load three-dimensional model construction.

Description

Technical Field

[0001] This invention relates to the field of model building technology, specifically to a method for constructing a three-dimensional model of overburden load in a mining area. Background Technology

[0002] In the underground mining of coal resources, the overburden load of the mining area is a key factor affecting the stability of the surrounding rock, the reliability of the support system, and the safe production of the mine. With the increase of mining depth and mining intensity, the stress environment of the overburden in the mining area becomes increasingly complex, so more detailed three-dimensional stress modeling of the overburden is required. Existing stress detection methods cannot detect stress at every location; they can only detect stress at a subset of points. Stress values ​​not at the detection points are calculated by linearly varying the stress between two detection points. However, stress variations differ between different types of rock strata and depths, leading to inaccurate stress calculations. In other words, existing model-building techniques have low accuracy in obtaining stress at locations not detected by the instrument, resulting in low accuracy in constructing three-dimensional models of overburden loads in mining areas. Summary of the Invention

[0003] This invention aims to at least partially solve one of the technical problems in the prior art. It involves obtaining a scatter plot of sample data based on sample data of the overburden to be detected; obtaining a first stress threshold and a second stress threshold based on the scatter plot; obtaining stress representation values ​​based on the first and second stress thresholds; obtaining a stress variation function based on the stress representation values; obtaining generated data points from the overburden to be detected; obtaining a first detection point based on the generated data points; obtaining a real-time input value based on the first detection point and the stress variation function; obtaining a function-fitted stress based on the stress variation function and the real-time input value; and generating a three-dimensional model of the overburden load in the stope based on the function-fitted stress. This addresses the problem of low accuracy in obtaining stress at locations not detected by the instrument in existing model construction techniques, leading to low accuracy in constructing three-dimensional models of the overburden load in the stope.

[0004] To achieve the above objectives, this application provides a method for constructing a three-dimensional model of overburden load in a stope, comprising the following steps: A scatter plot of sample data was obtained based on the sample data of the overlying rock to be detected; The first and second stress thresholds are obtained based on the scatter plot of the sample data. The first stress threshold and the second stress threshold are used to obtain the stress representation value; Obtain the stress variation function based on the stress representation value; Acquire the overlying rock to be detected and generate data points; The first detection point is obtained based on the generated data points; The real-time input value is obtained based on the first detection point and the stress change function; The stress is obtained by fitting the function based on the stress variation function and the real input value, and a three-dimensional model of the overburden load in the mining area is generated based on the stress fitted by the function.

[0005] Furthermore, obtaining a scatter plot of sample data based on the sample data of the overburden to be detected includes the following sub-steps: The first number of rock layer stresses were obtained in the same area, the same rock layer, and the same overburden depth, and marked as sample stresses. A Cartesian coordinate system is established with sample stress as the horizontal axis data and the quantity of sample stress as the vertical axis data, and this system is labeled as the sample data coordinate system. Obtain the sample stress and the corresponding quantity as the coordinate points of the x and y coordinates, and mark them as sample data coordinate points; Plot all the coordinate points of the sample data on the sample data coordinate system to obtain a scatter plot, and label it as the sample data scatter plot.

[0006] Furthermore, obtaining the first stress threshold and the second stress threshold based on the scatter plot of the sample data includes the following sub-steps: Get the maximum value of the ordinate among all sample data coordinate points and mark it as the sample data height; Set a length value and mark it as the rectangle length setting value; Create a rectangle on the horizontal axis of the sample data coordinate system with a height equal to the sample data height and a width equal to the set rectangle length, and be able to move left and right. Mark this rectangle as the set search rectangle. The total number of sample stresses corresponding to the coordinate points of the sample data included within the set search rectangle is marked as the rectangle search quantity; The range length of the sample stress is marked as the sample data length; Assuming the sample stress is uniformly distributed within the range of sample stress, the number of rectangle judgments at this time is obtained and marked as the average number of rectangles; Set a ratio, which is marked as the first set ratio; obtain the product of the average number of rectangles and the first set ratio, and mark it as the abnormal number threshold.

[0007] Furthermore, obtaining the first stress threshold and the second stress threshold based on the scatter plot of the sample data also includes the following sub-steps: In the scatter plot of the sample data, the set search rectangle is shifted from the leftmost side to the right. When the number of rectangles searched is greater than or equal to the anomaly threshold, the movement of the set search rectangle is stopped. The sample stress corresponding to the minimum x-coordinate of the set search rectangle at this time is obtained and marked as the first stress threshold. In the scatter plot of the sample data, the set search rectangle is shifted from the rightmost side to the left. When the number of rectangles searched is greater than or equal to the anomaly threshold, the movement of the set search rectangle is stopped. The sample stress corresponding to the largest x-coordinate of the set search rectangle at this time is obtained and marked as the second stress threshold.

[0008] Furthermore, obtaining the stress representation value from the first stress threshold and the second stress threshold includes the following sub-steps: The sample stress between the first stress threshold and the second stress threshold is marked as the screening stress; Obtain the mean value of the screening stress and mark it as the stress representation value; Obtain stress values ​​for the same region, the same stratum, and different overburden depths.

[0009] Furthermore, obtaining the stress change function based on the stress representation value includes the following sub-steps: A plane rectangular coordinate system is established with overburden depth as the horizontal axis data and stress as the vertical axis data, and it is marked as the stress change coordinate system. Obtain the overburden depth and the corresponding stress values ​​on the x-axis and y-axis, and mark them as stress change coordinate points; Plot all stress variation coordinate points in the stress variation coordinate system; The function is obtained by fitting a function to all the stress change coordinate points, and it is labeled as the stress change function.

[0010] Furthermore, obtaining the data points from the overlying rock to be detected includes the following sub-steps: Obtain the vertically downward detection line from the ground detection point and mark it as the real-time detection line; Obtain stress detection points on the real-time detection line and mark them as real-time detection points; Obtain a three-dimensional rock stratum model of the overlying rock to be detected and mark it as a real-time three-dimensional rock stratum model; Obtain any point within the real-time 3D model of the same rock stratum and mark it as a generated data point.

[0011] Furthermore, obtaining the first detection point based on the generated data points includes the following sub-steps: Get the nearest point on the real-time detection line with the same height as the generated data point and mark it as the line reference point; Obtain the real-time detection points adjacent to the line reference point on both sides and mark them as adjacent detection points; Obtain the stress values ​​of adjacent detection points and mark them as transformed stress; The adjacent detection points corresponding to the smaller transformation stress between the two transformation stresses are marked as the first detection point.

[0012] Furthermore, obtaining the real-time input value based on the first detection point and the stress change function includes the following sub-steps: Obtain the points of transformed stress on the stress change function and mark them as function reference points; Obtain the stress change function between two function reference points and mark it as the function reference segment; Get the length of the function reference segment and mark it as the function reference length; Obtain the length of two adjacent detection points and mark it as the detection interval length; The detection interval length is scaled proportionally to match the function reference length, and the distance from the first detection point to the straight line reference point is obtained and marked as the first offset distance. Obtain the x-coordinate of the function reference point with the smaller x-coordinate, and mark it as the first starting distance; Get the sum of the first offset distance and the first starting distance, and mark it as the real-time input value.

[0013] Furthermore, the stress is fitted using a function based on the stress variation function and real-time input values. The generation of a three-dimensional model of the overburden load in the stope based on this fitted stress includes the following sub-steps: Substitute the real-time input value as the horizontal axis value into the stress variation function to obtain the vertical value, which is then labeled as the function-fitted stress. The stress is fitted to the function to generate the stress at the data points; Obtain the stress of all generated data points within the real-time 3D rock strata model to generate a 3D model of overburden load in the mining area.

[0014] The beneficial effects of this invention are as follows: This invention obtains a scatter plot of sample data based on sample data of the overburden to be detected; obtains a first stress threshold and a second stress threshold based on the scatter plot; obtains a stress representation value based on the first stress threshold and the second stress threshold; obtains a stress change function based on the stress representation value; obtains data points of the overburden to be detected; obtains a first detection point based on the generated data points; obtains a real-time input value based on the first detection point and the stress change function; obtains a function-fitted stress based on the stress change function and the real-time input value; and generates a three-dimensional model of the overburden load in the stope based on the function-fitted stress. The advantage lies in improving the accuracy of stress acquisition at locations not detected by the instrument and improving the accuracy of constructing a three-dimensional model of the overburden load in the stope. This invention obtains the stress change function based on the stress representation value. Its advantage lies in obtaining the stress change relationship of different rock layers and improving the accuracy of constructing a three-dimensional model of the overburden load in the mining area. Attached Figure Description

[0015] Figure 1 This is a flowchart illustrating the steps of the method of the present invention; Figure 2 This is a schematic diagram of the first stress threshold and the second stress threshold of the present invention; Figure 3 This is a schematic diagram of the stress variation function of the present invention; Figure 4 This is a schematic diagram of the function reference points of the present invention. Detailed Implementation

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

[0017] Example 1, please refer to Figure 1 As shown, this application provides a method for constructing a three-dimensional model of overburden load in a stope, including the following steps: Step S1: Obtain a scatter plot of sample data based on the sample data of the overburden to be detected; Step S1 includes the following sub-steps: Step S101: Obtain a first number of rock layer stresses in the same area, the same rock stratum, and the same overburden depth, and mark them as sample stresses; the existing underground contains different rock strata, and the stress changes of different rock strata are different. At the same time, the data is collected at points, which makes it difficult to obtain the stress at all points; the stress is vertical stress, which indicates that the surrounding rock is not easy to yield and the support is simple; in order to obtain the range of sample stresses, the first number should not be too small, for example, the first number is 1000; Step S102: Establish a Cartesian coordinate system with sample stress as the horizontal axis data and the number of sample stresses as the vertical axis data, and mark it as the sample data coordinate system; Step S103: Obtain the sample stress and the corresponding quantity, with the coordinate points as the x-axis and y-axis respectively, and mark them as sample data coordinate points; Step S104: Plot all the coordinate points of the sample data on the sample data coordinate system to obtain a scatter plot, and mark it as the sample data scatter plot. For practical applications, please refer to Figure 2 As shown, this is a scatter plot of sample data obtained when the same overburden depth is 200m.

[0018] Step S2 involves obtaining a first stress threshold and a second stress threshold based on the scatter plot of the sample data. Step S2 includes the following sub-steps: Step S201: Obtain the maximum value of the ordinate among all sample data coordinate points and mark it as the sample data height; Step S202: Set a length value, marked as the rectangle length setting value; setting the rectangle length setting value is for constructing the search rectangle, and setting the search rectangle is for analyzing the stress distribution of the sample; the rectangle length setting value should not be too large, for example, the rectangle length setting value is 0.1; Step S203: Create a rectangle on the horizontal axis of the sample data coordinate system with a height equal to the sample data height and a width equal to the rectangle length setting value, and mark it as the set search rectangle; Step S204: The total number of sample stresses corresponding to the coordinate points of the sample data included within the set search rectangle is marked as the number of rectangles to search. Step S205: Mark the range length of the sample stress as the sample data length; Step S206: Assuming that the sample stress is uniformly distributed within the range of sample stress, obtain the number of rectangle judgments at this time and mark it as the average number of rectangles. In practical applications, the length of the obtained sample data is 3.2, and the average number of rectangles is: 1000 × (0.1 ÷ 3.2) = 31.25; where 1000 is the first quantity and 0.1 is the set value of the rectangle length. Step S207: Set a ratio, marked as the first set ratio; obtain the product of the average number of rectangles and the first set ratio, marked as the abnormal number threshold; the first set ratio is to filter the range of stress of samples with too small a distribution number, so the first set ratio should not be set too small, for example, the first set ratio is 0.2. In practical applications, the threshold for obtaining the number of anomalies is: 31.25 × 2 = 6.25.

[0019] Step S208: In the scatter plot of the sample data, the set search rectangle is shifted from the leftmost side to the right. When the number of rectangles searched is greater than or equal to the abnormal number threshold, the movement of the set search rectangle is stopped. The sample stress corresponding to the minimum horizontal coordinate of the set search rectangle at this time is obtained and marked as the first stress threshold. In order to filter out excessively small sample stresses and obtain the accurate minimum sample stress value. Step S209: In the scatter plot of the sample data, the set search rectangle is shifted from the rightmost side to the left. When the number of rectangles searched is greater than or equal to the abnormal number threshold, the movement of the set search rectangle is stopped. The sample stress corresponding to the largest horizontal coordinate of the set search rectangle at this time is obtained and marked as the second stress threshold. In order to filter out excessive sample stress, the accurate maximum value of sample stress is obtained. In practical applications, in a scatter plot of sample data, the search rectangle is shifted from the leftmost position to the right. The shift stops when the number of rectangles searched is greater than or equal to 6.25. (See [link to relevant documentation]). Figure 2As shown, stop moving the search rectangle to its current position; obtain the sample stress corresponding to the minimum x-coordinate of the search rectangle, which is 4.1, then the first stress threshold is 4.1; in the sample data scatter plot, shift the search rectangle from the rightmost position to the left. When the number of rectangles searched is greater than or equal to 6.25, stop moving the search rectangle. Please refer to [link to relevant documentation]. Figure 2 As shown, the sample stress corresponding to the largest horizontal coordinate of the search rectangle at this time is 7.1, and the second stress threshold is 7.1.

[0020] Step S3: Obtain stress representation values ​​from the first stress threshold and the second stress threshold; Step S3 includes the following sub-steps: Step S301: The sample stress between the first stress threshold and the second stress threshold is marked as the screening stress; Step S302: Obtain the mean value of the screening stress and mark it as the stress representation value; Step S303: Obtain stress values ​​for the same region, the same stratum, and different overburden depths; In practical applications, for example, when the same area, the same rock stratum, and the overburden depth is 200m, the sample stress between 4.1 and 7.1 is the screening stress.

[0021] Step S4: Obtain the stress change function based on the stress representation value; Step S4 includes the following sub-steps: Step S401: Establish a plane rectangular coordinate system with overburden depth as the horizontal axis data and stress as the vertical axis data, and mark it as the stress change coordinate system; Step S402: Obtain the overburden depth and the corresponding stress values ​​as the x-axis and y-axis coordinate points, and mark them as stress change coordinate points; Step S403: Plot all stress change coordinate points in the stress change coordinate system; Step S404: Fit a function to all stress change coordinate points to obtain a function, and mark it as the stress change function; obtain the stress change region with depth to facilitate the acquisition of strain at each point; For practical applications, please refer to Figure 3 As shown, a stress variation function was obtained.

[0022] Step S5: Obtain the overburden to be detected and generate data points; Step S5 includes the following sub-steps: Step S501: Obtain the straight line detected vertically downward from the ground detection point and mark it as the real-time detection line; Step S502: Obtain the stress detection points on the real-time detection line and mark them as real-time detection points; Step S503: Obtain the three-dimensional rock strata model of the overlying rock to be detected and mark it as the real-time three-dimensional rock strata model; Step S504: Obtain any point within the real-time 3D model of the same rock stratum and mark it as a generated data point. The generated data point is a 3D data point.

[0023] Step S6: Obtain the first detection point based on the generated data points; Step S6 includes the following sub-steps: Step S601: Obtain the point with the same height on the real-time detection line closest to the generated data point and mark it as the line reference point; convert the 3D data point into a point on the real-time detection line to facilitate data analysis; Step S602: Obtain the real-time detection points adjacent to the line control point on both sides and mark them as adjacent detection points; Step S603: Obtain the stress values ​​of adjacent detection points and mark them as transformed stress; Step S604: Mark the adjacent detection points corresponding to the smaller transformation stress between the two transformation stresses as the first detection point.

[0024] Step S7: Obtain the real-time input value based on the first detection point and the stress change function; Step S7 includes the following sub-steps: Step S701: Obtain the point of transformed stress on the stress change function and mark it as the function reference point; Step S702: Obtain the stress change function between two function reference points and mark it as the function reference segment; Step S703: Obtain the length of the function reference segment and mark it as the function reference length; the function reference length is the length of the horizontal axis. Step S704: Obtain the length of two adjacent detection points and mark it as the detection interval length; Step S705: Scale the detection interval length proportionally to match the function reference length, obtain the distance from the first detection point to the straight line reference point at this time, and mark it as the first offset distance; Step S706: Obtain the x-coordinate of the function reference point with the smaller x-coordinate and mark it as the first starting distance; Step S707: Obtain the sum of the first offset distance and the first starting distance, and mark it as the real-time input value; because the collected data is limited and there is data offset, the data change is obtained through the stress change function and converted into depth; For practical applications, please refer to Figure 4 As shown, the function reference point is obtained with a detection interval length of 50m. The detection interval length is scaled proportionally to match the function reference length. The distance from the first detection point to the straight line reference point is 25m, so the first offset distance is 25m. The x-coordinate of the function reference point with the smaller x-coordinate is 200m, so the first starting distance is 200m. The sum of the first offset distance and the first starting distance is 225m, so the actual input value is 225m.

[0025] Step S8 involves obtaining a function to fit stress based on the stress variation function and real-time input values, and generating a three-dimensional model of the overburden load in the stope based on the stress fitted by the function. Step S8 includes the following sub-steps: Step S801: Substitute the real input value as the horizontal axis value into the stress change function to obtain the vertical value, and mark it as the function fitting stress; Step S802: Fit the stress of the function to the stress of the generated data points; Step S803: Obtain the stress of all generated data points in the real-time three-dimensional rock strata model to generate a three-dimensional model of the overburden load in the mining area. For practical applications, please refer to Figure 3 As shown, 225m is substituted into the stress variation function as the horizontal axis value to obtain a vertical value of 6MPa. 6MPa is then used as the stress of this generated data point. Therefore, this method can be used to obtain the stress of each data point in the real-time three-dimensional rock strata model, so that even with fewer data points collected, the stress of all points in the real-time three-dimensional rock strata model can be obtained. The stress can be used to determine whether the overburden load of the mining area is operating safely.

[0026] Example 2: This application also provides an electronic device, which may include: a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus. The memory stores computer-readable instructions, and the processor can call the instructions in the memory. When the computer-readable instructions are executed by the processor, the steps of a method for constructing a three-dimensional model of overburden load in a mining area are performed to achieve the following functions: obtaining a scatter plot of sample data based on sample data of the overburden to be detected; obtaining a first stress threshold and a second stress threshold based on the scatter plot of sample data; obtaining stress representation values ​​based on the first stress threshold and the second stress threshold; obtaining a stress change function based on the stress representation values; obtaining generated data points of the overburden to be detected; obtaining a first detection point based on the generated data points; obtaining a real-time input value based on the first detection point and the stress change function; obtaining a function-fitted stress based on the stress change function and the real-time input value; and generating a three-dimensional model of overburden load in the mining area based on the function-fitted stress.

[0027] Furthermore, when the logical instructions in the aforementioned memory can be implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0028] Example 3: This application also provides a computer program product, which includes a computer program stored on a computer-readable storage medium. The computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute a method for constructing a three-dimensional model of overburden load in a stope, as provided by the methods described above. This method includes: obtaining a scatter plot of sample data based on sample data of the overburden to be detected; obtaining a first stress threshold and a second stress threshold based on the scatter plot of sample data; obtaining a stress representation value based on the first stress threshold and the second stress threshold; obtaining a stress variation function based on the stress representation value; obtaining generated data points of the overburden to be detected; obtaining a first detection point based on the generated data points; obtaining a real-time input value based on the first detection point and the stress variation function; obtaining a function-fitted stress based on the stress variation function and the real-time input value; and generating a three-dimensional model of overburden load in a stope based on the function-fitted stress.

[0029] Example 4: This application also provides a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, it performs the steps in the above-described method for constructing a three-dimensional model of overburden load in a stope, to achieve the following functions: obtaining a scatter plot of sample data based on sample data of the overburden to be detected; obtaining a first stress threshold and a second stress threshold based on the scatter plot of sample data; obtaining stress representation values ​​based on the first stress threshold and the second stress threshold; obtaining a stress variation function based on the stress representation values; obtaining generated data points of the overburden to be detected; obtaining a first detection point based on the generated data points; obtaining a real-time input value based on the first detection point and the stress variation function; obtaining a function-fitted stress based on the stress variation function and the real-time input value; and generating a three-dimensional model of overburden load in the stope based on the function-fitted stress.

[0030] Based on the above description of the embodiments, the embodiments of the present invention can be provided as methods, systems, or computer program products. Based on this understanding, the above technical solutions, in essence or in terms of their contribution to the prior art, 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., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or certain parts of the embodiments.

[0031] In the embodiments provided in this application, it should be understood that the disclosed system or method can be implemented in other ways. The embodiments described above are merely illustrative. For example, the division of modules or units is only a logical functional division, and there may be other division methods in actual implementation. Furthermore, multiple modules or units may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the coupling or direct coupling or communication connection shown or discussed may be through some communication interfaces. The indirect coupling or communication connection between systems, modules, and units may be electrical, mechanical, or other forms.

[0032] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of this application, and are not intended to limit them. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application.

Claims

1. A method for constructing a three-dimensional model of overburden load in a stope, characterized in that, Includes the following steps: A scatter plot of sample data was obtained based on the sample data of the overlying rock to be detected; The first and second stress thresholds are obtained based on the scatter plot of the sample data. The first stress threshold and the second stress threshold are used to obtain the stress representation value; Obtain the stress variation function based on the stress representation value; Acquire the overlying rock to be detected and generate data points; The first detection point is obtained based on the generated data points; The real-time input value is obtained based on the first detection point and the stress change function; The stress is obtained by fitting the function based on the stress variation function and the real input value, and a three-dimensional model of the overburden load in the mining area is generated based on the stress fitted by the function.

2. The method for constructing a three-dimensional model of overburden load in a stope according to claim 1, characterized in that, Obtaining a scatter plot of sample data based on the sample data of the overburden to be detected includes the following sub-steps: The first number of rock layer stresses were obtained in the same area, the same rock layer, and the same overburden depth, and marked as sample stresses. A Cartesian coordinate system is established with sample stress as the horizontal axis data and the quantity of sample stress as the vertical axis data, and this system is labeled as the sample data coordinate system. Obtain the sample stress and the corresponding quantity as the coordinate points of the x and y coordinates, and mark them as sample data coordinate points; Plot all the coordinate points of the sample data on the sample data coordinate system to obtain a scatter plot, and label it as the sample data scatter plot.

3. The method for constructing a three-dimensional model of overburden load in a stope according to claim 2, characterized in that, Obtaining the first and second stress thresholds based on the scatter plot of sample data includes the following sub-steps: Get the maximum value of the ordinate among all sample data coordinate points and mark it as the sample data height; Set a length value and mark it as the rectangle length setting value; Create a rectangle on the horizontal axis of the sample data coordinate system with a height equal to the sample data height and a width equal to the set rectangle length, and be able to move left and right. Mark this rectangle as the set search rectangle. The total number of sample stresses corresponding to the coordinate points of the sample data included within the set search rectangle is marked as the rectangle search quantity; The range length of the sample stress is marked as the sample data length; Assuming the sample stress is uniformly distributed within the range of sample stress, the number of rectangle judgments at this time is obtained and marked as the average number of rectangles; Set a ratio, which is marked as the first set ratio; obtain the product of the average number of rectangles and the first set ratio, and mark it as the abnormal number threshold.

4. The method for constructing a three-dimensional model of overburden load in a stope according to claim 3, characterized in that, Obtaining the first and second stress thresholds based on the scatter plot of sample data also includes the following sub-steps: In the scatter plot of the sample data, the set search rectangle is shifted from the leftmost side to the right. When the number of rectangles searched is greater than or equal to the anomaly threshold, the movement of the set search rectangle is stopped. The sample stress corresponding to the minimum x-coordinate of the set search rectangle at this time is obtained and marked as the first stress threshold. In the scatter plot of the sample data, the set search rectangle is shifted from the rightmost side to the left. When the number of rectangles searched is greater than or equal to the anomaly threshold, the movement of the set search rectangle is stopped. The sample stress corresponding to the largest x-coordinate of the set search rectangle at this time is obtained and marked as the second stress threshold.

5. The method for constructing a three-dimensional model of overburden load in a stope according to claim 4, characterized in that, Obtaining the stress representation value from the first stress threshold and the second stress threshold includes the following sub-steps: The sample stress between the first stress threshold and the second stress threshold is marked as the screening stress; Obtain the mean value of the screening stress and mark it as the stress representation value; Obtain stress values ​​for the same region, the same stratum, and different overburden depths.

6. The method for constructing a three-dimensional model of overburden load in a stope according to claim 5, characterized in that, Obtaining the stress change function based on the stress representation value includes the following sub-steps: A plane rectangular coordinate system is established with overburden depth as the horizontal axis data and stress as the vertical axis data, and it is marked as the stress change coordinate system. Obtain the overburden depth and the corresponding stress values ​​on the x-axis and y-axis, and mark them as stress change coordinate points; Plot all stress variation coordinate points in the stress variation coordinate system; The function is obtained by fitting a function to all the stress change coordinate points, and it is labeled as the stress change function.

7. The method for constructing a three-dimensional model of overburden load in a stope according to claim 6, characterized in that, Obtaining the overlying rock to be detected and generating data points includes the following sub-steps: Obtain the vertically downward detection line from the ground detection point and mark it as the real-time detection line; Obtain stress detection points on the real-time detection line and mark them as real-time detection points; Obtain a three-dimensional rock stratum model of the overlying rock to be detected and mark it as a real-time three-dimensional rock stratum model; Obtain any point within the real-time 3D model of the same rock stratum and mark it as a generated data point.

8. The method for constructing a three-dimensional model of overburden load in a stope according to claim 7, characterized in that, Obtaining the first detection point based on the generated data points includes the following sub-steps: Get the nearest point on the real-time detection line with the same height as the generated data point and mark it as the line reference point; Obtain the real-time detection points adjacent to the line reference point on both sides and mark them as adjacent detection points; Obtain the stress values ​​of adjacent detection points and mark them as transformed stress; The adjacent detection points corresponding to the smaller transformation stress between the two transformation stresses are marked as the first detection point.

9. A method for constructing a three-dimensional model of overburden load in a stope according to claim 8, characterized in that, Obtaining the real-time input value based on the first detection point and the stress change function includes the following sub-steps: Obtain the points of transformed stress on the stress change function and mark them as function reference points; Obtain the stress change function between two function reference points and mark it as the function reference segment; Get the length of the function reference segment and mark it as the function reference length; Obtain the length of two adjacent detection points and mark it as the detection interval length; The detection interval length is scaled proportionally to match the function reference length, and the distance from the first detection point to the straight line reference point is obtained and marked as the first offset distance. Obtain the x-coordinate of the function reference point with the smaller x-coordinate, and mark it as the first starting distance; Get the sum of the first offset distance and the first starting distance, and mark it as the real-time input value.

10. A method for constructing a three-dimensional model of overburden load in a stope according to claim 9, characterized in that, The process of obtaining a function-fitted stress based on a stress variation function and real-time input values, and generating a three-dimensional model of the overburden load in the stope based on the function-fitted stress, includes the following sub-steps: Substitute the real-time input value as the horizontal axis value into the stress variation function to obtain the vertical value, which is then labeled as the function-fitted stress. The stress is fitted to the function to generate the stress at the data points; Obtain the stress of all generated data points within the real-time 3D rock strata model to generate a 3D model of overburden load in the mining area.

Images