Non-destructive assessment method and device for the instability risk level of cave temples

Abstract

The present invention provides a nondestructive assessment method and device for the instability risk level of a grotto temple. The method comprises: determining a risk level standard corresponding to each evaluation indicator, the risk level standard including a value range corresponding to each risk level; determining an actual indicator value corresponding to each evaluation indicator of the grotto temple to be evaluated, normalizing each actual indicator value and each value range to obtain each normalized indicator value and each normalized value range; determining a subjective weight corresponding to each evaluation indicator based on the analytic hierarchy process, determining an objective weight corresponding to each evaluation indicator using the entropy weight method, and determining a comprehensive weight corresponding to each evaluation indicator based on the subjective and objective weights corresponding to each evaluation indicator; and determining a risk level assessment result corresponding to the grotto temple to be evaluated using a distance discriminant method based on each comprehensive weight, each normalized indicator value, and the normalized value range corresponding to each risk level. This application improves the accuracy of the assessment results of the instability risk of grotto temples.

Description

Technical Field

[0001] The present invention relates to the technical field of cave temple risk assessment, and in particular to a non-destructive assessment method and device for the instability risk level of a cave temple. Background Art

[0002] Grotto temples, a splendid legacy of Chinese culture, have not only witnessed the passage of history but also carry profound cultural and artistic value. However, under the long-term influence of natural forces and human activities, the grotto rock masses are prone to deterioration, such as cracking, pulverization, and flaking, which seriously threaten their long-term preservation. To more effectively address these challenges, in-depth and systematic conservation research is urgently needed, of which scientific assessment of the instability risk of grotto rock masses is particularly critical. Considering the diversity of grotto lithology, occurrence environments, and cave structures, as well as the variety of deterioration and instability characteristics, developing a method to assess the instability risk of grotto temples has significant theoretical and practical significance for guiding the development of targeted conservation measures and extending the life cycle of these cultural treasures.

[0003] At present, the evaluation system for grotto stability is relatively simple, mainly focusing on the grotto structure itself, and lacking more comprehensive considerations, which limits the accuracy and reliability of the evaluation results. my country has a vast territory, and the engineering geological environment in different regions varies significantly. This difference is the root cause of the regional characteristics of the deterioration and instability of grotto temples. The existing evaluation methods do not fully consider the comprehensive impact of key factors such as active structures, climate, hydrology and landforms on grotto temples, resulting in inaccurate evaluation results of grotto temple stability.

[0004] In summary, while existing cave temple instability risk assessment methods can reflect the stability of cave temples to a certain extent, they have limitations in terms of accuracy. Therefore, how to improve the accuracy of cave temple instability risk assessment results is a technical problem that needs to be solved urgently. Summary of the Invention

[0005] In view of this, an embodiment of the present invention provides a non-destructive assessment method and device for the instability risk level of a cave temple, so as to eliminate or improve one or more defects existing in the prior art.

[0006] One aspect of the present invention provides a non-destructive assessment method for the instability risk level of a cave temple, the method comprising:

[0007] Determine the risk level standards corresponding to each evaluation indicator, including rock durability, occurrence environment, cave structure stability, rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock mass integrity and degree of disease development. The risk level standards include the value range corresponding to each level of risk;

[0008] Determine the actual indicator value corresponding to each evaluation indicator of the grotto temple to be evaluated, and normalize the actual indicator value and each value interval to obtain each normalized indicator value and each normalized value interval;

[0009] Determine the subjective weight corresponding to each evaluation indicator based on the hierarchical analysis method, determine the objective weight corresponding to each evaluation indicator based on the normalized indicator value and each normalized value interval through the entropy weight method, and determine the comprehensive weight corresponding to each evaluation indicator based on the subjective weight and objective weight corresponding to each evaluation indicator;

[0010] Based on the comprehensive weights, the normalized index values and the normalized value intervals corresponding to the risks at each level, the risk level assessment result corresponding to the grotto temple to be evaluated is determined by the distance discrimination method.

[0011] In some embodiments of the present invention, determining the actual indicator value corresponding to each evaluation indicator of the grotto temple to be evaluated includes:

[0012] Determine the actual rock type, actual occurrence environment, and actual cave shape of the grotto temple to be evaluated, and determine the actual index values of the rock durability, occurrence environment, and cave structure stability evaluation indicators of the grotto temple to be evaluated based on the actual rock type, actual occurrence environment, and actual cave shape and the risk level standard;

[0013] Conduct on-site testing on the grotto temple to be evaluated to obtain actual index values corresponding to the rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock integrity and disease development degree evaluation indicators.

[0014] In some embodiments of the present invention, determining the subjective weight corresponding to each evaluation indicator based on the analytic hierarchy process includes:

[0015] Determine the relative importance of each evaluation indicator, and construct a judgment matrix based on the relative importance of each evaluation indicator;

[0016] If the judgment matrix meets the consistency test requirements, determining the maximum eigenvalue in the judgment matrix, and calculating the eigenvector corresponding to the maximum eigenvalue;

[0017] The subjective weight corresponding to each evaluation indicator is calculated based on the feature vector.

[0018] In some embodiments of the present invention, the objective weight corresponding to each evaluation indicator is determined by an entropy weight method based on each normalized indicator value and each normalized value interval, including:

[0019] Determine the upper limit of the normalized value interval corresponding to the risk level to which each normalized indicator value belongs;

[0020] Calculate the contribution degree of each evaluation indicator based on the value of each normalized indicator and the upper limit of each normalized value interval;

[0021] Calculating the entropy corresponding to each evaluation indicator based on the contribution degree corresponding to each evaluation indicator;

[0022] The objective weight corresponding to each evaluation indicator is calculated based on the entropy corresponding to each evaluation indicator.

[0023] In some embodiments of the present invention, the calculation formula of the comprehensive weight is:

[0024] w j =ξw j1 +(1-ξ)w j2 ;

[0025] Among them, w j1 represents the subjective weight of evaluation index j, w j2 represents the objective weight of the evaluation index j, ξ = 0.5.

[0026] In some embodiments of the present invention, the entropy is calculated as follows:

[0027]

[0028] Among them, f ij represents the contribution of the jth evaluation indicator in the i-th grotto temple to be evaluated, n represents the total number of grotto temples to be evaluated, V represents the number of risk levels, and V takes values of Ⅰ, Ⅱ, Ⅲ, Ⅳ or Ⅴ, f kj Indicates the contribution of the jth evaluation indicator in the kth level risk, C' ij represents the normalized index value corresponding to the jth evaluation index in the i-th grotto temple to be evaluated, B' kj It represents the upper limit of the normalized value range of the jth evaluation indicator in the kth level risk.

[0029] In some embodiments of the present invention, determining the risk level standard corresponding to each evaluation indicator includes:

[0030] Based on the assignment method, the rock durability, occurrence environment, and cave structure stability evaluation indicators are divided into five levels respectively;

[0031] Based on the empirical method, the rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock integrity and disease development degree evaluation indicators are divided into five levels respectively.

[0032] In some embodiments of the present invention, the risk level assessment result corresponding to the grotto temple to be evaluated is determined by a distance discrimination method based on the comprehensive weights, the normalized index values, and the normalized value intervals corresponding to the risk levels, including:

[0033] Based on the comprehensive weights, the normalized index values and the upper limits of the normalized value intervals, the two-dimensional Euclidean distance values between the grotto temple to be evaluated and each risk level are calculated;

[0034] Based on the two-dimensional Euclidean distance value, the risk level assessment result corresponding to the grotto temple to be evaluated is determined by the clamping criterion.

[0035] According to another aspect of the present invention, a non-destructive assessment system for the instability risk level of a cave temple is disclosed. The system includes a processor, a memory, and a computer program stored in the memory. The processor is used to execute the computer program. When the computer program is executed, the system implements the steps of the method described in any of the above embodiments.

[0036] According to yet another aspect of the present invention, a computer-readable storage medium is disclosed, on which a computer program is stored. When the computer program is executed by a processor, the steps of the method described in any of the above embodiments are implemented.

[0037] The non-destructive assessment method and device for the instability risk level of grotto temples disclosed in the above-mentioned embodiments of the present invention include evaluation indicators such as rock durability, occurrence environment, cave structure stability, rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock integrity and degree of disease development. It not only covers the main characteristics of the grotto temples, but also deeply considers the deterioration state of rocks and rock structures, as well as the occurrence environment in which the grotto temples are located. This multi-dimensional evaluation system can more accurately reflect the actual situation of the stability of the grotto temples, provide solid data support for the classification of risk levels, and thus improve the accuracy of the assessment results of the instability risk of the grotto temples.

[0038] In addition, this application adopts a comprehensive weighting method to determine the weight of each evaluation indicator, making the weight distribution more scientific and reasonable, avoiding the one-sidedness and arbitrariness of subjective weighting; at the same time, combined with the distance judgment method and the squeeze criterion, the risk level evaluation results are more accurate and reliable.

[0039] Additional advantages, objects, and features of the present invention will be set forth in part in the following description and will become apparent to those skilled in the art upon examination of the following or may be learned from practice of the present invention. The objects and other advantages of the present invention may be realized and obtained by the structures particularly pointed out in the description and drawings.

[0040] Those skilled in the art will understand that the purposes and advantages that can be achieved by the present invention are not limited to the above specific descriptions, and the above and other purposes that can be achieved by the present invention will be more clearly understood based on the following detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] The drawings described herein are intended to provide a further understanding of the present invention, constitute a part of this application, and do not constitute a limitation of the present invention. The components in the drawings are not drawn to scale, but are merely for the purpose of illustrating the principles of the present invention. To facilitate the illustration and description of certain portions of the present invention, corresponding portions in the drawings may be exaggerated, that is, may be larger than other components in an exemplary device actually manufactured according to the present invention. In the drawings:

[0042] Figure 1 The figure is a flow chart of a non-destructive assessment method for the instability risk level of a cave temple according to an embodiment of the present invention.

[0043] Figure 2 Schematic diagram of a flow chart of a non-destructive assessment method for the instability risk level of a cave temple according to another embodiment of the present invention. DETAILED DESCRIPTION

[0044] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the embodiments and the accompanying drawings. Here, the exemplary embodiments of the present invention and their descriptions are used to explain the present invention, but are not intended to limit the present invention.

[0045] It should also be noted that, in order to avoid obscuring the present invention due to unnecessary details, the accompanying drawings only show structures and / or processing steps closely related to the solutions according to the present invention, while other details that are not closely related to the present invention are omitted.

[0046] It should be emphasized that the term "include / comprises" when used herein refers to the existence of features, elements, steps or components, but does not exclude the existence or addition of one or more other features, elements, steps or components.

[0047] It should also be noted here that, unless otherwise specified, the term "connection" in this article can refer not only to a direct connection, but also to an indirect connection with an intermediary, and not only to a wired connection but also to a wireless connection, and the specific connection can be changed based on the actual application scenario.

[0048] Hereinafter, embodiments of the present invention will be described with reference to the accompanying drawings. In the accompanying drawings, the same reference numerals represent the same or similar components, or the same or similar steps.

[0049] Figure 1 FIG. 1 is a flow chart of a non-destructive assessment method for the instability risk level of a cave temple according to an embodiment of the present invention. Figure 1 As shown, the non-destructive assessment method for the instability risk level of a cave temple includes at least steps S10 to S40.

[0050] Step S10: Determine the risk level standards corresponding to each evaluation indicator, where the evaluation indicators include rock durability, occurrence environment, cave structure stability, rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock integrity and degree of disease development. The risk level standards include the value range corresponding to each level of risk.

[0051] In this step, a risk level standard corresponding to each evaluation indicator is constructed as a reference system for evaluating the risk level of grotto temples. The rock longitudinal wave velocity, rock heat absorption coefficient, and rock Leeb hardness evaluation indicators can also be collectively referred to as rock physical property evaluation indicators. Specifically, multiple evaluation indicators can be assigned the same number of risk levels. In other words, the risk level standard represents the corresponding relationship between each level of risk and the numerical value range of each evaluation indicator.

[0052] In some embodiments, rock durability mainly refers to the durability of sandstone, conglomerate, mudstone, limestone or crystalline rock, so the rock durability index can be graded based on the type of rock; the hosting environment includes climate, landform, hydrology and active structural characteristics, so the risk level corresponding to the hosting environment index can be graded based on the specific scores of climate, landform, hydrology and active structure; the stability of cave structure refers to the stability of flat-top caves, arched caves, gabled-top caves, flat-top central column caves or bucket-top caves, so the cave structure stability index can be graded based on the cave shape of the grotto temple; in addition, the longitudinal wave velocity, heat absorption coefficient and Leeb hardness in the rock physical properties can be graded based on the scores respectively; the rock integrity is characterized by the integrity coefficient Kv, so the rock integrity index can be graded based on the Kv value; the degree of disease development refers to the proportion of the development area of water erosion diseases, weathering diseases and biological diseases to the cave area, so the disease development degree index can be graded based on the proportion value.

[0053] Exemplarily, each evaluation indicator is divided into five levels: I, II, III, IV, and V, with the lowest risk of instability of the cave temple at level I and the highest risk at level V. Exemplarily, the total number of risk levels is five, and determining the risk level standard corresponding to each evaluation indicator may include the following steps: based on the assignment method, the rock durability, occurrence environment, and cave structure stability evaluation indicators are divided into five levels respectively; based on the empirical method, the rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock mass integrity, and disease development degree evaluation indicators are divided into five levels respectively.

[0054] Specifically, the expert assignment method (10-point scale) can be used to divide the rock durability index into five levels, among which crystalline rock, limestone, sandstone, conglomerate and mudstone are classified as levels I, II, III, IV and V respectively, and the corresponding value intervals of each level are: 10 to 8, 8 to 6, 6 to 4, 4 to 2 and 2 to 0 respectively. Similarly, the expert assignment method (10-point scale) is used to divide the cave structure stability index into five levels, among which flat-top central pillar cave, bucket-top cave, gable-top cave, arch cave and flat-top cave are classified as levels I, II, III, IV and V respectively, and the corresponding value intervals of each level are: 10 to 8, 8 to 6, 6 to 4, 4 to 2 and 2 to 0 respectively. The expert assignment method (10-point scale) is used to score the climate, landform, hydrology and active structure in the cave temple environment, and the values are recorded as F q 、F d 、F s and F h , then the comprehensive score of the environment is given by F q 、F d 、F s and F h Jointly determined, that is, comprehensive score = 0.25F q +0.25F d +0.25F s +0.25F h ; Based on this, according to expert opinions, the occurrence environment indicators are divided into five levels: I, II, III, IV, and V, and the corresponding value intervals for each level are: 10 to 8, 8 to 6, 6 to 4, 4 to 2, and 2 to 0. In addition, the grading standards for the integrity coefficient, longitudinal wave velocity, heat absorption coefficient, Leeb hardness, and degree of disease development can be determined based on expert experience. Similar to the above, each indicator is divided into levels I, II, III, IV, and V, and each level has a corresponding value interval. The following table is a grading reference table for the risk level standards of an embodiment of the present application:

[0055]

[0056] Step S20: determining the actual indicator values corresponding to the evaluation indicators of the grotto temple to be evaluated, and normalizing the actual indicator values and the value intervals to obtain normalized indicator values and normalized value intervals.

[0057] In this step, the actual indicator values of the grotto temple to be evaluated are obtained. Since the grading in step S10 is generally based on a 10-point scale, in order to better represent the evaluation indicators, the actual indicator values and the value intervals in the grading reference table are further normalized.

[0058] In one embodiment, determining the actual indicator values corresponding to each evaluation indicator of the grotto temple to be evaluated may specifically include: determining the actual rock type, actual occurrence environment, and actual cave shape of the grotto temple to be evaluated; and determining the actual indicator values of the rock durability, occurrence environment, and cave structure stability evaluation indicators of the grotto temple to be evaluated based on the actual rock type, actual occurrence environment, and actual cave shape and the risk level standard; and conducting on-site testing on the grotto temple to be evaluated to obtain the actual indicator values corresponding to the rock longitudinal wave velocity, rock heat absorption coefficient, rock hardness, rock integrity, and disease development degree evaluation indicators. In this embodiment, the rock type, occurrence environment, and cave shape of the grotto temple to be evaluated can be ascertained through on-site investigations, etc., and then, based on the investigation results, the risk level corresponding to the rock type, occurrence environment, and cave shape indicators of the grotto temple to be evaluated can be pre-defined. Then, based on the risk level corresponding to each indicator, a corresponding value range can be determined from a graded reference table. Finally, based on experience, experts select an appropriate value from the corresponding value range as the actual indicator value of the corresponding evaluation indicator. In addition, based on field tests, the rock physical properties, rock integrity and degree of disease development of the grotto temple to be evaluated are obtained, and the actual indicator values corresponding to each indicator are determined.

[0059] For example, the wave velocity, heat absorption coefficient and Leeb hardness of rock physical properties can be measured based on RSM-SY6(C) non-metallic ultrasonic detector, Fluke Ti480 Pro infrared thermal imager and TIME5350 Leeb hardness tester respectively; while the rock integrity coefficient Kv is obtained by the longitudinal wave velocity V of the rock mass. m and the longitudinal wave velocity V of the rock mass r Determined by the square of the ratio, that is, Kv=V m 2 / V r 2 ; P-wave velocity V m and the longitudinal wave velocity V r These are the rock mass and rock mass wave velocities of the grotto temple being evaluated, measured using the RSM-SY6(C) non-metallic ultrasonic detector. Furthermore, the extent of water erosion, weathering, and biological damage can be determined through on-site statistical measurements.

[0060] Specifically, when using the RSM-SY6(C) non-metallic ultrasonic detector to measure the wave velocity of the grotto temple to be evaluated, in order to ensure the reliability of the test results, a conical sensor probe is preferred, and the distance between the two probes is set to 25mm; before the test, the two probes are docked to measure the instrument delay, and this delay time is used as the zero time in the instrument parameter setting interface. Wave velocity measurement includes rock block wave velocity measurement and rock mass wave velocity measurement; when measuring the rock block wave velocity V rWhen measuring the rock mass wave velocity V, it is necessary to collect cylindrical sample rock blocks with a diameter of 50 mm and a height of 25 mm from the cave research area in advance and use a 50KHz plane probe to measure the wave velocity. m When measuring, it is necessary to select a cone-shaped sensor probe to measure directly on the rock surface. In addition, after measuring the rock block wave velocity and rock mass wave velocity of the grotto temple to be evaluated, the formula Kv=V m 2 / V r 2 The rock mass integrity coefficient Kv is calculated.

[0061] In addition, when obtaining the actual index value of the rock heat absorption coefficient of the grotto temple to be evaluated, a constant heat flux halogen heating lamp with adjustable power is selected to thermally excite the grotto rock mass, and a Fluke Ti480 Pro portable infrared thermal imager and a radiation heat flux meter are selected to obtain the temperature cloud map and surface heat flux density q during the rock mass heating process; the SmartView infrared image is used to extract the heating process curve and initial temperature T0 of the complete area; the heat flux density q, initial temperature T0 and time t of the complete area are substituted into Calculate the average heat absorption coefficient during the heating process ierfc(0) is the first integral of the Gaussian error complement function, with a value of 0.5642, T(0,t) refers to the temperature of the rock surface at time t, λ refers to the thermal coefficient, ρ refers to the density, and c refers to the specific heat capacity.

[0062] To obtain the actual Leeb hardness values for the grottoes to be evaluated, a TIME5350 Leeb hardness tester, preferably with a D-type universal probe, was used. This instrument is compact, lightweight, and simple to test, with minimal impact on the rock mass during testing. Furthermore, to obtain the actual values for the degree of damage to the grottoes to be evaluated, manual on-site measurements can be used to determine the extent of water erosion, weathering, and biological damage.

[0063] After obtaining the actual indicator values of the grotto temples to be evaluated, the actual indicator values and the value intervals corresponding to the risks at each level in the above-mentioned graded reference table are further normalized. Assuming that risk assessment is conducted on multiple grotto temples to be evaluated, the j-th indicator data of the i-th (i=0,…,n) grotto temple is recorded as C ij The upper limit of the jth indicator of the kth risk level (k=Ⅰ,…,V) is recorded as B kj That is, B kj Indicates the upper limit of the value interval corresponding to the kth level of the jth indicator, the normalized C ij and B kj Denoted as C' ij and B' kjFor example, the actual indicator values and the value intervals corresponding to each level of risk can be normalized based on the following calculation method, positive indicators: Negative indicators: Where X represents the value to be normalized, max(X) and min(X) represent the maximum and minimum values of the range of values of the evaluation index corresponding to the value to be normalized. For example, when normalizing the actual index value corresponding to rock durability, X represents the actual index value, and max(X) and min(X) are 10 and 0, respectively. It is also understandable that for positive indicators, the larger the normalized value, the better. Conversely, for negative indicators, the smaller the normalized value, the better.

[0064] For example, the above risk level standard classification reference table is normalized as follows:

[0065]

[0066] Step S30: Determine the subjective weight corresponding to each evaluation indicator based on the hierarchical analysis method, determine the objective weight corresponding to each evaluation indicator through the entropy weight method based on the normalized indicator values and the normalized value intervals, and determine the comprehensive weight corresponding to each evaluation indicator based on the subjective weight and objective weight corresponding to each evaluation indicator.

[0067] In this step, the comprehensive weighting method is used to determine the comprehensive weight corresponding to each evaluation indicator. The comprehensive weighting method includes the hierarchical analysis method and the entropy weight method. That is, the subjective weight corresponding to each evaluation indicator is determined based on the hierarchical analysis method, and the objective weight corresponding to each evaluation indicator is determined based on the entropy weight method.

[0068] Exemplarily, determining the subjective weight corresponding to each of the evaluation indicators based on the hierarchical analysis method may include the following steps: determining the relative importance of each evaluation indicator, and constructing a judgment matrix based on the relative importance; if the judgment matrix meets the consistency test requirements, determining the maximum eigenvalue in the judgment matrix, and calculating the eigenvector corresponding to the maximum eigenvalue; and calculating the subjective weight corresponding to each of the evaluation indicators based on the eigenvector.

[0069] Based on the normalized indicator values and the normalized value intervals, the objective weight corresponding to each evaluation indicator is determined by the entropy weight method, including: determining the upper limit value of the normalized value interval corresponding to the risk level to which each normalized indicator value belongs; calculating the contribution degree corresponding to each evaluation indicator based on each normalized indicator value and the upper limit value of each normalized value interval; calculating the entropy corresponding to each evaluation indicator based on the contribution degree corresponding to each evaluation indicator; and calculating the objective weight corresponding to each evaluation indicator based on the entropy corresponding to each evaluation indicator.

[0070] In a specific embodiment of determining subjective weights using the hierarchical analysis method, the following steps are included: S1: establishing a hierarchical structure. S2: selecting a 1-9 scale as a reference standard for assigning values, and comparing evaluation indicators at the same level in pairs; six experts in the fields of engineering geology and cultural relics protection can be invited to assign values to the relative importance of each indicator, and constructing a column judgment matrix E = (e vp ) 8×8 , where e vp Indicates the importance of the vth indicator relative to the pth indicator. S3: Perform consistency test on the judgment matrix and calculate the consistency index CI and consistency ratio CR; if CR < 0.1, the judgment matrix meets the consistency test requirements; the specific formulas of consistency index CI and consistency ratio CR are as follows: where λ max is the maximum eigenvalue of the judgment matrix, n is the order of the judgment matrix, and RI is the random consistency index. Step S4: Calculate the eigenvector corresponding to the maximum eigenvalue of each judgment matrix and normalize the eigenvector to obtain the weight coefficient of each factor at each level; multiply the weight coefficient of the same evaluation index at each level to obtain the subjective weight w j1 .

[0071] In addition, when conducting risk assessment on multiple grotto temples to be evaluated and using the hierarchical analysis method to determine the subjective weight, the normalized data C' can also be used. ij and B' kj Construct matrix A:

[0072] In a specific embodiment of determining the objective weight using the entropy weight method, the following steps are included: S1: According to the definition of entropy, the entropy H of each evaluation index is calculated. j : Among them, f ij represents the contribution of the jth evaluation indicator in the i-th grotto temple to be evaluated, n represents the total number of grotto temples to be evaluated, V represents the number of risk levels, and V takes values of Ⅰ, Ⅱ, Ⅲ, Ⅳ or Ⅴ, f kj Indicates the contribution of the jth evaluation indicator in the kth level risk, C' ij represents the normalized index value corresponding to the jth evaluation index in the i-th grotto temple to be evaluated, B' kj Indicates the upper limit of the normalized value range of the jth evaluation indicator in the kth level risk. Step S2: Calculate the entropy weight (objective weight) w of the jth evaluation indicator j2 :

[0073] Furthermore, the calculation formula of the comprehensive weight is: w j =ξwj1 +(1-ξ)w j2 ; Among them, w j1 represents the subjective weight of evaluation index j, w j2 represents the objective weight of the evaluation index j, ξ = 0.5.

[0074] Step S40: Based on the comprehensive weights, the normalized index values and the normalized value intervals corresponding to the risks at each level, a risk level assessment result corresponding to the grotto temple to be evaluated is determined by a distance discrimination method.

[0075] In this step, the risk level assessment result corresponding to the grotto temple to be evaluated is determined based on the distance discrimination method. Exemplarily, this step may include: calculating the two-dimensional Euclidean distance between the grotto temple to be evaluated and each risk level based on the comprehensive weights, the normalized index values, and the upper limits of the normalized value ranges; and determining the risk level assessment result corresponding to the grotto temple to be evaluated based on the two-dimensional Euclidean distance values using a clamping criterion.

[0076] In a specific embodiment, when determining the risk level assessment result corresponding to the grotto temple to be evaluated based on the distance discrimination method, the two-dimensional Euclidean distance D between the upper limit of the value range of the five risk levels and the i-th grotto temple to be evaluated is calculated respectively. k : C' ij and B' kj They are respectively the C mentioned in the above steps ij and B kj The value obtained after normalization. Further, based on D k The minimum value and the second-to-last minimum value in the risk level of the cave temple to be evaluated are determined by the squeeze criterion; for example, when the distance D Ⅱ When it is the minimum value, the risk level may be Ⅰ or Ⅱ. Further, when the distance D between the first-level risk Ⅰ When the distance D between the risk level and the fifth level is calculated, the final risk level is determined to be level I. Ⅴ When it is the minimum value, further calculation is required When D' is the second smallest value, it is determined that the cave temple belongs to the V level risk, and when D Ⅳ When it is the second smallest value, it belongs to Level IV risk.

[0077] In the non-destructive assessment method of the instability risk level of the grotto temple in the above embodiment, not only the intrinsic characteristics of the grotto temple are covered, but also the degradation state of the rock and rock structure, as well as the occurrence environment of the grotto temple are deeply considered; this multi-dimensional evaluation system can more accurately reflect the actual situation of the stability of the grotto temple and provide solid data support for the classification of risk levels.

[0078] Figure 2 FIG. 1 is a flow chart of a non-destructive assessment method for the instability risk level of a cave temple according to another embodiment of the present invention. Figure 2 As shown in the figure, the non-destructive assessment method specifically includes the following steps: S01: constructing an index system for evaluating the risk level of grotto temples based on rock durability, occurrence environment, cave structure stability, rock physical properties, rock integrity and disease development degree; step S02: dividing the above indicators into five levels, and determining a grading reference table; step S03: based on on-site and indoor investigations, ascertaining the rock type, occurrence environment and cave shape of the grotto temple to be evaluated, and determining the level and specific value of each indicator; step S04: obtaining rock physical properties, rock integrity and disease development degree based on on-site tests, and determining the level of each indicator; step S05: normalizing the index data and the grading reference table, and determining the weight of each indicator using a comprehensive weighting method; step S06: determining the risk level of the grotto temple to be evaluated based on the indicator weight and the distance discrimination method.

[0079] Through the above embodiments, it can be found that the non-destructive assessment method of the instability risk level of the grotto temple of the present invention adopts advanced instruments such as wave velocity tester, portable Leeb hardness tester, infrared thermal imager, etc. when obtaining the actual index value of the grotto temple. The above instruments have the characteristics of non-destructive, flexible, convenient, fast and real-time. In the actual testing process, there is no need to perform destructive operations on the grotto temple, avoiding secondary damage to the grotto temple; this non-destructive testing method not only protects the integrity of the grotto temple, but also greatly improves the testing efficiency and safety. In addition, the application adopts a comprehensive weighting method to determine the weight of each evaluation index, making the weight distribution more scientific and reasonable, avoiding the one-sidedness and arbitrariness of subjective weighting; at the same time, combined with the distance discrimination method and the squeeze criterion, the risk level evaluation result is more accurate and reliable. The non-destructive assessment method of the instability risk level of the grotto temple of the present application not only provides strong support for the protection of the grotto temple, but also provides new ideas and methods for research and practice in related fields.

[0080] According to another aspect of the present invention, a non-destructive assessment system for the instability risk level of cave temples is also provided. The system includes a processor, a memory, and a computer program stored in the memory. The processor is used to execute the computer program. When the computer program is executed, the system implements the steps of the method described in any of the above embodiments.

[0081] Embodiments of the present invention further provide a computer-readable storage medium and a computer program product, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the method described in any of the above embodiments are implemented. The computer-readable storage medium can be a tangible storage medium, such as a random access memory (RAM), a memory, a read-only memory (ROM), an electrically programmable ROM, an electrically erasable programmable ROM, a register, a floppy disk, a hard disk, a removable storage disk, a CD-ROM, or any other form of storage medium known in the art.

[0082] It should be understood by those skilled in the art that the various exemplary components, systems and methods described in conjunction with the embodiments disclosed herein can be implemented in hardware, software or a combination of the two. Whether it is specifically performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention. When implemented in hardware, it can be, for example, an electronic circuit, an application specific integrated circuit (ASIC), appropriate firmware, a plug-in, a function card, etc. When implemented in software, the elements of the present invention are programs or code segments that are used to perform the required tasks. The program or code segment can be stored in a machine-readable medium, or transmitted on a transmission medium or a communication link via a data signal carried in a carrier.

[0083] It should be understood that the present invention is not limited to the specific configurations and processes described above and illustrated in the figures. For the sake of brevity, a detailed description of known methods is omitted. In the above embodiments, several specific steps are described and illustrated as examples. However, the method of the present invention is not limited to the specific steps described and illustrated. Those skilled in the art may make various changes, modifications, and additions, or change the order of the steps after understanding the spirit of the present invention.

[0084] In the present invention, features described and / or illustrated for one embodiment may be used in the same or similar manner in one or more other embodiments, and / or combined with or replace features of other embodiments.

[0085] The foregoing description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Those skilled in the art will readily appreciate that various modifications and variations to the present invention are possible. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present invention are intended to be within the scope of protection of the present invention.

Claims

1. A non-destructive assessment method for the instability risk level of a cave temple, characterized by: The method comprises: Determine the risk level standards corresponding to each evaluation indicator, including rock durability, occurrence environment, cave structure stability, rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock mass integrity and degree of disease development. The risk level standards include the value range corresponding to each level of risk; Determine the actual indicator value corresponding to each evaluation indicator of the grotto temple to be evaluated, and normalize the actual indicator value and each value interval to obtain each normalized indicator value and each normalized value interval; Determine the subjective weight corresponding to each evaluation indicator based on the hierarchical analysis method, determine the objective weight corresponding to each evaluation indicator based on the normalized indicator value and each normalized value interval through the entropy weight method, and determine the comprehensive weight corresponding to each evaluation indicator based on the subjective weight and objective weight corresponding to each evaluation indicator; Determine the risk level assessment result corresponding to the grotto temple to be evaluated by a distance discrimination method based on the comprehensive weights, the normalized index values and the normalized value intervals corresponding to the risk levels; Determining the objective weight corresponding to each evaluation indicator based on the normalized indicator values and the normalized value intervals by the entropy weight method includes: Determine the upper limit of the normalized value interval corresponding to the risk level to which each normalized indicator value belongs; Calculate the contribution degree of each evaluation indicator based on the value of each normalized indicator and the upper limit of each normalized value interval; Calculating the entropy corresponding to each evaluation indicator based on the contribution degree corresponding to each evaluation indicator; Calculating the objective weight corresponding to each evaluation indicator based on the entropy corresponding to each evaluation indicator; The calculation formula of the entropy is: Among them, f ij represents the contribution of the jth evaluation indicator in the i-th grotto temple to be evaluated, n represents the total number of grotto temples to be evaluated, V represents the number of risk levels, and V takes values of I, II, III, IV or V, f kj Indicates the contribution of the jth evaluation indicator in the kth level risk, C′ ij represents the normalized index value corresponding to the jth evaluation index in the i-th grotto temple to be evaluated, B' kj It represents the upper limit of the normalized value range of the jth evaluation indicator in the kth level risk; The objective weight w of the jth evaluation index j2 : Based on the comprehensive weights, the normalized index values, and the normalized value intervals corresponding to the risks at each level, a risk level assessment result corresponding to the grotto temple to be evaluated is determined by a distance discrimination method, including: Based on the comprehensive weights, the normalized index values and the upper limit of the normalized value interval, the two-dimensional Euclidean distance between the grotto temple to be evaluated and each risk level is calculated; the two-dimensional Euclidean distance D k : w j is the comprehensive weight; Based on the two-dimensional Euclidean distance value, the risk level assessment result corresponding to the grotto temple to be evaluated is determined by the squeeze criterion; based on D k The minimum value and the second-to-last minimum value in the risk level assessment result of the cave temple to be evaluated are determined by using the squeeze criterion. When the distance D between the two risks is II is the minimum value, and the distance D between it and the first-level risk I When it is the second minimum value, the final risk level is determined to be Level I.

2. The non-destructive assessment method for the instability risk level of a cave temple according to claim 1 is characterized in that: Determine the actual indicator values corresponding to each evaluation indicator of the grotto temple to be evaluated, including: Determine the actual rock type, actual occurrence environment, and actual cave shape of the grotto temple to be evaluated, and determine the actual index values of the rock durability, occurrence environment, and cave structure stability evaluation indicators of the grotto temple to be evaluated based on the actual rock type, actual occurrence environment, and actual cave shape and the risk level standard; Conduct on-site testing on the grotto temple to be evaluated to obtain actual index values corresponding to the rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock integrity and disease development degree evaluation indicators.

3. The non-destructive assessment method for the instability risk level of a cave temple according to claim 1 is characterized in that: The subjective weight corresponding to each evaluation indicator is determined based on the hierarchical analysis method, including: Determine the relative importance of each evaluation indicator, and construct a judgment matrix based on the relative importance of each evaluation indicator; If the judgment matrix meets the consistency test requirements, determining the maximum eigenvalue in the judgment matrix, and calculating the eigenvector corresponding to the maximum eigenvalue; The subjective weight corresponding to each evaluation indicator is calculated based on the feature vector.

4. The non-destructive assessment method for the instability risk level of a cave temple according to claim 3 is characterized in that: The calculation formula of the comprehensive weight is: w j =ξw j1 +(1-ξ)w j2 ; Among them, w j1 represents the subjective weight of evaluation index j, w j2 represents the objective weight of the evaluation index j, ξ = 0.

5.

5. The non-destructive assessment method for the instability risk level of a cave temple according to claim 1 is characterized in that: Determine the risk level standards corresponding to each evaluation indicator, including: Based on the assignment method, the rock durability, occurrence environment, and cave structure stability evaluation indicators are divided into five levels respectively; Based on the empirical method, the rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock integrity and disease development degree evaluation indicators are divided into five levels respectively.

6. A non-destructive assessment system for the instability risk level of a grotto temple, characterized by: The system includes a processor, a memory, and a computer program stored in the memory, wherein the processor is configured to execute the computer program, and when the computer program is executed, the system implements the steps of the method according to any one of claims 1 to 5.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 5 are implemented.

Images