Lossless assessment method and device for instability risk level of stone cave temple

Abstract

The invention provides a lossless assessment method and device for a grotto temple instability risk grade, and the method comprises the steps: determining a risk grade standard corresponding to each evaluation index, and the risk grade standard comprises a value interval corresponding to each grade risk; determining an actual index value corresponding to each evaluation index of the stone cave temple to be evaluated, and performing normalization processing on each actual index value and each value interval to obtain each normalized index value and each normalized value interval; determining a subjective weight corresponding to each evaluation index based on an analytic hierarchy process, determining an objective weight corresponding to each evaluation index through an entropy weight method, and determining a comprehensive weight corresponding to each evaluation index based on the subjective weight and the objective weight corresponding to each evaluation index; and determining a risk grade evaluation result corresponding to the to-be-evaluated grotto temple through a distance discrimination method based on each comprehensive weight, each normalized index numerical value and the normalized value interval corresponding to each grade risk. According to the method, the accuracy of the evaluation result of the grotto temple instability risk is improved.

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, as a brilliant heritage of Chinese culture, have not only witnessed the changes of history, but also carried profound cultural and artistic values. However, under the long-term influence of natural forces and human activities, the rock mass of the grottoes is prone to deterioration such as cracking, pulverization, and peeling, which seriously threatens the long-term preservation of the grotto temples. In order to more effectively respond to these challenges, it is urgent to carry out in-depth and systematic protection research, among which the scientific assessment of the instability risk of the grotto rock mass is particularly critical. Considering the diversity of grotto rock properties, occurrence environments, and cave forms, as well as the diversity of deterioration and instability characteristics, exploring a method that can assess the instability risk of grotto temples has great theoretical value and practical significance for guiding us to formulate targeted protection measures and extend 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, which makes the evaluation results of the stability of grotto temples inaccurate.

[0004] In summary, although the existing assessment methods for the instability risk of grotto temples can reflect the stability of grotto temples to a certain extent, they have limitations in accuracy. Therefore, how to improve the accuracy of the assessment results of the instability risk of grotto temples 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 standard corresponding to each evaluation index, wherein the evaluation index includes rock durability, occurrence environment, cave structure stability, rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock mass integrity and disease development degree, and the risk level standard includes 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 values ​​and the value intervals to obtain normalized indicator values ​​and normalized value intervals;

[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 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;

[0010] Based on the comprehensive weights, normalized indicator values ​​and normalized value intervals corresponding to each level of risk, 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 index value corresponding to each evaluation index 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 index 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 field tests 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 hierarchical analysis method 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, determine the maximum eigenvalue in the judgment matrix, and calculate the eigenvector corresponding to the maximum eigenvalue;

[0017] The subjective weight corresponding to each of the evaluation indicators 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 normalized indicator value and the upper limit of each normalized value interval;

[0021] Calculate 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 index 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 is Ⅰ, Ⅱ, Ⅲ, Ⅳ 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 risks of each level, 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 between the grotto temple to be evaluated and each risk level is 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 cave 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 cave temples, but also deeply considers the degradation state of rocks and rock structures, as well as the occurrence environment in which the cave temples are located. This multi-dimensional evaluation system can more accurately reflect the actual situation of the stability of the cave temples, provide solid data support for the division of risk levels, and thus improve the accuracy of the assessment results of the instability risk of the cave 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, purposes, and features of the present invention will be described in part in the following description, and will become apparent to those skilled in the art after studying the following, or may be learned from the practice of the present invention. The purposes and other advantages of the present invention may be achieved and obtained by the structures specifically indicated in the specification and the accompanying drawings.

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

[0041] The drawings described herein are used 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 only for illustrating the principles of the present invention. In order to facilitate the illustration and description of some parts of the present invention, the corresponding parts in the drawings may be enlarged, that is, they may become larger relative to other components in the exemplary device actually manufactured according to the present invention. In the drawings:

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

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

[0044] In order to make the purpose, technical solution 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 illustrative 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, only structures and / or processing steps closely related to the solutions according to the present invention are shown in the accompanying drawings, 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 presence of features, elements, steps or components, but does not exclude the presence 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, 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 standard corresponding to each evaluation index, wherein the evaluation index includes 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, and the risk level standard includes the value range corresponding to each level of risk.

[0051] In this step, the risk level standard corresponding to each evaluation index is constructed as the evaluation reference system of the risk level of the cave temple. The rock longitudinal wave velocity, rock heat absorption coefficient and rock Leeb hardness evaluation indexes can also be collectively referred to as rock physical property evaluation indexes. Specifically, the number of risk levels divided into multiple evaluation indicators can be the same, that is, 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 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 index is divided into five levels: I, II, III, IV, and V, wherein the risk of instability of the cave temple is the lowest under the I level state and the risk of instability of the cave temple is the highest under the V level state. Exemplarily, the total number of risk levels is five, and the risk level standard corresponding to each evaluation index is determined, which may specifically 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 level 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. 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 column cave, bucket-top cave, gable-top cave, arch cave and flat-top cave are classified as level 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 storage environment of the cave temple, and the values ​​are recorded as F, F, F, F and F. q 、F d 、F s and F h , then the comprehensive score of the habitat 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: Ⅰ, Ⅱ, Ⅲ, Ⅳ, and Ⅴ. The corresponding value ranges of 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 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 Ⅰ, Ⅱ, Ⅲ, Ⅳ, and Ⅴ, and each level has a corresponding value range. The following table is a grading reference table of risk level standards for an embodiment of the present application:

[0055]

[0056] Step S20: Determine the actual indicator value corresponding to each evaluation indicator of the grotto temple to be evaluated, and normalize the actual indicator values ​​and the value intervals to obtain normalized indicator values ​​and normalized value intervals.

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

[0058] In one embodiment, determining the actual index value corresponding to each evaluation index of the grotto temple to be evaluated may specifically include: determining the actual rock type, actual occurrence environment and actual cave form of the grotto temple to be evaluated, and determining the actual index value of the rock durability, occurrence environment and cave structure stability evaluation index of the grotto temple to be evaluated based on the actual rock type, actual occurrence environment and actual cave form and the risk level standard; conducting on-site testing on the grotto temple to be evaluated to obtain the actual index values ​​corresponding to the rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock integrity and disease development degree evaluation index. In this embodiment, the rock type, occurrence environment and cave form of the grotto temple to be evaluated can be ascertained through on-site investigation and other methods, and then the risk level corresponding to the rock type, occurrence environment and cave form of the grotto temple to be evaluated is pre-defined based on the investigation results, and then the corresponding value range is determined from the graded reference table based on the risk level corresponding to each index, and finally the expert selects a suitable value from the corresponding value range based on experience as the actual index value of the corresponding evaluation index. 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; and the rock integrity coefficient Kv is determined by the longitudinal wave velocity V of the rock mass. m 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 They are the wave velocity of the rock mass and the wave velocity of the rock block of the grotto temple to be evaluated, measured by the RSM-SY6(C) non-metallic ultrasonic detector. In addition, the development area of ​​water erosion, weathering and biological diseases in the degree of disease development can be obtained through on-site statistical measurement.

[0060] Specifically, when using the RSM-SY6 (C) non-metallic ultrasonic detector to measure the wave velocity of the cave 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 connected to measure the delay of the instrument, and this delay time is used as the zero calibration 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 sensor probe to measure directly on the rock surface. In addition, after the rock block wave velocity and rock mass wave velocity of the grotto temple to be evaluated are measured, 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 T of the complete area. 0 ; The heat flux density q and initial temperature T of the entire area 0 Substitute time t 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] When obtaining the actual index value of the rock Leeb hardness index of the grotto temple to be evaluated, the TIME5350 Leeb hardness tester is used for measurement, and the D-type universal probe is preferred; the instrument is small in size, light in weight, and simple to test, and has little impact on the surface of the grotto rock mass during the test. In addition, when obtaining the actual index value of the disease development degree index of the grotto temple to be evaluated, the development area of ​​water erosion disease, weathering disease and biological disease of the grotto temple can be determined through manual on-site measurement.

[0063] After obtaining the actual index values ​​of the grotto temples to be evaluated, the actual index 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 jth index 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 = I, ..., V) is recorded as B kj That is, B kj represents 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' kj For example, the actual indicator value and the value interval corresponding to each level of risk can be normalized based on the following calculation method, positive indicator: 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 the actual index value corresponding to rock durability is normalized, X represents the actual index value, and max(X) and min(X) are 10 and 0, respectively. It can also be understood that for positive indicators, the larger the normalized value, the better, on the contrary, for negative indicators, the smaller the normalized value, the better.

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

[0065]

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

[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 of the evaluation indicators, 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 weights corresponding to the evaluation indicators are determined by the entropy weight method, including: determining the upper limit value of the normalized value interval corresponding to the risk level to which the normalized indicator value belongs; calculating the contribution degree corresponding to each evaluation indicator based on the normalized indicator values ​​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 the 1-9 scale method as a reference standard for assigning values, and comparing the 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 requirements of consistency test; 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 weights, the normalized data C' ij and B' kj Construct the 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 index 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 is Ⅰ, Ⅱ, Ⅲ, Ⅳ 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 represents the upper limit of the normalized value interval 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: 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.

[0074] Step S40: 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 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 value between the grotto temple to be evaluated and each risk level based on each comprehensive weight, each normalized index value and each normalized value interval upper limit value; and determining the risk level assessment result corresponding to the grotto temple to be evaluated by the clamping criterion based on the two-dimensional Euclidean distance value.

[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 value of the value interval of the five risk levels and the i-th grotto temple to be evaluated is calculated respectively. k : C' ij and B' kj 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 clamping criterion; for example, when the distance D Ⅱ When it is the minimum value, the risk level may be Ⅰ or Ⅱ. When the distance D between the first-level risk Ⅰ is the second minimum value, the final risk level is determined to be level I. It should also be noted that when the distance D between the calculated and level 5 risk is Ⅴ 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. Ⅳ 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 temples in the above-mentioned embodiment, not only the main characteristics of the grotto temples are covered, but also the degradation state of the rock and rock structure, as well as the existence environment of the grotto temples are deeply considered; this multi-dimensional evaluation system can more accurately reflect the actual situation of the stability of the grotto temples, 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, 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: finding out the rock type, occurrence environment and cave shape of the grotto temple to be evaluated based on on-site and indoor investigations, 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 index weight and 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 test 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 test efficiency and safety. In addition, the application adopts a comprehensive weighting method to determine the weights 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 clamping criterion, the evaluation result of the risk level 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 a cave temple 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] The embodiments of the present invention also 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 technical field.

[0082] It should be understood by those skilled in the art that the 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 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 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 clear that the present invention is not limited to the specific configuration and processing described above and shown in the figures. For the sake of simplicity, a detailed description of the known method is omitted here. In the above embodiments, several specific steps are described and shown as examples. However, the method process of the present invention is not limited to the specific steps described and shown, and those skilled in the art can make various changes, modifications and additions, or change the order between 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 features of other embodiments or replace features of other embodiments.

[0085] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the embodiments of the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A non-destructive assessment method for the risk level of instability of a cave temple, characterized in that: The method comprises: Determine the risk level standard corresponding to each evaluation index, wherein the evaluation index includes rock durability, occurrence environment, cave structure stability, rock longitudinal wave velocity, rock heat absorption coefficient, rock Leeb hardness, rock mass integrity and disease development degree, and the risk level standard includes 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 values ​​and the value intervals to obtain normalized indicator values ​​and normalized value intervals; 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; Based on the comprehensive weights, normalized indicator values ​​and normalized value intervals corresponding to each level of risk, the risk level assessment result corresponding to the grotto temple to be evaluated is determined by the distance discrimination method.

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 index values ​​corresponding to each evaluation index 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 index 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 field tests 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 weights corresponding to the evaluation indicators are 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, determine the maximum eigenvalue in the judgment matrix, and calculate the eigenvector corresponding to the maximum eigenvalue; The subjective weight corresponding to each of the evaluation indicators 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: Based on the normalized index values ​​and the normalized value intervals, the objective weights corresponding to the evaluation indicators are determined by the entropy weight method, including: 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 normalized indicator value and the upper limit of each normalized value interval; Calculate the entropy corresponding to each evaluation indicator based on the contribution degree corresponding to each evaluation indicator; The objective weight corresponding to each evaluation indicator is calculated based on the entropy corresponding to each evaluation indicator.

5. The non-destructive assessment method for the instability risk level of a cave temple according to claim 4 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 the evaluation index, w j2 Represents the objective weight of the evaluation index, ξ=0.

5.

6. The non-destructive assessment method for the instability risk level of a cave temple according to claim 4 is characterized in that: The calculation formula of the entropy is: Among them, f ij represents the contribution of the ith evaluation index in the ith 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 can be I, II, III, IV or V, f kj Indicates the contribution of the ,th evaluation index in the kth level risk, C ij represents the normalized index value corresponding to the ith evaluation index in the ith grotto temple to be evaluated, B' kj It represents the upper limit of the normalized value range of the evaluation indicator in the k-th level risk.

7. 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.

8. The non-destructive assessment method for the instability risk level of a cave temple according to claim 1 is characterized in that: 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 a distance discrimination method, including: Based on the comprehensive weights, the normalized index values ​​and the upper limits of the normalized value intervals, the two-dimensional Euclidean distance between the grotto temple to be evaluated and each risk level is calculated; 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.

9. A non-destructive assessment system for the risk level of instability of cave temples, characterized in that: The system comprises a processor, a memory and a computer program stored in the memory, wherein the processor is used 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 8.

10. 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 8 are implemented.

Images