CN111462450A - Mountain torrent early warning method considering rainfall spatial heterogeneity - Google Patents

Mountain torrent early warning method considering rainfall spatial heterogeneity Download PDF

Info

Publication number
CN111462450A
CN111462450A CN202010047148.6A CN202010047148A CN111462450A CN 111462450 A CN111462450 A CN 111462450A CN 202010047148 A CN202010047148 A CN 202010047148A CN 111462450 A CN111462450 A CN 111462450A
Authority
CN
China
Prior art keywords
rainfall
composite
early warning
warned
drainage basin
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
CN202010047148.6A
Other languages
Chinese (zh)
Other versions
CN111462450B (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.)
Huazhong University of Science and Technology
Bureau of Hydrology Changjiang Water Resources Commission
Original Assignee
Huazhong University of Science and Technology
Bureau of Hydrology Changjiang Water Resources Commission
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 Huazhong University of Science and Technology, Bureau of Hydrology Changjiang Water Resources Commission filed Critical Huazhong University of Science and Technology
Priority to CN202010047148.6A priority Critical patent/CN111462450B/en
Publication of CN111462450A publication Critical patent/CN111462450A/en
Application granted granted Critical
Publication of CN111462450B publication Critical patent/CN111462450B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/02Alarms for ensuring the safety of persons
    • G08B21/10Alarms for ensuring the safety of persons responsive to calamitous events, e.g. tornados or earthquakes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F23/00Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/18Status alarms
    • G08B21/182Level alarms, e.g. alarms responsive to variables exceeding a threshold
    • 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
    • Y02A50/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather

Abstract

The invention discloses a mountain torrent early warning method considering rainfall spatial heterogeneity, belonging to the field of hydrologic prediction and comprising the following steps: acquiring soil saturation and composite rainfall in a set time period before occurrence of multi-field historical flood peaks of a drainage basin to be pre-warned; the composite rainfall gives consideration to the rainfall and the rainfall spatial distribution; constructing a mountain torrent early warning model according to the soil saturation and the composite rainfall; the method comprises the steps of collecting the soil saturation of the drainage basin to be pre-warned, obtaining a critical composite rainfall index corresponding to the current soil saturation according to a mountain torrent pre-warning model, and judging whether the rainfall-generated flood is the super-warning flood or not by comparing the current composite rainfall of the drainage basin to be pre-warned with the critical composite rainfall index. The method can specifically quantify the unevenness of rainfall spatial distribution and the distribution condition of rainstorm centers, and the calculated critical composite rainfall index takes the rainfall into account as well as the rainfall spatial distribution, so that the accuracy of mountain torrents early warning can be improved.

Description

Mountain torrent early warning method considering rainfall spatial heterogeneity
Technical Field
The invention belongs to the technical field of hydrologic forecasting, and particularly relates to a mountain torrent early warning method considering rainfall spatial heterogeneity.
Background
Under the influence of global warming, the local heavy rainfall frequency in the hilly area is increased, and along with the aggravation of human activities in the hilly area, the mountain flood disaster gradually becomes a short board for flood prevention and disaster reduction in China. Therefore, national mountain torrent disaster investigation and evaluation work is carried out in 2010-2016 in China, great progress is made in key technologies of mountain torrent disaster prevention such as mountain torrent disaster risk zoning, forecasting and early warning, and a mountain torrent disaster prevention system is initially built. On one hand, however, the prevention and treatment amount of mountain torrent disasters in China is wide, and the treatment difficulty is high; on the other hand, the existing defense technology still has some 'neck clamping' problems to be broken through urgently, such as accuracy and timeliness of mountain torrents early warning and the like.
The early warning of the mountain torrents can be generally realized by establishing a high-precision flood forecasting model and carrying out early warning according to the fact that whether the forecast flood peak exceeds the warning flow of a river channel or not. Another feasible way is to establish early warning indexes closely related to mountain torrents and judge whether to perform early warning according to the early warning indexes. The critical rainfall is one of the commonly used early warning indexes, and is directly related to the rainfall, and is also influenced by the water content of soil, the heterogeneity of rainfall space and the like. The spatial heterogeneity of rainfall means that rainfall exhibits different distributions in different spatial regions of the flow field.
Quantitative indexes of rainfall spatial heterogeneity comprise a variation coefficient, an uneven coefficient and the like; the variation coefficient reflects the spatial discrete degree of rainfall data, and the non-uniform coefficient reflects the point-surface reduction degree of rainfall. Although the indexes can reflect the spatial difference of rainfall, the distribution condition of a rainstorm center cannot be reflected, the information contained in the indexes is not comprehensive enough, and the critical rainfall cannot be directly calculated, so that the accuracy of mountain torrent early warning is influenced.
Disclosure of Invention
Aiming at the defects or improvement requirements of the prior art, the invention provides a mountain torrent early warning method considering rainfall spatial heterogeneity, aiming at accurately quantifying the unevenness of rainfall spatial distribution and the distribution situation of rainstorm centers, and calculating a critical composite rainfall index which can consider the rainfall and the rainfall spatial distribution situation so as to improve the accuracy of mountain torrent early warning.
In order to achieve the purpose, the invention provides a mountain torrent early warning method considering rainfall spatial heterogeneity, which comprises the following steps:
s1: acquiring soil saturation and composite rainfall in a set time period before occurrence of multi-field historical flood peaks of a drainage basin to be pre-warned; the composite rainfall is calculated according to the surface average rainfall and the rainfall spatial heterogeneous index; the rainfall spatial heterogeneity index is obtained by calculating the average rainfall of the surface of the watershed to be pre-warned, the relative area of the equal-flow time surface and the average rainfall of the equal-flow time surface in a set time period;
s2: constructing a mountain torrent early warning model according to the soil saturation and the composite rainfall to obtain critical composite rainfall indexes corresponding to different soil saturations;
s3: the method comprises the steps of collecting the soil saturation of the drainage basin to be pre-warned, obtaining a critical composite rainfall index corresponding to the current soil saturation according to a mountain torrent pre-warning model, and judging whether the rainfall-generated flood is the super-warning flood or not by comparing the current composite rainfall of the drainage basin to be pre-warned with the critical composite rainfall index.
Further, the calculation method of the rainfall spatial heterogeneity index specifically comprises the following steps:
01: calculating the area of a Thiessen polygon corresponding to each rainfall site in the area to be early-warned, and calculating the average rainfall of the drainage basin surface to be early-warned according to the area of the Thiessen polygons corresponding to each rainfall site and the historical rainfall;
02: calculating the relative area and the average rainfall of each equal-flow time surface of the area to be pre-warned;
03: calculating the rainfall spatial heterogeneity index according to the following formula:
Figure BDA0002369834030000021
wherein SHIP represents rainfall spatial heterogeneity index, RkThe surface average rainfall of the k-th equal-flow surface is shown.
Further, the calculation formula of the average rainfall of the drainage basin surface to be pre-warned is as follows:
Figure BDA0002369834030000031
wherein the content of the first and second substances,
Figure BDA0002369834030000032
for the average rainfall of the drainage basin surface to be pre-warned, lambdaiThe proportion of the area of the Thiessen polygon corresponding to the ith rainfall station to the total area of the drainage basin is represented, wherein i is 1,2 …, n represents the total number of the rainfall stations in the drainage basin to be pre-warned; piIndicating the historical rainfall of the ith rainfall site.
Further, the relative area calculation formula of the equal flow time surface is as follows:
Figure BDA0002369834030000033
wherein f isjDenotes the relative area of the jth equal flow time plane, j is 1,2, …, m, m denotes the total number of equal flow time planes, ri,jAnd the proportion of the area of the Thiessen polygon corresponding to the rainfall station i contained in the jth equal-flow time plane to the total area of the drainage basin is shown.
Further, the surface average rainfall of the equal-flow time surface is calculated by the following formula:
Figure BDA0002369834030000034
wherein R isjThe surface average rainfall of the jth equal-flow surface is shown.
Further, the soil saturation S is calculated by the formula:
Figure BDA0002369834030000035
wherein, PaInfluence of rainfall in the early stages of flood peak occurrence, WmThe maximum soil moisture content in the drainage basin.
Further, the compound rainfall PsComprises the following steps:
Figure BDA0002369834030000036
further, the step S2 of constructing the mountain torrent early warning model according to the above indexes specifically includes the following steps:
01. data set D { (x)l,yl) 1,2, …, L } as input, yl1 as output; wherein x isl=[Ps,S],PsThe composite rainfall is represented by y which is 1 and represents an overtemperature, y which is 0 and represents no overtemperature, and L represents the number of actually measured historical peaks;
02. based on the classification principle of the support vector machine, the following objective function is obtained to enable the classification interval of the super-alarm flood and the non-super-alarm flood to be maximum:
Figure BDA0002369834030000041
wherein w ═ w1,w2]TOptimal hyperplane w for maximizing classification interval of hypersonic floods and non-hypersonic floodsTxl+ b is a normal vector of 0, b is a displacement term; c>0 is a penalty parameter, so as to make misclassification pointξ in as small a number as possiblelThe relaxation variable is more than or equal to 0, so that the influence of outliers on the classification model is reduced;
03. solving the objective function by adopting a Lagrange multiplier method and a sequence minimum optimization algorithm to obtain an optimal classification hyperplane wTxlThe + b is the parameter w and b of 0, and the critical composite rainfall index P corresponding to different soil saturation degrees is obtainedcComprises the following steps:
Figure BDA0002369834030000042
in general, the above technical solutions contemplated by the present invention can achieve the following advantageous effects compared to the prior art.
Compared with the traditional mountain torrent early warning model, the method and the system have the advantages that the influence of rainfall spatial heterogeneity on the critical rainfall of the mountain torrent is considered, the unevenness of rainfall spatial distribution and the distribution situation of a rainstorm center can be accurately quantized, the critical composite rainfall index which can consider the rainfall and the rainfall spatial distribution situation is obtained, and accordingly the accuracy of flood early warning is effectively improved.
Drawings
Fig. 1 is a flowchart of a mountain torrent early warning method considering rainfall spatial heterogeneity according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a rainfall spatial distribution curve;
fig. 3(a) is a schematic diagram of flood warning results without considering spatial heterogeneity;
fig. 3(b) is a schematic diagram of flood warning results considering spatial heterogeneity.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. In addition, the technical features involved in the embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
As shown in fig. 1, the present invention provides a mountain torrent early warning method considering rainfall spatial heterogeneity, including:
s1: acquiring soil saturation and composite rainfall in a set time period before occurrence of multi-field historical flood peaks of a drainage basin to be pre-warned; the composite rainfall is calculated according to the surface average rainfall and the rainfall spatial heterogeneous index; the rainfall spatial heterogeneity index is obtained by calculating the average rainfall of the surface of the watershed to be pre-warned, the relative area of the equal-flow time surface and the average rainfall of the equal-flow time surface in a set time period;
specifically, in the embodiment of the present invention, in order to obtain the rainfall spatial distribution of the multi-field flood peak which occurs in the history, the soil saturation and the composite rainfall corresponding to 3 hours before each flood peak occurs and with the maximum continuous rainfall amount are calculated, and in a specific application, data corresponding to 6 hours, 12 hours, and 24 hours before the flood peak occurs and with the maximum continuous rainfall amount can also be obtained.
Wherein, the calculation formula of the soil saturation S is as follows:
Figure BDA0002369834030000051
wherein, PaInfluence of rainfall in the early stages of flood peak occurrence, WmThe maximum soil water content of the drainage basin;
wherein the compound rainfall PsThe calculation formula is as follows:
Figure BDA0002369834030000052
wherein, SHIP is rainfall spatial heterogeneity index,
Figure BDA0002369834030000053
the average rainfall of the drainage basin surface is obtained;
the method for calculating the rainfall spatial heterogeneity index comprises the following steps:
01: drawing a Thiessen polygon corresponding to each rainfall site in the area to be early-warned to obtain the area of the Thiessen polygon corresponding to each rainfall site, and calculating the average rainfall of the drainage basin surface to be early-warned according to the area of the Thiessen polygon corresponding to each rainfall site and the historical rainfall, wherein the specific calculation formula is as follows:
Figure BDA0002369834030000061
wherein the content of the first and second substances,
Figure BDA0002369834030000062
for the average rainfall of the drainage basin surface to be pre-warned, lambdaiThe proportion of the area of the Thiessen polygon corresponding to the ith rainfall station to the total area of the drainage basin is represented, wherein i is 1,2 …, n represents the total number of the rainfall stations in the drainage basin to be pre-warned; piRepresenting the historical rainfall of the ith rainfall station;
02: drawing an equal flow time line of a region to be pre-warned to obtain a plurality of equal flow time surfaces, and calculating the relative area and surface average rainfall of each equal flow time surface; the relative area calculation formula of the equal flow time surface is as follows:
Figure BDA0002369834030000063
wherein f isjDenotes the relative area of the jth equal flow time plane, j is 1,2, …, m, m denotes the total number of equal flow time planes, ri,jAnd the proportion of the area of the Thiessen polygon corresponding to the rainfall station i contained in the jth equal-flow time plane to the total area of the drainage basin is shown.
The surface average rainfall of the surface during equal flow is calculated according to the formula:
Figure BDA0002369834030000064
wherein R isjThe surface average rainfall of the jth equal-flow surface is shown.
03: calculating the rainfall spatial heterogeneity index according to the following formula:
Figure BDA0002369834030000065
wherein R iskSurface average of k-th equal flow surfaceAnd (4) rainfall. With equal flow time surface relative area fjThe accumulated value of (A) is an abscissa, and the surface average rainfall R of the corresponding equal-flow time surface isjThe accumulated value of (2) is a vertical coordinate, a rainfall spatial distribution curve is drawn, the result shown in fig. 2 is obtained, when the SHIP is positive, the rainstorm center is located at the downstream, such as a curve A in fig. 2; when SHIP is negative, it indicates that the storm centre is upstream, as shown by curve B in FIG. 2.
S2: constructing a mountain torrent early warning model according to the indexes, and calculating dynamic critical composite rainfall indexes under different soil saturation degrees;
specifically, S2 includes the steps of:
2.1, dividing actually measured L flood peaks into two types of super-alarm and non-super-alarm according to the warning flow, respectively representing that y is 1 and y is 0, and combining two indexes [ P ] corresponding to each flood peaks,S]Construct dataset D { (x)l,yl)|l=1,2,…,L},xl=[Ps,S]As an input; y islAs output {0,1 }.
2.2 based on the classification principle of the support vector machine, there is an optimal classification hyperplane
wTxl+b=0
The classification interval (2/| w | |) of the flood with the excess flood and the flood without the excess flood is maximized. Wherein w ═ w1,w2]TB is a displacement term for optimally classifying the normal vector of the hyperplane. Thus, the problem of finding the optimal hyperplane translates into the following optimization problem:
Figure BDA0002369834030000071
in the formula, C>0 is a penalty parameter, the number of misclassification points can be reduced as much as possible by introducing C ξ l0 or more as a relaxation variable ξlThe introduction of (2) can reduce the influence of outliers on the classification model.
2.3, obtaining the dual problem of the optimization problem by using a Lagrange multiplier method, solving by using a sequence minimum optimization algorithm to obtain parameters w and b of the optimal hyperplane equation, wherein the hyperplane equation is a critical early warning plane equation, and thus obtaining the corresponding relation between the critical composite rainfall index and the soil saturation, namely the corresponding relation between the critical composite rainfall index and the soil saturation
Figure BDA0002369834030000072
Wherein, PcIs a critical composite rainfall index.
S3: collecting the soil saturation of the drainage basin to be pre-warned, obtaining a critical composite rainfall index corresponding to the current soil saturation according to the mountain torrent pre-warning model, and judging whether the flood generated by rainfall is over-warning flood or not by comparing the current composite rainfall of the drainage basin to be pre-warned with the critical rainfall index.
When the current average rainfall of the watershed to be pre-warned is greater than the critical composite rainfall index, judging that the watershed to be pre-warned is over-warning flood; otherwise, the flood is not over-police.
In order to verify the effectiveness of the method, the embodiment of the invention takes the river basin of Hanjiang river as an example to carry out mountain torrent early warning. Specifically, hydrological data of each rainfall site of the constant river basin in 1980-2017 are obtained, 30 floods of the basin are selected for model training, and 350m is measured3The flood of more than s is taken as early warning flood, and the soil saturation and the composite rainfall corresponding to 3 hours with the maximum continuous rainfall within 24 hours before the flood peak of each field occurs are only calculated in the embodiment of the invention.
The critical composite rainfall indexes corresponding to different soil saturation degrees of the constant river basin are obtained by the method, and are shown in table 1:
TABLE 1
Figure BDA0002369834030000081
In consideration of the spatial distribution characteristics of rainfall and watershed topographic features, the heavy rainfall process near the downstream of the rainstorm center often causes the downstream river channel to have steep rising flood, so that the mountain torrents rise very quickly, while the rainfall process near the upstream of the rainstorm center needs to pass through the long regulation and storage effect of the watershed, so that the mountain torrents rise slowly. In the invention, the influence of spatial non-uniformity of rainfall on the flood in the drainage basin is considered, critical composite rainfall is used as an early warning index, the composite rainfall depends on the product of (1+ SHIP) and surface average rainfall, when the rainfall center position is at the upstream of the drainage basin, the SHIP value is a negative value, and the drainage basin can reach the early warning value of the flood only by larger rainfall; when the rainfall center is positioned at the downstream of the drainage basin, the SHIP value is a positive value, and the drainage basin can reach the early warning value of flood only by small rainfall. The critical composite rainfall index is more consistent with the actual result because the spatial heterogeneity of rainfall and the surface average rainfall are simultaneously considered.
The mountain torrent early warning result considering spatial heterogeneity is shown in fig. 3(a), and the early warning accuracy reaches 80%; the mountain torrent early warning result without considering spatial heterogeneity is shown in fig. 3(b), the early warning accuracy is only 66.67%, and it can be seen that the method of the invention can effectively improve the accuracy of flood early warning.
It will be understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, and that any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (8)

1. A mountain torrent early warning method considering rainfall spatial heterogeneity is characterized by comprising the following steps:
s1: acquiring soil saturation and composite rainfall in a set time period before occurrence of multi-field historical flood peaks of a drainage basin to be pre-warned; the composite rainfall is calculated according to the surface average rainfall and the rainfall spatial heterogeneous index; the rainfall spatial heterogeneity index is obtained by calculation according to the average rainfall of the surface of the watershed to be pre-warned in a set time period, the relative area of the equal-flow time surface and the average rainfall of the equal-flow time surface;
s2: constructing a mountain torrent early warning model according to the soil saturation and the composite rainfall to obtain critical composite rainfall indexes corresponding to different soil saturations;
s3: collecting the soil saturation of the drainage basin to be pre-warned, obtaining a critical composite rainfall index corresponding to the current soil saturation according to the mountain torrent pre-warning model, and judging whether the rainfall-generated flood is the super-warning flood or not by comparing the current composite rainfall of the drainage basin to be pre-warned with the critical composite rainfall index.
2. The mountain torrent early warning method considering rainfall spatial heterogeneity as claimed in claim 1, wherein the calculation method of the rainfall spatial heterogeneity index specifically comprises:
01: calculating the area of a Thiessen polygon corresponding to each rainfall site in the area to be early-warned, and calculating the average rainfall of the drainage basin surface to be early-warned according to the area of the Thiessen polygons corresponding to each rainfall site and the historical rainfall;
02: calculating the relative area and the average rainfall of each equal-flow time surface of the area to be pre-warned;
03: calculating the rainfall spatial heterogeneity index according to the following formula:
Figure FDA0002369834020000011
wherein, SHIP represents a rainfall spatial heterogeneity index,
Figure FDA0002369834020000012
for average rainfall of the watershed surface to be warned, fjDenotes the relative area of the jth equal flow time surface, m denotes the total number of equal flow time surfaces, RkThe surface average rainfall of the k-th equal-flow surface is shown.
3. The mountain torrent early warning method considering rainfall spatial heterogeneity as claimed in claim 2, wherein the calculation formula of the average rainfall of the drainage basin to be early warned is as follows:
Figure FDA0002369834020000021
wherein λ isiThe proportion of the area of the Thiessen polygon corresponding to the ith rainfall station to the total area of the drainage basin is represented, wherein i is 1,2 …, n represents the total number of the rainfall stations in the drainage basin to be pre-warned; piIndicating the historical rainfall of the ith rainfall site.
4. The mountain torrent early warning method considering rainfall spatial heterogeneity as claimed in claim 2 or 3, wherein the formula for calculating the relative area of the equal flow time surface is as follows:
Figure FDA0002369834020000022
wherein r isi,jAnd the proportion of the area of the Thiessen polygon corresponding to the rainfall station i contained in the jth equal-flow time plane to the total area of the drainage basin is shown.
5. The mountain torrent early warning method considering rainfall spatial heterogeneity as claimed in any one of claims 2 to 4, wherein a surface average rainfall calculation formula of the equal flow time surface is as follows:
Figure FDA0002369834020000023
wherein R isjThe surface average rainfall of the jth equal-flow surface is shown.
6. The mountain torrent early warning method considering rainfall spatial heterogeneity as claimed in any one of claims 1 to 5, wherein a soil saturation S calculation formula is:
Figure FDA0002369834020000024
wherein, PaInfluence of rainfall in the early stages of flood peak occurrence, WmThe maximum soil moisture content in the drainage basin.
7. The mountain torrent early warning method considering rainfall spatial heterogeneity as claimed in any one of claims 1 to 6, wherein composite rainfall PsComprises the following steps:
Figure FDA0002369834020000025
8. the method as claimed in any one of claims 1 to 7, wherein the step of constructing the mountain torrent early warning model according to the above indexes in step S2 comprises the following steps:
01. data set D { (x)l,yl) 1, 2.., L } as input, yl1 as output; wherein x isl=[Ps,S],PsThe composite rainfall is represented by y which is 1 and represents an overtemperature, y which is 0 and represents no overtemperature, and L represents the number of actually measured historical peaks;
02. based on the classification principle of the support vector machine, the following objective function is obtained to enable the classification interval of the super-alarm flood and the non-super-alarm flood to be maximum:
Figure FDA0002369834020000031
wherein w ═ w1,w2]TOptimal hyperplane w for maximizing classification interval of hypersonic floods and non-hypersonic floodsTxl+ b is a normal vector of 0, b is a displacement term; c>0 is penalty parameter to minimize the number of misclassification points ξlThe relaxation variable is more than or equal to 0, so that the influence of outliers on the classification model is reduced;
03. solving the objective function by adopting a Lagrange multiplier method and a sequence minimum optimization algorithm to obtain an optimal classification hyperplane wTxlThe + b is the parameter w and b of 0, and the critical composite rainfall index P corresponding to different soil saturation degrees is obtainedcComprises the following steps:
Figure FDA0002369834020000032
CN202010047148.6A 2020-01-16 2020-01-16 Mountain torrent early warning method considering rainfall spatial heterogeneity Active CN111462450B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010047148.6A CN111462450B (en) 2020-01-16 2020-01-16 Mountain torrent early warning method considering rainfall spatial heterogeneity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010047148.6A CN111462450B (en) 2020-01-16 2020-01-16 Mountain torrent early warning method considering rainfall spatial heterogeneity

Publications (2)

Publication Number Publication Date
CN111462450A true CN111462450A (en) 2020-07-28
CN111462450B CN111462450B (en) 2021-06-11

Family

ID=71679936

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010047148.6A Active CN111462450B (en) 2020-01-16 2020-01-16 Mountain torrent early warning method considering rainfall spatial heterogeneity

Country Status (1)

Country Link
CN (1) CN111462450B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112067043A (en) * 2020-08-14 2020-12-11 常州机电职业技术学院 Defective degree detecting system of timber structure ancient building
CN112216061A (en) * 2020-09-29 2021-01-12 浪潮云信息技术股份公司 Rainwater condition monitoring and early warning method and system
CN112396297A (en) * 2020-11-03 2021-02-23 华中科技大学 Method and system for analyzing encounter time and magnitude occurrence rule in flood process
CN113687448A (en) * 2021-08-26 2021-11-23 中水珠江规划勘测设计有限公司 Precipitation center position and change determination method and device and electronic equipment

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106529176A (en) * 2016-11-11 2017-03-22 中国水利水电科学研究院 Dual-core dual-drive flood forecast method
CN108597189A (en) * 2018-04-24 2018-09-28 河海大学 Small watershed geological disaster and flood warning method in distribution based on Critical Rainfall
CN109633790A (en) * 2019-01-18 2019-04-16 三峡大学 The method of sub-basin rainfall spatial and temporal distributions is determined in natural basin partitioning
CN109920213A (en) * 2019-03-13 2019-06-21 河海大学 The method that Critical Rainfall based on rainfall rainfall process carries out real-time mountain torrents early warning
KR102009373B1 (en) * 2019-05-22 2019-08-12 (주)현이엔씨 Estimation method of flood discharge for varying rainfall duration

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106529176A (en) * 2016-11-11 2017-03-22 中国水利水电科学研究院 Dual-core dual-drive flood forecast method
CN108597189A (en) * 2018-04-24 2018-09-28 河海大学 Small watershed geological disaster and flood warning method in distribution based on Critical Rainfall
CN109633790A (en) * 2019-01-18 2019-04-16 三峡大学 The method of sub-basin rainfall spatial and temporal distributions is determined in natural basin partitioning
CN109920213A (en) * 2019-03-13 2019-06-21 河海大学 The method that Critical Rainfall based on rainfall rainfall process carries out real-time mountain torrents early warning
KR102009373B1 (en) * 2019-05-22 2019-08-12 (주)현이엔씨 Estimation method of flood discharge for varying rainfall duration

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
罗倩等: "降雨量空间分布对山洪临界雨量的影响", 《人民长江》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112067043A (en) * 2020-08-14 2020-12-11 常州机电职业技术学院 Defective degree detecting system of timber structure ancient building
CN112216061A (en) * 2020-09-29 2021-01-12 浪潮云信息技术股份公司 Rainwater condition monitoring and early warning method and system
CN112396297A (en) * 2020-11-03 2021-02-23 华中科技大学 Method and system for analyzing encounter time and magnitude occurrence rule in flood process
CN112396297B (en) * 2020-11-03 2021-06-29 华中科技大学 Method and system for analyzing encounter time and magnitude occurrence rule in flood process
CN113687448A (en) * 2021-08-26 2021-11-23 中水珠江规划勘测设计有限公司 Precipitation center position and change determination method and device and electronic equipment
CN113687448B (en) * 2021-08-26 2023-11-10 中水珠江规划勘测设计有限公司 Precipitation center position and variation determining method and device thereof and electronic equipment

Also Published As

Publication number Publication date
CN111462450B (en) 2021-06-11

Similar Documents

Publication Publication Date Title
CN111462450B (en) Mountain torrent early warning method considering rainfall spatial heterogeneity
CN112070286B (en) Precipitation forecast and early warning system for complex terrain river basin
CN103729550B (en) Multiple-model integration Flood Forecasting Method based on propagation time cluster analysis
CN111582755A (en) Mountain torrent disaster comprehensive risk dynamic assessment method based on multi-dimensional set information
CN102509173B (en) A kind of based on markovian power system load Accurate Prediction method
CN110047291A (en) A kind of Short-time Traffic Flow Forecasting Methods considering diffusion process
CN108597022A (en) A kind of method of small watershed inland river road width in estimation
CN113705931B (en) Method for predicting runoff elements by using K nearest neighbor method
CN112633645B (en) Social and economic benefit accounting method for water resource mining and efficient utilization effects of river source arid region
CN115115262B (en) Flood risk disaster assessment method
CN112686426A (en) Incoming water quantity early warning method and system based on hydropower station basin key points
CN110543660B (en) Low-impact development simulation method, system and related device
CN114723283A (en) Ecological bearing capacity remote sensing evaluation method and device for urban group
CN116596303A (en) Drought risk assessment and zoning method, system, medium, equipment and terminal
CN106022507B (en) Optimization method and system for farmland connected piece remediation
CN110442988A (en) A kind of city overland runoff based on cellular automata flows to calculation method and device
CN108615092A (en) A method of the sewage treatment plant inflow amount prediction based on exponential smoothing model
CN113128811A (en) Power grid system geological disaster risk assessment method and system based on strong precipitation
CN112329969A (en) Building intelligent engineering investment prediction method based on support vector machine
CN115358587A (en) Regional multi-department collaborative infrastructure planning method and system
CN115453664A (en) Rainfall runoff forecasting method suitable for data-free areas
CN112528563B (en) Urban waterlogging early warning method based on SVM algorithm
Teng et al. Early warning index of flash flood disaster: a case study of Shuyuan watershed in Qufu City
CN111815121B (en) Operation management evaluation method of drainage deep tunnel system based on SWMM model
CN116628411B (en) High-precision flow online monitoring intelligent method based on full-sense fusion

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