CN113704849A - Evaluation method for power transmission and transformation project slope instability risk degree - Google Patents

Evaluation method for power transmission and transformation project slope instability risk degree Download PDF

Info

Publication number
CN113704849A
CN113704849A CN202110986635.3A CN202110986635A CN113704849A CN 113704849 A CN113704849 A CN 113704849A CN 202110986635 A CN202110986635 A CN 202110986635A CN 113704849 A CN113704849 A CN 113704849A
Authority
CN
China
Prior art keywords
index
sub
data
indexes
power transmission
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110986635.3A
Other languages
Chinese (zh)
Other versions
CN113704849B (en
Inventor
李熙
陈垚
江世雄
李昕妍
王重卿
车艳红
翁孙贤
程慧青
陈震平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd
Original Assignee
Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd filed Critical Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd
Priority to CN202110986635.3A priority Critical patent/CN113704849B/en
Publication of CN113704849A publication Critical patent/CN113704849A/en
Application granted granted Critical
Publication of CN113704849B publication Critical patent/CN113704849B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

The invention relates to an evaluation method of the instability risk degree of a side slope of a power transmission and transformation project, which comprises the following steps of setting indexes for evaluating the instability risk degree of the side slope of the power transmission and transformation project, and collecting data of each index; when data of each index is collected, if the index is a descriptive index, scoring the index by an expert to obtain the data of the index; carrying out standardization processing on the data of each index to eliminate the difference among the data of each index; screening each index, and rejecting indexes with small information amount and high repeated information; determining the weight coefficient corresponding to each index and weighting to obtain a comprehensive index score; forming a corresponding relation between the comprehensive index score and the instability risk level; and obtaining the comprehensive index score of the target power transmission and transformation project slope and the instability risk level thereof. The invention evaluates the instability risk of the side slope by a plurality of indexes related to the stability of the side slope and an index processing method. The method can be used for slope instability early warning or site selection in the power transmission and transformation project according to comprehensive evaluation of slope indexes.

Description

Evaluation method for power transmission and transformation project slope instability risk degree
Technical Field
The invention relates to an evaluation method for slope instability risk degree of power transmission and transformation engineering, and belongs to the technical field of slope instability risk early warning.
Background
The central element of the geological environment of the power transmission and transformation project is the power transmission line project or a transformer substation, and for the development of the power grid technology, various geological conditions which influence the activities of the power transmission and transformation project belong to the geological environment of the power transmission and transformation project. The power transmission and transformation project is most easily influenced by geological environment factors due to the characteristics of long wiring and wide span. According to the relevant data, the natural environmental factors causing the failure of the transmission line are mainly as follows: high temperature, low temperature, lightning, ice and snow, storm wind, earthquake, landslide, debris flow, and the like. The geological disasters directly impacting the power transmission line tower and the foundation mainly comprise ice disasters, strong storms, earthquakes, floods, landslide and other geological disasters. In the prior art, an expert is adopted to score indexes, and then comprehensive evaluation is carried out according to scoring results. The method has the disadvantages that the scoring is greatly influenced by subjective factors, and the instability risk degree of each side slope cannot be transversely compared. And is not beneficial to forming a slope risk management method and site selection basis of a system.
Application document with the application number of CN202011392937.X discloses a method for evaluating technical states of slope engineering. And evaluating the slope state by collecting indexes related to the slope state. The method has the disadvantages that the stability of the side slope is not evaluated, and the side slope instability early warning effect is achieved through evaluation.
Disclosure of Invention
In order to overcome the problems, the invention provides an evaluation method for the instability risk degree of a side slope in power transmission and transformation engineering. The method can be used for slope instability early warning or site selection in the power transmission and transformation project according to comprehensive evaluation of slope indexes.
The technical scheme of the invention is as follows:
a method for evaluating the instability risk degree of a side slope of a power transmission and transformation project comprises the following steps,
setting a plurality of indexes for evaluating the instability risk degree of the power transmission and transformation project slope, and acquiring data of each index; the indexes comprise a landform index, a lithology index, a structure index, a meteorological condition index, a bad geological action index and an artificial activity index; when the data of each index is collected, if the index is a descriptive index, the index is scored by an expert to be used as the data of the index;
respectively carrying out standardization processing on the data of each index to eliminate the difference among the data of each index;
screening each index, and eliminating indexes with small information quantity and high repeated information;
determining a weight coefficient corresponding to each index and weighting to obtain a comprehensive index score;
repeating the steps to obtain comprehensive index scores of a plurality of power transmission and transformation project side slopes;
grading the instability risk degree of the power transmission and transformation project side slope according to each comprehensive index score to form a corresponding relation between the comprehensive index score and the instability risk grade;
and collecting data of each index of the target power transmission and transformation project slope to obtain a comprehensive index score and a instability risk level thereof.
Preferably, the topographic and geomorphic indexes comprise three sub-indexes of slope height, slope and vegetation coverage; the lithology and structure index comprises two sub-indexes of a rock-soil body characteristic and a geological structure; the meteorological condition indexes comprise three sub-indexes of precipitation, icing thickness and wind speed; the adverse geological action index comprises two sub-indexes of earthquake and landslide; the manual activity index comprises a manual excavation degree sub-index.
Further, the normalizing the data of each index includes the following steps:
step S1, selecting a sub-index which is not subjected to standardization processing to carry out standardization processing;
step S2, determining whether the sub-indicator is a numerical indicator, if yes, executing step S3, if no, considering that the sub-indicator has been standardized, and executing step S5;
step S3, judging whether the data of the sub-index has an outlier, if so, correcting the outlier, and executing step S4, otherwise, executing step S4;
step S4, standardizing the data of the sub-index to obtain the data score of the sub-index;
in step S5, it is determined whether or not all sub-indicators have been normalized, and if yes, the normalization process is terminated, and if not, step S1 is executed.
Further, the outlier is corrected in step S2 according to the following formula:
Figure BDA0003230920160000021
wherein x isiThe ith data of the sub-index; q1 and Q3 are the first quartile and the third quartile of the data of the sub-index, respectively; IQR is the interquartile range.
Further, the step of normalizing the data of the sub-index in step S3 includes the following steps:
step T1, establishing a plurality of index scoring functions; each index scoring function corresponds to different data distribution types and positive and negative respectively;
step T2, arranging the data of the sub-indexes in ascending order, establishing a coordinate system by taking the horizontal axis as a data serial number and the vertical axis as a data value, and fitting to obtain a data distribution function of the sub-indexes;
and T3, determining a matched index scoring function according to the data distribution function and the positive and negative of the index to obtain a basic scoring function, and obtaining a basic score corresponding to each data of the sub-index through the basic scoring function.
Further, the method further comprises the step of correcting the basic score, specifically comprising the following steps:
manually scoring each data of the sub-index to obtain a manual score and a manual score interval of each data of the sub-index;
and mapping the basic score to the manual score interval to obtain a corrected basic score.
Further, the screening of each of the indexes includes screening of indexes based on information amount, and specifically includes:
calculating the coefficient of variation CV of the sub-index, and the formula is as follows:
Figure BDA0003230920160000031
where s is the standard deviation of the sub-indices, xiIs the ith data of the sub-indicator, m is the data amount of the sub-indicator,
Figure BDA0003230920160000032
is the mean of the sub-indices;
setting a threshold T1The coefficient of variation CV is smaller than the threshold value T1Sub-indices of (1).
Further, the screening of each of the indexes includes correlation-based index screening, and specifically includes:
selecting any two sub-indexes, and arranging the data of the two sub-indexes respectively according to an ascending order;
calculating the correlation coefficient rho of the two sub-indexes, wherein the formula is as follows:
Figure BDA0003230920160000033
wherein, giAnd tiThe ith data of the two sub-indices,
Figure BDA0003230920160000034
and
Figure BDA0003230920160000035
respectively the mean values of the two sub-indexes;
repeating the steps until the correlation coefficient rho among all the sub-indexes is obtained;
setting a threshold T2The absolute value of the rejection correlation coefficient rho is more than T2One of the two sub-indices.
Further, the determining and weighting the weight coefficient corresponding to each index specifically includes:
g1, sorting the sub-indexes according to the importance of the sub-indexes on the instability of the power transmission and transformation project slope to obtain an importance sequence;
g2, manually endowing two adjacent sub-indexes in the importance sequence with relative importance riExpressed as follows:
Figure BDA0003230920160000036
wherein the content of the first and second substances,
Figure BDA0003230920160000037
and
Figure BDA0003230920160000038
respectively the importance degrees of the ith and the (i-1) th sub-indexes in the importance sequence;
step G3, calculating the importance of the sub-index with the lowest importance
Figure BDA0003230920160000039
n is the sub-index number, and the formula is as follows:
Figure BDA0003230920160000041
g4, according to the relative importance r of two adjacent sub indexesiAnd calculating to obtain the weight coefficient omega of the ith sub-indexi
Step G5, weighting to obtain a comprehensive index score Y, and according to the formula:
Figure BDA0003230920160000042
wherein x isjAnd ωjAnd respectively the basic score of the jth sub-index of a certain power transmission and transformation project slope and the weight coefficient corresponding to the sub-index.
Further, according to each of the comprehensive index scores, grade division is performed on the instability risk degree of the power transmission and transformation project side slope, and a corresponding relation between the comprehensive index score and the instability risk grade is formed, and the method comprises the following steps:
d1, determining a grading number K, and arranging the comprehensive index grades Y in an ascending order;
step D2, calculating the sum SDAM of the square deviation of the composite index score Y, wherein the formula is as follows:
Figure BDA0003230920160000043
wherein, YjFor the j-th composite index score,
Figure BDA0003230920160000044
the mean value of each comprehensive index score; m' is the number of the comprehensive index scores;
step D3, grading the comprehensive index score Y;
step D4, calculating the sum SDBC of the square deviation of the interstage comprehensive index score Y, wherein the formula is as follows:
Figure BDA0003230920160000045
wherein, YkiThe ith composite index score of k level,
Figure BDA0003230920160000046
is the mean value of the k-level comprehensive index scores, mkScoring the number of k-level comprehensive indexes;
step D5, calculating the goodness of fit of variance GVF, wherein the formula is as follows:
Figure BDA0003230920160000047
and D6, repeating the steps D3 to D5 until the variance goodness of fit under all grading conditions is obtained, and grading the slope instability risk degree of each power transmission and transformation project according to the grading scheme corresponding to the minimum variance goodness of fit when the variance goodness of fit is taken as the minimum value.
The invention has the following beneficial effects:
1. the evaluation method selects 5 indexes with the maximum relevance to the slope instability of the power transmission and transformation project and sub-indexes under each index to evaluate the slope instability risk degree, and early warning is carried out on the slope instability according to the evaluation result. The loss caused by the instability of the side slope of the power transmission and transformation project is reduced.
2. Compared with the common minimum-maximum standardization method, the evaluation method is changed to establish various index scoring functions. And matching the index scoring function according to the positive and negative of the index and the distribution function of the data. And obtaining a basic scoring function by adjusting parameters of the index scoring function. The basic scoring function is better in fit with the relation between the indexes and the risk degree.
3. Compared with an expert scoring method, the scoring range is easily overlarge through a basic scoring function, the scoring interval of the index is obtained through the expert scoring by the evaluation method, the basic scoring is mapped to the scoring interval, and the influence on comprehensive evaluation caused by the enlargement of the scoring range is reduced.
4. The evaluation method reduces the influence of outliers on the scoring structure by correcting the outliers.
5. The evaluation method eliminates indexes with small information quantity and high correlation, and improves the accuracy of evaluation.
Drawings
FIG. 1 is a schematic flow chart of the method of the present invention.
FIG. 2 is a schematic diagram of the index of the present invention.
FIG. 3 is a schematic diagram of a normalization process according to the present invention.
Detailed Description
The invention is described in detail below with reference to the figures and the specific embodiments.
Example one
Referring to fig. 1, a method for evaluating the degree of risk of side slope instability in power transmission and transformation engineering includes the following steps,
setting a plurality of indexes for evaluating the instability risk degree of the power transmission and transformation project slope, and acquiring data of each index; the indexes comprise a landform index, a lithology index, a structure index, a meteorological condition index, a bad geological action index and an artificial activity index; when the data of each index is collected, if the index is a descriptive index, the index is scored by an expert to be used as the data of the index;
because the actual numerical values of the indexes have differences, the specific condition of the power transmission and transformation project slope is not easy to know through the numerical values of the indexes, so that the data of the indexes need to be standardized respectively, and the differences among the data of the indexes are eliminated;
the indexes have the problems of small information quantity and high correlation with other indexes, if a certain index is always highly expressed in a research area, namely the information quantity is small, the index loses significance on the risk degree, the calculated amount is increased, and if the two indexes have certain correlation, repeated information is amplified, so that the accuracy of an evaluation result is influenced. Therefore, each index needs to be screened, and the indexes with small information quantity and high repeated information are removed;
determining a weight coefficient corresponding to each index and weighting to obtain a comprehensive index score;
repeating the steps to obtain comprehensive index scores of a plurality of power transmission and transformation project side slopes;
grading the instability risk degree of the power transmission and transformation project side slope according to each comprehensive index score to form a corresponding relation between the comprehensive index score and the instability risk grade;
and collecting data of each index of the target power transmission and transformation project slope to obtain a comprehensive index score and a instability risk level thereof.
Referring to fig. 2, in the present embodiment, the topographic indexes include three sub-indexes of slope height, slope and vegetation coverage; landslide is more likely to occur as the slope height is higher, so that great influence is generated on the construction of a transformer substation and the stability of a foundation, and the ground distance of the power transmission line is increased as the slope height is higher, so that the sensitivity of the power transmission line to meteorological conditions such as wind speed and ice coating is increased; the gradient is an important factor influencing the geological environment of the power transmission and transformation project. The slope of a side slope where a foundation of the power transmission and transformation project is located is large, and foundation damage is easy to occur; the vegetation coverage is divided into positive and negative aspects to the influence of power transmission and transformation engineering, and on the one hand, the vegetation coverage is high to be favorable to soil and water to keep, has reduced the probability that natural disasters such as landslide, mud-rock flow take place, and on the other hand, because the particularity of power transmission and transformation engineering, the vegetation coverage is high to be unfavorable for the construction of engineering and the maintenance of later stage. The lithology and structure index comprises two sub-indexes of a rock-soil body characteristic and a geological structure; the main influence factors of the characteristics of the rock-soil body on the power transmission and transformation engineering include the mechanical property of a foundation soil body, the stability of a rock stratum, the content and the influence degree of a special soil body; the address structure can influence the integrity of the rock-soil body structure, and the rock-soil body with poor integrity often causes slope instability, thereby influencing the safety of power transmission and transformation engineering. The meteorological condition indexes comprise three sub-indexes of precipitation, icing thickness and wind speed; the precipitation can determine the height of the water level of the underground water and the corrosivity of the underground water, and the influence on the stability of the foundation is large; the larger the icing thickness is, the larger the load loaded on the power transmission line is, extra tension can be generated on a pole tower, so that the instability of the pole tower is caused, and the ice coating of the line can cause the distance of the line to the ground to be reduced, so that the safe operation of the power transmission line is influenced; for the overhead transmission line with higher height, the wind load generated by strong wind acts on the transmission line, which is equal to the dynamic load applied to a plurality of electric lines, and the damage such as foundation pull-up is easy to cause. The adverse geological action index comprises two sub-indexes of earthquake and landslide; the earthquake can generate destructive damage to the power transmission and transformation project; the volume, speed and distance from the slope of the landslide all affect the stability of the foundation. The manual activity index comprises a manual excavation degree sub-index; the manual excavation degree directly influences the stability of the regional address environment of the power transmission and transformation engineering construction site, and the slope risk degree is higher in the region with the excessive manual excavation degree.
Referring to fig. 3, in the present embodiment, the normalizing the data of each index includes the following steps:
step S1, selecting a sub-index which is not subjected to standardization processing to carry out standardization processing;
step S2, determining whether the sub-indicator is a numerical indicator, if yes, executing step S3, if no, considering that the sub-indicator has been standardized, and executing step S5;
step S3, judging whether the data of the sub-index has an outlier, if so, correcting the outlier, and executing step S4, otherwise, executing step S4; compared with the expert scoring, the actual numerical value of the index may have an outlier, and the outlier can bring serious influence on the scoring structure, so that the score needs to be corrected;
step S4, standardizing the data of the sub-index to obtain the data score of the sub-index;
in step S5, it is determined whether or not all sub-indicators have been normalized, and if yes, the normalization process is terminated, and if not, step S1 is executed.
Example two
On the basis of the first embodiment, the method for evaluating the degree of slope instability risk of the power transmission and transformation project corrects the outlier in the step S2, and the formula is as follows:
Figure BDA0003230920160000071
wherein x isiIs the sub-indexThe ith data; q1 and Q3 are the first quartile and the third quartile of the data of the sub-index, respectively; IQR is the interquartile range. The correction method corrects the outlier to the maximum value or the minimum value of the normal interval, and reduces the influence of the outlier on the evaluation structure.
EXAMPLE III
On the basis of the first embodiment, the method for evaluating the instability risk degree of the side slope of the power transmission and transformation project standardizes the data of the sub-indexes in the step S3, and specifically comprises the following steps:
step T1, establishing a plurality of index scoring functions; each index scoring function corresponds to different data distribution types and positive and negative respectively; in reality, natural distribution of data mainly comprises three types of small value, large value and normal distribution, and six index scoring functions are needed by combining the natural distribution with the positive and negative of the data; in this embodiment, the index scoring function is as follows:
first metric scoring function:
Figure BDA0003230920160000072
second index scoring function:
Figure BDA0003230920160000073
third metric scoring function:
Figure BDA0003230920160000074
fourth metric scoring function:
Figure BDA0003230920160000081
fifth index scoring function:
Figure BDA0003230920160000082
sixth metric scoring function:
Figure BDA0003230920160000083
wherein, ynAnd xnRespectively the data of a certain sub-index of the nth power transmission and transformation project slope and the score, x, of the dataminAnd xmaxRespectively the minimum value and the maximum value of the data of the sub-index, alpha is a first curvature adjusting coefficient, gamma is a second curvature adjusting coefficient, beta is a horizontal direction adjusting coefficient, and epsilon and eta are a function vertical direction adjusting coefficient and a function vertical direction expansion coefficient respectively;
step T2, arranging the data of the sub-indexes in ascending order, establishing coordinates by taking the horizontal axis as a data serial number and the vertical axis as a data value, and fitting to obtain a data distribution function of the sub-indexes; obtaining the change condition of the data according to the distribution function, wherein the change condition comprises that the data is buffered first and then is fast, and the data is buffered first and then is fast and then is slow; the data of the small numerical value area is more after being buffered, the data of the large numerical value area is more after being buffered, and the data of the middle value is more after being buffered; these three distributions reflect the fact that data is concentrated in smaller values, in larger values and in normal distributions.
And T3, determining a matched index scoring function according to the data distribution function, adjusting parameters of the index scoring function according to the data distribution of the sub-indexes to obtain a basic scoring function, and obtaining basic scores corresponding to the data of the sub-indexes through the basic scoring function. When the index scoring function is determined, whether the sub-indexes are positively or negatively correlated with the risk degree, namely the positive and negative of the data, is judged, and then the matched index scoring function is determined according to the data distribution function of the sub-indexes. The distribution function of the positive indexes is a first index scoring function matched with the first index in a slow-first and fast-second mode, a second index scoring function matched with the first index scoring function in a fast-first and slow-second mode, a fifth index scoring function matched with the first index scoring function in a slow-first and fast-second mode, a third index scoring function matched with the negative indexes in a slow-first and fast-second mode, a fourth index scoring function matched with the first index scoring function in a fast-first and slow-second mode, and a sixth index scoring function matched with the first index scoring function in a slow-first and fast-second mode.
In this embodiment, the basic scores obtained through the basic score function are vertically distributed in 1 to 10, and compared with the expert scores, the distribution range of the basic scores is widened, and the actual scores should be concentrated in a certain area, so that the basic scores need to be corrected, in this embodiment, the basic scores are further corrected, specifically:
manually scoring each data of the sub-index to obtain a manual score and a manual score interval of each data of the sub-index;
and mapping the basic score to the manual score interval to obtain a corrected basic score.
Example four
On the basis of the third embodiment, the method for evaluating the instability risk degree of the power transmission and transformation project slope screens each index, including the index screening based on the information quantity, specifically comprises the following steps:
calculating the coefficient of variation CV of the sub-index, and the formula is as follows:
Figure BDA0003230920160000091
where s is the standard deviation of the sub-indices, xiIs the ith data of the sub-indicator, m is the data amount of the sub-indicator,
Figure BDA0003230920160000092
is the mean of the sub-indices;
setting a threshold T1The coefficient of variation CV is smaller than the threshold value T1Sub-indices of (1). The smaller the variation coefficient is, the more concentrated the value of the sub-index is, the smaller the amount of information contained. If the comprehensive evaluation includes too many indexes with small information amount, the accuracy of an evaluation system is influenced.
The screening of each index comprises correlation-based index screening, and specifically comprises the following steps:
selecting any two sub-indexes, and arranging the data of the two sub-indexes respectively according to an ascending order;
calculating the correlation coefficient rho of the two sub-indexes, wherein the formula is as follows:
Figure BDA0003230920160000093
wherein, giAnd tiThe ith data of the two sub-indices,
Figure BDA0003230920160000094
and
Figure BDA0003230920160000095
respectively the mean values of the two sub-indexes;
repeating the steps until the correlation coefficient rho among all the sub-indexes is obtained;
setting a threshold T2The absolute value of the rejection correlation coefficient rho is more than T2One of the two sub-indices. The correlation analysis method comprises Pearson correlation coefficient, Spearman correlation coefficient and Kendall correlation coefficient. The Pearson first relation number requires that the variable is normally distributed and is in a linear relation with the target, most of the influencing factors of the instability risk degree cannot meet the normal distribution, and the linear relation between the index and the risk degree part exists. The application object of the Kendall correlation coefficient is a categorical variable, and the Spearman correlation coefficient is applied to a numerical variable, so that the index is screened based on the Spearman correlation coefficient.
EXAMPLE five
On the basis of the fourth embodiment, the method for evaluating the degree of risk of slope instability of the power transmission and transformation project determines and weights the weight coefficients corresponding to the indexes, and specifically comprises the following steps:
g1, sorting the sub-indexes according to the importance of the sub-indexes on the instability of the power transmission and transformation project slope to obtain an importance sequence;
g2, artificially assigning two adjacent ones of the importance sequencesRelative importance of sub-indices riExpressed as follows:
Figure BDA0003230920160000101
wherein the content of the first and second substances,
Figure BDA0003230920160000102
and
Figure BDA0003230920160000103
respectively the importance degrees of the ith and the (i-1) th sub-indexes in the importance sequence;
step G3, calculating the importance of the sub-index with the lowest importance
Figure BDA0003230920160000104
n is the sub-index number, and the formula is as follows:
Figure BDA0003230920160000105
g4, according to the relative importance r of two adjacent sub indexesiAnd calculating to obtain the weight coefficient omega of the ith sub-indexi
Step G5, weighting to obtain a comprehensive index score Y, and according to the formula:
Figure BDA0003230920160000106
wherein x isjAnd ωjAnd respectively the basic score of the jth sub-index of a certain power transmission and transformation project slope and the weight coefficient corresponding to the sub-index.
The method comprises the following steps of grading the instability risk degree of the power transmission and transformation project side slope according to each comprehensive index score to form a corresponding relation between the comprehensive index score and the instability risk grade, wherein the grading comprises the following steps:
d1, determining a grading number K, and arranging the comprehensive index grades Y in an ascending order;
step D2, calculating the sum SDAM of the square deviation of the composite index score Y, wherein the formula is as follows:
Figure BDA0003230920160000107
wherein, YjFor the j-th composite index score,
Figure BDA0003230920160000108
the mean value of each comprehensive index score; m' is the number of the comprehensive index scores;
step D3, grading the comprehensive index score Y;
step D4, calculating the sum SDBC of the square deviation of the interstage comprehensive index score Y, wherein the formula is as follows:
Figure BDA0003230920160000109
wherein, YkiThe ith composite index score of k level,
Figure BDA00032309201600001010
is the mean value of the k-level comprehensive index scores, mkScoring the number of k-level comprehensive indexes;
step D5, calculating the goodness of fit of variance GVF, wherein the formula is as follows:
Figure BDA0003230920160000111
and D6, repeating the steps D3 to D5 until the variance goodness of fit under all grading conditions is obtained, and grading the slope instability risk degree of each power transmission and transformation project according to the grading scheme corresponding to the minimum variance goodness of fit when the variance goodness of fit is taken as the minimum value.
After the grade division is completed, the risk degree of the new power transmission and transformation project slope instability is judged, and only the comprehensive index grading band grading scheme is needed, and the grade does not need to be divided again. The grading scheme can also be updated through new slope data.
The above description is only an embodiment of the present invention, and not intended to limit the scope of the present invention, and all equivalent structures made by using the contents of the specification and the drawings of the present invention or directly or indirectly applied to other related technical fields are included in the scope of the present invention.

Claims (10)

1. A method for evaluating the instability risk degree of a side slope of a power transmission and transformation project is characterized by comprising the following steps,
setting a plurality of indexes for evaluating the instability risk degree of the power transmission and transformation project slope, and acquiring data of each index; the indexes comprise a landform index, a lithology index, a structure index, a meteorological condition index, a bad geological action index and an artificial activity index; when the data of each index is collected, if the index is a descriptive index, the index is scored by an expert to be used as the data of the index;
respectively carrying out standardization processing on the data of each index to eliminate the difference among the data of each index;
screening each index, and eliminating indexes with small information quantity and high repeated information;
determining a weight coefficient corresponding to each index and weighting to obtain a comprehensive index score;
repeating the steps to obtain comprehensive index scores of a plurality of power transmission and transformation project side slopes;
grading the instability risk degree of the power transmission and transformation project side slope according to each comprehensive index score to form a corresponding relation between the comprehensive index score and the instability risk grade;
and collecting data of each index of the target power transmission and transformation project slope to obtain a comprehensive index score and a instability risk level thereof.
2. The method for evaluating the instability risk degree of the power transmission and transformation project slope according to claim 1, wherein the topographic indicators comprise three sub-indicators of slope height, slope gradient and vegetation coverage; the lithology and structure index comprises two sub-indexes of a rock-soil body characteristic and a geological structure; the meteorological condition indexes comprise three sub-indexes of precipitation, icing thickness and wind speed; the adverse geological action index comprises two sub-indexes of earthquake and landslide; the manual activity index comprises a manual excavation degree sub-index.
3. The method for evaluating the slope instability risk degree of the power transmission and transformation project according to claim 1, wherein the step of standardizing the data of each index comprises the following steps:
step S1, selecting a sub-index which is not subjected to standardization processing to carry out standardization processing;
step S2, determining whether the sub-indicator is a numerical indicator, if yes, executing step S3, if no, considering that the sub-indicator has been standardized, and executing step S5;
step S3, judging whether the data of the sub-index has an outlier, if so, correcting the outlier, and executing step S4, otherwise, executing step S4;
step S4, standardizing the data of the sub-index to obtain the data score of the sub-index;
in step S5, it is determined whether or not all sub-indicators have been normalized, and if yes, the normalization process is terminated, and if not, step S1 is executed.
4. The method for evaluating the risk degree of slope instability of power transmission and transformation project according to claim 3, wherein the outlier is corrected in step S2 according to the following formula:
Figure FDA0003230920150000011
wherein x isiThe ith data of the sub-index; q1 and Q3 are the first quartile and the third quartile of the data of the sub-index, respectively; IQR is the interquartile range.
5. The method for evaluating the slope instability risk degree of the power transmission and transformation project according to claim 3, wherein the step S3 of normalizing the data of the sub-indexes comprises the following steps:
step T1, establishing a plurality of index scoring functions; each index scoring function corresponds to different data distribution types and positive and negative respectively;
step T2, arranging the data of the sub-indexes in ascending order, establishing a coordinate system by taking the horizontal axis as a data serial number and the vertical axis as a data value, and fitting to obtain a data distribution function of the sub-indexes;
and T3, determining a matched index scoring function according to the data distribution function and the positive and negative of the index to obtain a basic scoring function, and obtaining a basic score corresponding to each data of the sub-index through the basic scoring function.
6. The method for evaluating the risk degree of slope instability of power transmission and transformation project according to claim 5, further comprising correcting the basic score, specifically:
manually scoring each data of the sub-index to obtain a manual score and a manual score interval of each data of the sub-index;
and mapping the basic score to the manual score interval to obtain a corrected basic score.
7. The method for evaluating the slope instability risk degree of the power transmission and transformation project according to claim 6, wherein the screening of each index comprises screening of indexes based on information quantity, and specifically comprises the following steps:
calculating the coefficient of variation CV of the sub-index, and the formula is as follows:
Figure FDA0003230920150000021
where s is the standard deviation of the sub-indices, xiI-th data of sub-index, m is data of sub-indexThe amount of the compound (A) is,
Figure FDA0003230920150000022
is the mean of the sub-indices;
setting a threshold T1The coefficient of variation CV is smaller than the threshold value T1Sub-indices of (1).
8. The method for evaluating the slope instability risk degree of the power transmission and transformation project according to claim 7, wherein the screening of each index comprises correlation-based index screening, and specifically comprises the following steps:
selecting any two sub-indexes, and arranging the data of the two sub-indexes respectively according to an ascending order;
calculating the correlation coefficient rho of the two sub-indexes, wherein the formula is as follows:
Figure FDA0003230920150000031
wherein, giAnd tiThe ith data of the two sub-indices,
Figure FDA0003230920150000032
and
Figure FDA0003230920150000033
respectively the mean values of the two sub-indexes;
repeating the steps until the correlation coefficient rho among all the sub-indexes is obtained;
setting a threshold T2The absolute value of the rejection correlation coefficient rho is more than T2One of the two sub-indices.
9. The method for evaluating the slope instability risk degree of the power transmission and transformation project according to claim 8, wherein the determining and weighting of the weighting coefficient corresponding to each index specifically comprises:
g1, sorting the sub-indexes according to the importance of the sub-indexes on the instability of the power transmission and transformation project slope to obtain an importance sequence;
g2, manually endowing two adjacent sub-indexes in the importance sequence with relative importance riExpressed as follows:
Figure FDA0003230920150000034
wherein the content of the first and second substances,
Figure FDA0003230920150000035
and
Figure FDA0003230920150000036
respectively the importance degrees of the ith and the (i-1) th sub-indexes in the importance sequence;
step G3, calculating the importance of the sub-index with the lowest importance
Figure FDA0003230920150000037
n is the sub-index number, and the formula is as follows:
Figure FDA0003230920150000038
g4, according to the relative importance r of two adjacent sub indexesiAnd calculating to obtain the weight coefficient omega of the ith sub-indexi
Step G5, weighting to obtain a comprehensive index score Y, and according to the formula:
Figure FDA0003230920150000039
wherein x isjAnd ωjAnd respectively the basic score of the jth sub-index of a certain power transmission and transformation project slope and the weight coefficient corresponding to the sub-index.
10. The method for evaluating the instability risk level of the power transmission and transformation project slope according to claim 9, wherein the step of grading the instability risk level of the power transmission and transformation project slope according to each of the composite index scores to form a corresponding relationship between the composite index score and the instability risk level comprises the following steps:
d1, determining a grading number K, and arranging the comprehensive index grades Y in an ascending order;
step D2, calculating the sum SDAM of the square deviation of the composite index score Y, wherein the formula is as follows:
Figure FDA0003230920150000041
wherein, YjFor the j-th composite index score,
Figure FDA0003230920150000042
the mean value of each comprehensive index score; m' is the number of the comprehensive index scores;
step D3, grading the comprehensive index score Y;
step D4, calculating the sum SDBC of the square deviation of the interstage comprehensive index score Y, wherein the formula is as follows:
Figure FDA0003230920150000043
wherein, YkiThe ith composite index score of k level,
Figure FDA0003230920150000044
is the mean value of the k-level comprehensive index scores, mkScoring the number of k-level comprehensive indexes;
step D5, calculating the goodness of fit of variance GVF, wherein the formula is as follows:
Figure FDA0003230920150000045
and D6, repeating the steps D3 to D5 until the variance goodness of fit under all grading conditions is obtained, and grading the slope instability risk degree of each power transmission and transformation project according to the grading scheme corresponding to the minimum variance goodness of fit when the variance goodness of fit is taken as the minimum value.
CN202110986635.3A 2021-08-26 2021-08-26 Evaluation method for power transmission and transformation project slope instability risk degree Active CN113704849B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110986635.3A CN113704849B (en) 2021-08-26 2021-08-26 Evaluation method for power transmission and transformation project slope instability risk degree

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110986635.3A CN113704849B (en) 2021-08-26 2021-08-26 Evaluation method for power transmission and transformation project slope instability risk degree

Publications (2)

Publication Number Publication Date
CN113704849A true CN113704849A (en) 2021-11-26
CN113704849B CN113704849B (en) 2024-03-01

Family

ID=78655015

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110986635.3A Active CN113704849B (en) 2021-08-26 2021-08-26 Evaluation method for power transmission and transformation project slope instability risk degree

Country Status (1)

Country Link
CN (1) CN113704849B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115100819A (en) * 2022-05-23 2022-09-23 深圳市北斗云信息技术有限公司 Landslide hazard early warning method and device based on big data analysis and electronic equipment

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930348A (en) * 2012-10-19 2013-02-13 广东电网公司电力科学研究院 Assessment method for rainstorm disaster risks of sectional power transmission line pole-tower foundation slopes
CN103268420A (en) * 2013-05-24 2013-08-28 河海大学 Method for evaluating risks of high rock slope
KR101580062B1 (en) * 2014-12-05 2015-12-28 연세대학교 산학협력단 Method and system for real-time prediction and warning of landslides
US20190051146A1 (en) * 2017-08-09 2019-02-14 Institute Of Mountain Hazards And Environment, Chinese Academy Of Sciences Three-dimensional multi-point multi-index early warning method for risk at power grid tower in landslide section
CN109726902A (en) * 2018-12-19 2019-05-07 天津理工大学 A kind of Slope Stability Evaluation method
CN112526104A (en) * 2020-11-06 2021-03-19 马鞍山矿山研究总院股份有限公司 Slope stability monitoring and early warning method, system and medium

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102930348A (en) * 2012-10-19 2013-02-13 广东电网公司电力科学研究院 Assessment method for rainstorm disaster risks of sectional power transmission line pole-tower foundation slopes
CN103268420A (en) * 2013-05-24 2013-08-28 河海大学 Method for evaluating risks of high rock slope
KR101580062B1 (en) * 2014-12-05 2015-12-28 연세대학교 산학협력단 Method and system for real-time prediction and warning of landslides
US20190051146A1 (en) * 2017-08-09 2019-02-14 Institute Of Mountain Hazards And Environment, Chinese Academy Of Sciences Three-dimensional multi-point multi-index early warning method for risk at power grid tower in landslide section
CN109726902A (en) * 2018-12-19 2019-05-07 天津理工大学 A kind of Slope Stability Evaluation method
CN112526104A (en) * 2020-11-06 2021-03-19 马鞍山矿山研究总院股份有限公司 Slope stability monitoring and early warning method, system and medium

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘乃盛;李慧慧;廉斌;: "运河工程坝坡失稳模糊风险评估研究", 人民黄河, no. 02 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115100819A (en) * 2022-05-23 2022-09-23 深圳市北斗云信息技术有限公司 Landslide hazard early warning method and device based on big data analysis and electronic equipment
CN115100819B (en) * 2022-05-23 2023-12-01 深圳市北斗云信息技术有限公司 Landslide hazard early warning method and device based on big data analysis and electronic equipment

Also Published As

Publication number Publication date
CN113704849B (en) 2024-03-01

Similar Documents

Publication Publication Date Title
Stephenson et al. A new approach to flood vulnerability assessment for historic buildings in England
CN1922508B (en) Method and system for automated location dependent probabilistic tropical cyclone forecast
CN108960599B (en) Power transmission line rainstorm disaster refined prediction method and system based on inversion algorithm
CN111260187A (en) Intelligent mine geological environment information evaluation system and evaluation method
CN111445156B (en) Bias tunnel construction safety evaluation method based on variable weight fuzzy comprehensive evaluation
CN102629294A (en) Probability evaluation method of failure caused by typhoon to power transmission line
CN112132470A (en) Landslide susceptibility assessment method based on weighted information quantity method
CN112966722A (en) Regional landslide susceptibility prediction method based on semi-supervised random forest model
CN114442198B (en) Forest fire weather grade forecasting method based on weighting algorithm
CN111144637A (en) Regional power grid geological disaster forecasting model construction method based on machine learning
CN107403004B (en) Remote-measuring rainfall site suspicious numerical inspection method based on terrain data
CN113704849B (en) Evaluation method for power transmission and transformation project slope instability risk degree
CN113723824A (en) Rainstorm disaster risk assessment method
CN113987813A (en) Landslide sensitivity mapping model based on multi-parameter decision and analytic hierarchy process
CN114861451A (en) Typhoon-resistant design method for tower in coastal region
CN115564288A (en) Power transmission line geological disaster risk assessment method and system based on data fusion
CN116434475A (en) Landslide early warning method based on rainfall intensity and earth surface displacement change rate
CN116842435A (en) SWAP-based field drought three-dimensional measurement and time-space evolution visualization method
CN111027881A (en) Transmission tower natural disaster loss evaluation method based on improved grey correlation degree
Negasa GIS-based irrigation potential assessment for surface irrigation: The case of birbir River Watershed, oromia, Ethiopia
CN114093133A (en) Regional geological disaster weather forecast early warning method
CN116910041B (en) Daily correction method for remote sensing precipitation product based on scale analysis
CN113128811A (en) Power grid system geological disaster risk assessment method and system based on strong precipitation
CN116523291A (en) Disk buckle type full framing risk assessment method based on physical element extension
CN113516418A (en) Apple planting area drought disaster risk assessment method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant