Monitoring and evaluating system for geologic structure instability induced by geothermal drilling construction

Abstract

The invention discloses a monitoring and evaluating system for geologic structure instability induced by geothermal drilling construction, and relates to the technical field of geologic monitoring. The rock fracture average value of the multiple areas, the settlement values detected by the n settlement sensors and the crack values detected by the m displacement sensors are collected for monitoring, and y soil compression resistance values detected through a pressure testing method are obtained; the data evaluation and analysis module is used for evaluating and analyzing the data which does not reach the standard thresholds in the comparison process of the data values acquired by the data acquisition module and the standard thresholds to obtain the instability data of the rock fracture difference value, the topographic change value and the soil looseness difference value, and then generating the instability index, so that the situation that a large amount of data needs to be monitored is avoided; meanwhile, through comprehensive evaluation of multiple groups of data, it is guaranteed that the geological instability condition is accurately monitored and evaluated, and calculation for a large number of dip angles and geometric dimensions is avoided.

Description

Technical Field

[0001] The present invention relates to the field of geological monitoring technology, and in particular to a monitoring and evaluation system for geological structure instability induced by geothermal drilling construction. Background Art

[0002] Unfavorable geological bodies include isolated rocks with the risk of toppling, dangerous rock masses on steep slopes, and landslides. With the rapid development of large-scale infrastructure in recent years, civil engineering projects are often built near mountains, and more and more disasters of landslides caused by the instability of unfavorable geological bodies are encountered, resulting in serious loss of life and property. Therefore, the research and development of monitoring and early warning systems for unfavorable geological bodies is of great significance.

[0003] A poor geological body stability monitoring system and method disclosed in the patent application with reference publication number CN109520474B is used to quickly and accurately monitor the state of the poor geological body. A poor geological body stability monitoring method provided by the present invention includes: obtaining the inclination angle information of the poor geological body; transmitting the inclination angle information to the monitoring and early warning system; the inclination angle information includes the inclination angle and geometric dimensions of the poor geological body. The state of the poor geological body is determined by monitoring the inclination angle, which effectively improves the accuracy of the instability early warning of the poor geological body.

[0004] In the above scheme, the geological monitoring displacement should be mainly based on the change of inclination angle. During the monitoring process, single inclination angle monitoring is prone to errors. In addition, the existing technology focuses on the direct changes in the terrain of the surface process, and cannot accurately represent the instability changes in the deep geology that are hidden in the terrain and cannot be directly observed. It is impossible to accurately monitor and evaluate the geological instability, and the data is not intuitive enough for the calculation of a large number of inclination angles and geometric dimensions.

[0005] Therefore, it is necessary to provide a monitoring and evaluation system for geological structure instability induced by geothermal drilling construction to solve the above technical problems. Summary of the invention

[0006] The purpose of the present invention is to provide a monitoring and evaluation system for geological structure instability induced by geothermal drilling construction, so as to solve the problems of defects of the prior art mentioned in the above background technology.

[0007] Based on the above ideas, the present invention provides the following technical solutions: a monitoring and evaluation system for geological structure instability induced by geothermal drilling construction, the system comprising:

[0008] Data acquisition module, select area to divide into multiple Sub-areas, collect multiple The average rock fracture value in the region {PL zb1 ,PL zb2 ,...,PL zbw}, the settlement values ​​detected by n settlement sensors, the crack values ​​detected by m displacement sensors, and the y soil compressive strength values ​​detected by the pressure test method;

[0009] The data evaluation and analysis module compares the data values ​​collected by the data acquisition module with the standard thresholds, and evaluates and analyzes the data that does not reach the standard threshold to obtain the rock fracture difference CPL zb , Terrain change value DX bh and soil relaxation difference CTR sc The instability data is then used to generate the instability index SW zs ;

[0010] The data early warning module will generate the instability index SW through the data evaluation and analysis module zs and the preset instability index threshold SW zsd Compare and generate the instability difference SW zsc , and determine whether to generate an instability warning prompt. If an instability warning prompt is generated, determine the warning level;

[0011] The module is divided. When the data warning module determines whether to generate an instability warning prompt, a warning level is generated, and then measures to deal with the instability of the geological structure are implemented according to the generated warning level;

[0012] The database storage module is used to store the data running in the data collection module, the data evaluation and analysis module, the data early warning module and the division module.

[0013] As a further solution of the present invention: the area is divided into multiple The method of sub-region is as follows: by controlling the UAV to fly over the boundary of the selected area, the surveying image of the selected area is obtained, and the boundary of the surveying image is marked;

[0014] Mark multiple inherent walking paths in the surveying image;

[0015] Using multiple walking paths as dividing lines, the mapping image is divided into sub-image, each sub-image corresponds to Sub-area.

[0016] As a further solution of the present invention: the rock fracture difference CPL zb The acquisition methods include:

[0017] X-ray tomography In the sub-area, w rocks are randomly scanned with the same size standard;

[0018] Scanning the rock to obtain a scanned image;

[0019] Mark the pixel values ​​of all pixels in the rock scanning image and count the total number of pixels;

[0020] Compare the pixel values ​​of all pixels with the preset lower limit pixel value and the preset upper limit pixel value one by one;

[0021] The pixel points whose pixel values ​​are between the preset lower limit pixel value and the preset upper limit pixel value are recorded as fractured pixel points, and the number of fractured pixel points in the rock KX is counted d ;

[0022] The number of cracked pixels in the rock KX d The total number of pixels Z in rock scanning d Compared with the above, the percentage of broken pixels PL is obtained. zb ;

[0023] Rupture pixel ratio PL zb The expression is:

[0024]

[0025] Will After the w rock fracture proportion values ​​of the sub-areas are accumulated, the average rock fracture value ZPL is obtained zb ;

[0026] Rock fracture average ZPL zb The expression is:

[0027]

[0028] Where ZPL zb for The average rock fracture value of the sub-area, PL zba for The average value of the a-th rock fracture in the sub-area;

[0029] The rock fracture mean ZPL zb Compared with the standard rock fracture average BPL zb Compare and obtain the rock fracture difference CPL zb ;

[0030] Rock fracture difference CPL zb The expression is:

[0031] CPL zb =BPL zb -ZPL zb ;

[0032] As a further solution of the present invention: the terrain change value DX bh How to obtain:

[0033] exist n sensors are randomly arranged in each sub-area;

[0034] Through n sedimentation sensors, The settlement value in the sub-area and the m displacement sensors in time t The crack values ​​in the sub-areas are monitored;

[0035] The sedimentation values ​​of n sensors within t time are accumulated and averaged to obtain the average sedimentation value;

[0036] The expression for the mean sedimentation value is:

[0037]

[0038] In the formula, ZCJ bh for Sub-region subsidence change value, ZCJ b for The settlement value of the bth sensor in the sub-area at time t;

[0039] The crack values ​​of m sensors within t time are accumulated and averaged to obtain the crack average value ZLF bh ;

[0040] Crack average ZLF bh The expression is:

[0041]

[0042] In the formula, ZLF d is the crack value of the dth sensor in the sub-area after t time;

[0043] pass The average settlement value ZCJ in the sub-area bh Add the crack average ZLF bh , get The terrain change value DX in the sub-area bh ;

[0044] The terrain change value DX bh The expression is:

[0045] DX bh =ZLF bh +ZCJ bh .

[0046] As a further solution of the present invention: the soil relaxation difference CTR sc The methods for obtaining include:

[0047] By In the sub-area, dig at the same depth to obtain the same volume of y soils;

[0048] Detect the compression resistance of y soils by using the pressure test method;

[0049] Remove the highest and lowest values ​​of the y soil compressive strength values, and add up the average value to get the soil relaxation value TR sc ; Soil relaxation value TR sc The formula is:

[0050]

[0051] In the formula, TR kC For the The Cth soil compression value of the sub-area;

[0052] The soil relaxation value TR sc Compared with the standard soil relaxation value TR scb Compare and obtain the soil relaxation difference CTR SC ; Soil relaxation difference CTR sc The expression is:

[0053] CTR sc =TR scb -TR sc .

[0054] As a further solution of the present invention: the pressure testing method comprises:

[0055] By placing y pieces of soil in a container, setting a pressure sensor on a plate, and pressing the y pieces of soil down a distance x, the pressure sensor detects the pressure value when the soil is pressed down a further distance x, thereby detecting the compression resistance value of the y pieces of soil.

[0056] As a further solution of the present invention: the instability index SW zs The expression is:

[0057] SW zs =θ1×CPL zb +θ2×DX bh +θ3×CTR sc ;

[0058] Where SW zs is the geological instability index, θ1 is the rock fracture difference CPL zb The weight coefficient of θ2 is the terrain change value DX bh The weight coefficient of θ3 is the soil relaxation difference CTR sc ’s weight coefficient, and θ1+θ2+θ3=1.

[0059] As a further solution of the present invention: the instability difference SW zsc The expression is:

[0060] SW zsc =SW zsd -SW zs ;

[0061] Where SW zsd is the preset instability index threshold;

[0062] When the instability difference SW zsc When it is greater than or equal to 0, no geological instability warning will be given;

[0063] When the instability difference SW zsc When it is less than 0, a geological instability warning is issued.

[0064] As a further solution of the present invention: the warning level includes: a first warning level and a second warning level;

[0065] The generation method of the first and second warning levels is as follows:

[0066] The instability difference SW zsc With preset level threshold YS z Compare, YS z Less than 0;

[0067] When SW ZSC Greater than YS z If a first-level warning is generated, construction will be suspended and drilling parameters will be optimized;

[0068] When SW ZSC Less than or equal to YS z If a level 2 warning is generated, construction will be stopped immediately and personnel and equipment will be evacuated.

[0069] As a further solution of the present invention: the rock fracture average value, settlement value, soil compressive strength value and crack value collected by the data collection module are compared with the rock fracture average value threshold, settlement value threshold, soil compressive strength threshold and crack threshold in the historical instability data one by one:

[0070] If the value of one of the data items is greater than or equal to its item threshold, an early warning alarm is generated immediately;

[0071] If all data are less than the threshold values, the monitoring data will be sent to the data evaluation and analysis module for analysis and evaluation.

[0072] Compared with the prior art, the present invention has the following beneficial effects:

[0073] 1. Divide the selected area into multiple sub-areas, collect the average rock fracture values ​​of multiple areas, the settlement values ​​detected by n settlement sensors, the crack values ​​detected by m displacement sensors for monitoring, and the y soil compressive strength values ​​detected by the pressure test method, and obtain the instability data of rock fracture difference, terrain change value and soil relaxation difference through monitoring and analysis; then generate the instability index to avoid monitoring a large amount of data. At the same time, through the comprehensive evaluation of multiple groups of data, it is guaranteed to monitor and evaluate the geological instability accurately, avoiding the calculation of a large number of inclination angles and geometric dimensions. BRIEF DESCRIPTION OF THE DRAWINGS

[0074] The present invention is further described below in conjunction with the accompanying drawings and embodiments.

[0075] Figure 1 It is a schematic diagram of the overall structure of the present invention;

[0076] Figure 2 It is a schematic diagram of the process structure of the present invention. DETAILED DESCRIPTION

[0077] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0078] like Figure 1 to Figure 2 As shown, a monitoring and evaluation system for geological structure instability induced by geothermal drilling construction includes the following embodiments:

[0079] Embodiment 1: Figure 1 to Figure 2 As shown, the system includes:

[0080] Data acquisition module, select area to divide into multiple Sub-areas, collect multiple The average rock fracture value in the region {PL zb1 ,PL zb2 ,...,PL zbw}, the settlement values ​​detected by n settlement sensors, the crack values ​​detected by m displacement sensors, and the y soil compressive strength values ​​detected by the pressure test method;

[0081] The method for dividing a region into multiple sub-regions is specifically as follows: by controlling a drone to fly over the boundary of a selected region, a surveying image of the selected region is obtained by shooting, and the boundary of the surveying image is marked; multiple inherent walking paths are marked in the surveying image; the surveying image is divided into sub-images using the multiple walking paths as dividing lines, and each sub-image corresponds to a sub-region, respectively. By dividing the selected region into multiple small regions, data collection, analysis and processing can be facilitated;

[0082] The data evaluation and analysis module compares the data values ​​collected by the data acquisition module with the standard thresholds, and evaluates and analyzes the data that does not reach the standard threshold to obtain the rock fracture difference CPL zb , Terrain change value DX bh and soil relaxation difference CTR sc The instability data is then used to generate the instability index SW zs ;

[0083] The rock fracture average value, settlement value, soil compressive strength value and crack value collected by the data collection module are compared with the rock fracture average value threshold, settlement value threshold, soil compressive strength threshold and crack threshold in the historical instability data:

[0084] If the value of one of the data is greater than or equal to its threshold, an early warning alarm is generated immediately; if the values ​​of all the data are less than the thresholds, the monitoring data is sent to the data evaluation and analysis module for analysis and evaluation. In the actual monitoring process, when any of the average rock fracture value, settlement value, compressive strength value and crack value reaches the preset threshold in the setting, it means that the geological instability has become irreversible. If the monitored data temporarily does not reach the threshold, an early warning cannot be issued. However, it is impossible to make a timely judgment on the instability situation through comprehensive data. Therefore, the accuracy of monitoring is improved by further evaluating and analyzing the average rock fracture value, settlement value, compressive strength value and crack value that do not reach the threshold.

[0085] Example 2: The rock fracture difference CPL zb The acquisition methods include:

[0086] X-ray tomography In the sub-area, w rocks are randomly scanned with the same size standard;

[0087] Scanning the rock to obtain a scanned image;

[0088] The specific operations are as follows

[0089] Conduct macroscopic observation of rock samples and select representative core samples; determine the sample size according to the scanning resolution requirements to ensure that the sample does not exceed the detector receiving range during the 360° scanning process; turn on the instrument power, preheat the X-ray source and detector, and set relevant parameters; place the sample on the sample stage and fix it, adjust the sample position so that the sample image is within the detector field of view; adjust the X-ray source and detector position according to the two-dimensional preview image of the sample on the detector to ensure that the scanned image does not exceed the detector field of view;

[0090] Mark the pixel values ​​of all pixels in the rock scan image and count the total number of pixels; compare the pixel values ​​of all pixels with the preset lower limit pixel value and the preset upper limit pixel value one by one; record the pixel points whose pixel values ​​are between the preset lower limit pixel value and the preset upper limit pixel value as fractured pixels, and count the number of fractured pixels in the rock KX d ; The number of cracked pixels in the rock KX d The total number of pixels Z in rock scanning d Compared with the above, the percentage of broken pixels PL is obtained. zb ;

[0091] Rupture pixel ratio PL zb The expression is:

[0092]

[0093] Will After the w rock fracture proportion values ​​of the sub-areas are accumulated, the average rock fracture value ZPL is obtained zb ;

[0094] Rock fracture average ZPL zb The expression is:

[0095]

[0096] Where ZPL zb for The average rock fracture value of the sub-area, PL zba for The average rock fracture value of the ath sub-area; the average rock fracture value ZPL zb Compared with the standard rock fracture average BPL zb Compare and obtain the rock fracture difference CPL zb ,

[0097] Rock fracture difference CPL zb The expression is:

[0098] CPL zb =BPL zb -ZPL zb ;

[0099] When the rock breaks, the average ZPL zb The larger the value, the The rock fracture conditions in the sub-area are used to determine the texture change, and the rock fracture difference CPL is used to determine the texture change. zb Intuitively show the geological changes, rock fracture difference CPL zb The larger it is, the greater the geological changes are and the higher the possibility of geological instability is.

[0100] The terrain change value DX bh How to obtain:

[0101] exist n sensors are randomly arranged in each sub-area;

[0102] Through n sedimentation sensors, The settlement value in the sub-area and the m displacement sensors in time t The crack values ​​in the sub-areas are monitored;

[0103] The sedimentation values ​​of n sensors within t time are accumulated and averaged to obtain the average sedimentation value;

[0104] The expression for the mean sedimentation value is:

[0105]

[0106] In the formula, ZCJ bh for Sub-region subsidence change value, ZCJ b for The settlement value of the bth sensor in the sub-area at time t;

[0107] The crack values ​​of m sensors within t time are accumulated and averaged to obtain the crack average value ZLF bh ;

[0108] Crack average ZLF bh The expression is:

[0109]

[0110] In the formula, ZLF d is the crack value of the dth sensor in the sub-area after t time;

[0111] pass The average settlement value ZCJ in the sub-area bh Add the crack average ZLF bh , get The terrain change value DX in the sub-area bh ;

[0112] The terrain change value DX bh The expression is:

[0113] DX bh =ZLF bh +ZCJ bh .

[0114] The soil relaxation difference CTR sc The methods for obtaining include:

[0115] By In the sub-area, dig at the same depth to obtain the same volume of y soils;

[0116] Detect the compression resistance of y soils by using the pressure test method;

[0117] The pressure testing method comprises:

[0118] By placing y pieces of soil in a container, setting a pressure sensor on a plate, and pressing the y pieces of soil down by a distance x, the pressure sensor detects the pressure value when the soil is pressed down by the distance x again, thereby detecting the compression resistance value of the y pieces of soil;

[0119] Remove the highest and lowest values ​​of the y soil compressive strength values ​​to avoid inaccurate data caused by excessive data deviation, and add up the average value to obtain the soil relaxation value TR sc ;

[0120] Soil relaxation value TR sc The formula is:

[0121]

[0122] In the formula, TR kC For the The Cth soil compression value of the sub-area;

[0123] The soil relaxation value TR sc Compared with the standard soil relaxation value TR scb Compare and obtain the soil relaxation difference CTR sc ; Soil relaxation difference CTR sc The expression is:

[0124] CTR sc =TR scb -TR sc .

[0125] As a further solution of the present invention: the instability index SW zs The expression is:

[0126] SWzs =θ1×CPL zb +θ2×DX bh +θ3×CTR sc ;

[0127] Where SW zs is the geological instability index, θ1 is the rock fracture difference CPL zb The weight coefficient of θ2 is the terrain change value DX bh The weight coefficient of θ3 is the soil relaxation difference CTR sc ’s weight coefficient, and θ1+θ2+θ3=1.

[0128] It should be noted that the size of the weight factor is a specific value obtained by quantifying each data to facilitate subsequent comparison. The size of the weight factor depends on the amount of unstable data and the preliminary setting of the corresponding weight factor for each set of comprehensive pest data by technical personnel in this field.

[0129] The data early warning module will generate the instability index SW through the data evaluation and analysis module zs and the preset instability index threshold SW zsd Compare and generate the instability difference SW zsc , and determine whether to generate an instability warning prompt. If an instability warning prompt is generated, determine the warning level;

[0130] The module is divided. When the data warning module determines whether to generate an instability warning prompt, a warning level is generated, and then measures to deal with the instability of the geological structure are implemented according to the generated warning level;

[0131] The database storage module is used to store the data running in the data collection module, the data evaluation and analysis module, the data early warning module and the division module.

[0132] The instability difference SW zsc The expression is:

[0133] SW zsc =SW zsd -SW zs ;

[0134] Where SW zsd is the preset instability index threshold;

[0135] When the instability difference SW zsc When it is greater than or equal to 0, no geological instability warning will be given;

[0136] When the instability difference SW zsc When it is less than 0, a geological instability warning is issued.

[0137] In this scheme, the selected area is divided into multiple sub-areas, and the average rock fracture values ​​of multiple areas, the settlement values ​​detected by n settlement sensors, the crack values ​​detected by m displacement sensors, and the y soil compressive strength values ​​detected by the pressure test method are collected for monitoring. The rock fracture difference CPL is obtained through monitoring analysis. zb , Terrain change value DX bh and soil relaxation difference CTR sc The instability data is then generated by the instability index SW zs , avoiding the monitoring of a large amount of data. Through the comprehensive evaluation of multiple sets of data, accurate monitoring and evaluation of geological instability can be ensured, avoiding the calculation of a large number of inclination angles and geometric dimensions.

[0138] The warning levels include: level one warning level and level two warning level;

[0139] The generation method of the first and second warning levels is as follows:

[0140] The instability difference SW zsc With preset level threshold YS z Compare, YS z Less than 0; instability difference SW zsc The larger the value, the farther away from the unstable state. zsc The smaller the value, the closer it is to the unstable state. zsc When it is less than 0, it means that the state of instability warning production has been reached.

[0141] When the instability difference SW zsc Greater than YS z If a first-level warning is generated, construction will be suspended and drilling parameters will be optimized;

[0142] When the instability difference SW zsc Less than or equal to YS z If a level 2 warning is generated, construction will be stopped immediately and personnel and equipment will be evacuated.

[0143] Meanwhile, the contents not described in detail in this specification belong to the prior art known to those skilled in the art.

[0144] It should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device.

[0145] Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the appended claims and their equivalents.

Claims

1. A monitoring and evaluation system for geological structural instability induced by geothermal drilling construction, characterized in that the system include: Data acquisition module, select area to divide into multiple Sub-areas, collect multiple The average rock fracture value in the region {PL zb1 ,PL zb2 ,...,PL zbw }, the settlement values ​​detected by n settlement sensors, the crack values ​​detected by m displacement sensors, and the y soil compressive strength values ​​detected by the pressure test method; The data evaluation and analysis module compares the data values ​​collected by the data acquisition module with the standard thresholds, and evaluates and analyzes the data that does not reach the standard threshold to obtain the rock fracture difference CPL zb , Terrain change value DX bh and soil relaxation difference CTR sc The instability data is then used to generate the instability index SW zs ; The data early warning module will generate the instability index SW through the data evaluation and analysis module zs and the preset instability index threshold SW zsd Compare and generate the instability difference SW zsc , and determine whether to generate an instability warning prompt. If an instability warning prompt is generated, determine the warning level; The module is divided. When the data warning module determines whether to generate an instability warning prompt, a warning level is generated, and then measures to deal with the instability of the geological structure are implemented according to the generated warning level; The database storage module is used to store the data running in the data collection module, the data evaluation and analysis module, the data early warning module and the division module.

2. A monitoring and evaluation system for geological structure instability induced by geothermal drilling construction according to claim 1, characterized in that: The area is divided into multiple The method of sub-region is as follows: by controlling the UAV to fly over the boundary of the selected area, the surveying image of the selected area is obtained, and the boundary of the surveying image is marked; Mark multiple inherent walking paths in the surveying image; Using multiple walking paths as dividing lines, the mapping image is divided into sub-image, each sub-image corresponds to Sub-area.

3. A monitoring and evaluation system for geological structure instability induced by geothermal drilling construction according to claim 2, characterized in that: The rock failure difference CPL ab The acquisition methods include: X-ray tomography In the sub-area, w rocks are randomly scanned with the same size standard; Scanning the rock to obtain a scanned image; Mark the pixel values ​​of all pixels in the rock scanning image and count the total number of pixels; Compare the pixel values ​​of all pixels with the preset lower limit pixel value and the preset upper limit pixel value one by one; The pixel points whose pixel values ​​are between the preset lower limit pixel value and the preset upper limit pixel value are recorded as fractured pixel points, and the number of fractured pixel points in the rock KX is counted d ; The number of cracked pixels in the rock KX d The total number of pixels Z in rock scanning d Compared with the above, the percentage of broken pixels PL is obtained. zb ; Rupture pixel ratio PL zb The expression is: Will After the w rock fracture proportion values ​​of the sub-areas are accumulated, the average rock fracture value ZPL is obtained zb ; Rock fracture average ZPL zb The expression is: Where ZPL zb for The average rock fracture value of the sub-area, PL zba for The average value of the a-th rock fracture in the sub-area; The rock fracture mean ZPL zb Compared with the standard rock fracture average BPL zb Compare and obtain the rock fracture difference CPL zb ; Rock fracture difference CPL zb The expression is: CPL zb =BPL zb -ZPL zb 。 4. A monitoring and evaluation system for geological structure instability induced by geothermal drilling construction according to claim 3, characterized in that: The terrain change value DX bh How to obtain: exist n sensors are randomly arranged in each sub-area; Through n sedimentation sensors, The settlement value in the sub-area and the m displacement sensors in time t The crack values ​​in the sub-areas are monitored; The sedimentation values ​​of n sensors within t time are accumulated and averaged to obtain the average sedimentation value; The expression for the mean sedimentation value is: In the formula, ZCJ bh for Sub-region subsidence change value, ZCJ b for The settlement value of the bth sensor in the sub-area at time t; The crack values ​​of m sensors within t time are accumulated and averaged to obtain the crack average value ZLF bh ; Crack average ZLF bh The expression is: In the formula, ZLF d is the crack value of the dth sensor in the sub-area after t time; pass The average settlement value ZCJ in the sub-area bh Add the crack average ZLF bh , get The terrain change value DX in the sub-area bh ; The terrain change value DX bh The expression is: DX bh =ZLF bh +ZCJ bh 。 5. A monitoring and evaluation system for geological structure instability induced by geothermal drilling construction according to claim 4, characterized in that: The soil relaxation difference CTR sc The methods for obtaining include: By In the sub-area, dig at the same depth to obtain the same volume of y soils; Detect the compression resistance of y soils by using the pressure test method; Remove the highest and lowest values ​​of the y soil compressive strength values, and add up the average value to get the soil relaxation value TR sc ; Soil relaxation value TR sc The formula is: In the formula, TR kC For the The Cth soil compression value of the sub-area; The soil relaxation value TR sc Compared with the standard soil relaxation value TR scb Compare and obtain the soil relaxation difference CTR SC ; Soil relaxation difference CTR sc The expression is: CTR sc =TR scb -TR sc 。 6. A monitoring and evaluation system for geological structure instability induced by geothermal drilling construction according to claim 5, characterized in that: The pressure testing method comprises: By placing y pieces of soil in a container, setting a pressure sensor on a plate, and pressing the y pieces of soil down a distance x, the pressure sensor detects the pressure value when the soil is pressed down a further distance x, thereby detecting the compression resistance value of the y pieces of soil.

7. A monitoring and evaluation system for geological structure instability induced by geothermal drilling construction according to claim 6, characterized in that: The instability index SW zs The expression is: SW zs =θ1×CPL zb +θ2×DX bh +θ3×CTR sc ; Where SW zs is the geological instability index, θ1 is the rock fracture difference CPL zb The weight coefficient, θ2 is the terrain change value DX bh The weight coefficient of θ3 is the soil relaxation difference CTR sc ’s weight coefficient, and θ1+θ2+θ3=1.

8. A monitoring and evaluation system for geological structure instability induced by geothermal drilling construction according to claim 7, characterized in that: The instability difference SW zsc The expression is: SW zsc =SW zsd -SW zs ; Where SW zsd is the preset instability index threshold; When the instability difference SW zsc When it is greater than or equal to 0, no geological instability warning will be given; When the instability difference SW zsc When it is less than 0, a geological instability warning is issued.

9. A monitoring and evaluation system for geological structure instability induced by geothermal drilling construction according to claim 8, characterized in that: The warning levels include: level one warning level and level two warning level; The generation method of the first and second warning levels is as follows: The instability difference SW zsc With preset level threshold YS z Compare, YS z Less than 0; When the instability difference SW zsc Greater than YS z If a first-level warning is generated, construction will be suspended and drilling parameters will be optimized; When the instability difference SW zsc Less than or equal to YS z If a level 2 warning is generated, construction will be stopped immediately and personnel and equipment will be evacuated.

10. A monitoring and evaluation system for geological structure instability induced by geothermal drilling construction according to claim 1, characterized in that: The rock fracture average value, settlement value, soil compressive strength value and crack value collected by the data collection module are compared with the rock fracture average value threshold, settlement value threshold, soil compressive strength threshold and crack threshold in the historical instability data: If the value of one of the data items is greater than or equal to its item threshold, an early warning alarm is generated immediately; If all data are less than the threshold values, the monitoring data will be sent to the data evaluation module for analysis.

Images