[go: up one dir, main page]

CN113961623A - Rainfall characteristic index statistical method based on short-duration data - Google Patents

Rainfall characteristic index statistical method based on short-duration data Download PDF

Info

Publication number
CN113961623A
CN113961623A CN202111225876.2A CN202111225876A CN113961623A CN 113961623 A CN113961623 A CN 113961623A CN 202111225876 A CN202111225876 A CN 202111225876A CN 113961623 A CN113961623 A CN 113961623A
Authority
CN
China
Prior art keywords
rainfall
station
area
intensity
total
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.)
Pending
Application number
CN202111225876.2A
Other languages
Chinese (zh)
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.)
Yellow River Engineering Consulting Co Ltd
Original Assignee
Yellow River Engineering Consulting 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 Yellow River Engineering Consulting Co Ltd filed Critical Yellow River Engineering Consulting Co Ltd
Priority to CN202111225876.2A priority Critical patent/CN113961623A/en
Publication of CN113961623A publication Critical patent/CN113961623A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24—Querying
    • G06F16/245—Query processing
    • G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2462—Approximate or statistical queries
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/14—Rainfall or precipitation gauges
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24—Querying
    • G06F16/245—Query processing
    • G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2477—Temporal data queries
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29—Geographical information databases
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Probability & Statistics with Applications (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Fuzzy Systems (AREA)
  • Software Systems (AREA)
  • General Business, Economics & Management (AREA)
  • Mathematical Physics (AREA)
  • Marketing (AREA)
  • Environmental & Geological Engineering (AREA)
  • Computational Linguistics (AREA)
  • Development Economics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Hydrology & Water Resources (AREA)
  • Atmospheric Sciences (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Remote Sensing (AREA)
  • Educational Administration (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Sewage (AREA)

Abstract

The invention provides a rainfall characteristic index statistical method based on short-duration data, which comprises the following steps: collecting and calculating precipitation data of each rainfall station in the flow field, and simultaneously determining the control area of each rainfall station; calculating the rainfall intensity of each rainfall record according to the beginning and ending time of each rainfall excerpted data of each rainfall station in a preset month and the rainfall; screening rainfall records of a single rainfall station according to different rainfall intensity standards, accumulating the rainfall of all the rainfall records reaching specific rainfall intensity, and multiplying the rainfall records by the control area of the rainfall station to obtain a first total rainfall amount reaching the specific rainfall intensity corresponding to the single rainfall station; and superposing the total rainfall amount reaching the specific rainfall intensity corresponding to each rainfall station in the river basin to obtain a second total rainfall amount reaching the specific rainfall intensity corresponding to the river basin. The method has the advantages of simple steps, clear concept, simple and convenient calculation and easy operation, and is an effective tool for researching rainfall change characteristics of the drainage basin.

Description

Rainfall characteristic index statistical method based on short-duration data
Technical Field
The invention relates to the technical field of hydrological weather, in particular to a rainfall characteristic index statistical method based on short-duration data.
Background
The rainfall change characteristic research is an important field of climate change research and is also an important means for analyzing the change reasons of runoff, flood and sediment in a drainage basin. For areas mainly based on the super-percolation runoff, the changes of runoff, flood and silt in the drainage basin are simultaneously influenced by rainfall and rainfall intensity, and comprehensive rainfall indexes capable of reflecting the rainfall and the rainfall intensity simultaneously are needed to identify the change characteristics and the change trend of rainfall. On the other hand, the reliable short-duration rainfall data observed by the current water conservancy department is generally published in the hydrological statistics yearbook and appears in the form of a rainfall extraction table. The precipitation extraction data only record flood season data generally, and are mainly concentrated in 6-9 months, wherein the data of 7-8 months are basically complete. The precipitation extraction data generally comprises three parts of starting time, stopping time and precipitation, wherein the minimum precision of the starting time and the stopping time is 1min, the minimum precision of the precipitation is 0.1mm, and the recorded time interval is different from 1min to 24 h. At present, rainfall change research is carried out on the basis of data extracted from rainfall, which takes hours as a scale, and the observation precision of partial rainfall is lost.
Aiming at the situations, the rainfall and the rainfall duration of each rainfall record are counted one by one on the basis of the short-duration rainfall extraction data, the observation time precision (the maximum precision is 1min) is fully utilized, and a comprehensive rainfall characteristic index capable of reflecting the rainfall intensity and the rainfall change situation of the drainage basin at the same time is constructed. The method has the advantages of simple steps, clear concept, simple and convenient calculation and easy operation, is an effective tool for researching rainfall change characteristics of the drainage basin, and particularly has positive significance for researching rainfall change in the area with the super-seepage as a main area.
Disclosure of Invention
The invention provides a rainfall characteristic index statistical method based on short-duration data, which is used for screening the single rainfall total amount of each rainfall station reaching specific rainfall intensity in a short-duration month according to the result of rainfall station measurement so as to obtain the total rainfall amount reaching the specific rainfall amount in a basin.
The invention provides a rainfall characteristic index statistical method based on short-duration data, which comprises the following steps:
step 1, collecting and calculating precipitation data of each rainfall station in a flow domain, and simultaneously determining the control area of each rainfall station;
step 2, calculating the rainfall intensity of each rainfall record according to the beginning and ending time of each rainfall excerpted data of each rainfall station in a preset month and the rainfall;
step 3, screening the rainfall records of the single rainfall station according to different rainfall intensity standards, accumulating the rainfall of all the rainfall records reaching the specific rainfall intensity, and multiplying the rainfall records by the control area of the rainfall station to obtain a first total rainfall amount reaching the specific rainfall intensity corresponding to the single rainfall station;
and 4, overlapping the total rainfall amount reaching the specific rainfall intensity corresponding to each rainfall station in the river basin to obtain a second total rainfall amount reaching the specific rainfall intensity corresponding to the river basin.
In a possible implementation manner, step 4, overlapping the total rainfall amounts reaching the specific rainfall intensity corresponding to each rainfall station in the drainage basin to obtain a second total rainfall amount reaching the specific rainfall intensity corresponding to the drainage basin, including:
calculating a second total rainfall amount corresponding to the drainage basin reaching a specific rainfall intensity according to the following formula:
Figure BDA0003313983530000021
in the formula: wrA second total rainfall amount corresponding to the drainage basin and reaching a specific rainfall intensity (r mm per minute) or more, wherein the unit is billionth of cubic meters; j is the jth rainfall station with the intensity more than r millimeters per minute in the drainage basin; m is the total number of the rainfall stations reaching the strength of more than r millimeters per minute in the drainage basin; n is the total rainfall record number of the rainfall station corresponding to the intensity of more than r millimeters per minute;
Figure BDA0003313983530000022
recording corresponding rainfall capacity in millimeter for the ith rainfall with the rainfall intensity of more than r millimeter per minute in the preset month flow domain; a is the control area of the corresponding rainfall station, and the unit is square kilometer,
Figure BDA0003313983530000023
and the j-th rainfall station is expressed to reach a first rainfall total amount corresponding to the intensity of more than r millimeters per minute in the preset month.
In a possible implementation manner, before collecting and calculating precipitation data of each rainfall station in the flow domain, step 1 further includes:
determining the current position of each rainfall station, constructing a historical weather model in a control area corresponding to the control area of each rainfall station according to the weather change condition of a preset month in a historical year, and simulating the rainfall distribution of the control area according to the historical weather model;
constructing regional terrain of a control region corresponding to each rainfall station, and carrying out terrain type division on the regional terrain;
carrying out region matching on the rainfall distribution and the divided terrain results one by one, and determining an independent rainfall measurement region and a dependent rainfall measurement region in each control region;
acquiring a distribution mode according to the regional attribute of the dependent rainfall measurement region, and distributing corresponding auxiliary rainfall stations according to the distribution mode, wherein the auxiliary rainfall stations are used for monitoring and recording the rainfall of the corresponding dependent rainfall measurement region;
the auxiliary area is in the control area, and the number of the auxiliary rainfall stations is at least one.
In a possible implementation manner, step 4, after obtaining a second total rainfall amount corresponding to the drainage basin reaching a specific rainfall intensity, further includes: and performing correction adjustment on the second total rainfall, which comprises the following steps:
marking a first rainfall station which is initially set;
acquiring a control area corresponding to the first rainfall station and taking the control area as an initial area Aj；
Acquiring monitoring information of each first rainfall station reaching a specific rainfall intensity in a preset month, wherein the monitoring information comprises: a first time point when the specific rainfall is reached, a second time set of continuous rainfall before the specific rainfall is reached, and a third time set of continuous rainfall after the specific rainfall is reached;
acquiring intra-domain weather change information of a sky area of the size corresponding to the initial area of each first rainfall station from a weather observation database corresponding to the preset month, and constructing a weather change model corresponding to a time period according to the intra-domain weather change information, wherein the intra-domain weather change information is related to cloud coverage area, cloud thickness and cloud moving speed;
according to the weather change model, the weather change layers corresponding to the specific rainfall intensity related to the first rainfall station at different time points are divided and output, and are placed in an overlapping mode to obtain a stability maintaining area and a residual area of the first rainfall station, wherein the areas of the stability maintaining area and the residual area form corresponding initial areas;
meanwhile, according to overlapped color results after overlapping, performing region division on the remaining region to obtain a plurality of sub-regions, calling and monitoring the rainfall of the corresponding sub-regions based on the fact that a second rainfall station is arranged in the sub-regions in advance, and obtaining the optimized rainfall corresponding to the initial area by combining the rainfall of the stable region monitored by the corresponding first rainfall station;
based on the optimized rainfall, the second total rainfall is corrected and adjusted to obtain a third total rainfall, the third total rainfall is compared with the second total rainfall, and if the absolute value of the difference value of the third total rainfall and the second total rainfall is within a preset difference value range, rainfall data monitored by the first rainfall station is continuously used for estimating the total rainfall;
otherwise, estimating the total rainfall of the rainfall data monitored by the first rainfall station and the second rainfall station;
in a possible implementation manner, obtaining an optimized rainfall corresponding to an initial area includes:
Figure BDA0003313983530000041
wherein, YjRepresenting the optimized rainfall of the original corresponding jth initial area; g represents the total number of the sub-areas of the area corresponding to the jth initial area, a second rainfall station is started in each sub-area, and the total number of the second rainfall stations is also G; fg1,iIndicating that the recorded rainfall corresponding to the g1 th second rainfall station when the recorded area corresponding to the j initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute; sg1,iRepresenting the area of the corresponding g1 th sub-area when the recorded area corresponding to the jth initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute; fj,iRecording rainfall measured by a corresponding first rainfall station when the area corresponding to the jth initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute; sj,iRepresenting the area of the corresponding stable region when the recorded corresponding region of the jth initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute; deltag1,iPresentation noteRecording effective measurement coefficients corresponding to the rainfall station when the corresponding area of the jth initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute, wherein the value range is (0.8, 1); deltaj,iThe effective measurement coefficient of the corresponding rainfall station is represented when the recorded corresponding area of the jth initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute, and the value range is (0.8, 1);
wherein when is deltag1,iWhen the content is equal to 0, the content,
Figure BDA0003313983530000042
when deltaj,iWhen the content is equal to 0, the content,
Figure BDA0003313983530000051
in one possible implementation manner, performing a correction adjustment on the second total rainfall amount based on the optimized rainfall amount to obtain a third total rainfall amount includes:
Figure BDA0003313983530000052
wherein, WsIs the third total amount of rainfall present,
Figure BDA0003313983530000053
the ratio of the area corresponding to the jth initial area to the second rainfall station when the ith area reaches the rainfall with the rainfall intensity of more than r millimeters per minute; when in use
Figure BDA0003313983530000054
When the ratio of (a) is greater than or equal to 1/2, a takes the value of 1; when in use
Figure BDA0003313983530000055
When the ratio of (a) is less than 1/2, a is 2.
In one possible implementation, the dimensionally stable region is inclusive of an independent rainfall measurement region, and the dependent rainfall measurement region is inclusive of a remaining region;
the second rainfall stations are obtained by screening in auxiliary rainfall stations arranged in the corresponding controlled area, and the number of the auxiliary rainfall stations is larger than that of the corresponding second rainfall stations.
In a possible implementation manner, before calculating the rainfall intensity of each rainfall record, the method further includes:
determining the rainfall time period of each rainfall station according to the starting time and the ending time of each rainfall excerpted data of each rainfall station in a preset month, and simultaneously acquiring weather data corresponding to each rainfall time period;
calculating a first matching degree P1 between the rainfall time length of the rainfall time period corresponding to the rainfall station and the estimated rainfall time length corresponding to the weather data according to the following formula;
Figure BDA0003313983530000056
wherein F1 represents the total rainfall time period in the preset month measured by the same rainfall station; p1f1A rainfall period representing the f1 th rainfall period measured at the same rainfall station; p2f1Indicating an estimated rainfall duration determined based on weather data when measuring the rainfall duration of the f1 th rainfall time period of the same rainfall station;
calculating a second matching degree P2 between the rainfall capacity of the rainfall time period corresponding to the rainfall station and the estimated rainfall capacity corresponding to the weather data according to the following formula;
Figure BDA0003313983530000061
wherein F1 represents the total rainfall time period in the preset month measured by the same rainfall station; y is1f1Represents the rainfall of the f1 th rainfall period measured by the same rainfall station; y is2f1Represents an estimated rainfall determined based on weather data when measuring the rainfall at the same rainfall station at the f1 th rainfall period; df1Representing weather-based numbersCorrecting parameters according to weather, wherein the value range is (0.1, 0.2);
when the first matching degree P1 is greater than the corresponding first preset matching degree and the second matching degree P2 is greater than the corresponding second preset matching degree, judging that the corresponding rainfall station is qualified, and controlling the qualified rainfall station to continue rainfall measurement;
otherwise, collecting the operation information of the unqualified rainfall station in a preset month, and analyzing the unqualified type of the unqualified rainfall station;
when the unqualified type is irrelevant to equipment upgrading, alarming and reminding are carried out, meanwhile, a standby rainfall station is started to carry out measurement work, meanwhile, a corresponding operation information is reversely analyzed based on a historical estimation model of the unqualified rainfall station, and the estimated rainfall of the unqualified rainfall station is estimated to replace a rainfall record in unqualified operation;
and when the unqualified type is related to equipment upgrading, upgrading the corresponding unqualified type of the corresponding unqualified rainfall station.
Compared with the prior art, the invention has the following beneficial effects:
1. the rainfall amount excerpted data with short duration is fully utilized, the rainfall amount and the rainfall duration of each rainfall record are counted one by one, the observation time precision is utilized to the maximum extent, and the comprehensive rainfall characteristic index capable of reflecting the rainfall intensity and the rainfall amount change condition of the drainage basin is constructed. The method has the advantages of simple steps, clear concept, simple and convenient calculation and easy operation, is an effective tool for researching rainfall change characteristics of the drainage basin, and particularly has positive significance for researching rainfall change in the area with the super-seepage as a main area.
2. According to the invention, the accurate rainfall of the sub-area set in the control area corresponding to each first rainfall station is obtained, the optimized rainfall of all the sub-areas in the control area is further determined according to the accurate rainfall, the rainfall of the control area under different weather conditions at the same time point is conveniently and accurately obtained, the total rainfall of the control area is further accurately obtained, finally, the second total rainfall is corrected and adjusted, the third total rainfall is effectively obtained, the subsequent effective comparison of the rainfall is convenient, the reasonable division of labor is carried out for the area size of the control area, and the measurement accuracy is further ensured.
3. The rainfall station is comprehensively determined whether corresponding rainfall stations have problems or not through the rainfall duration and the rainfall, wherein the rainfall stations are comprehensively determined through comparison of the calculated matching degree and the preset matching degree, the accuracy of the recorded rainfall in the initial process is guaranteed, the unqualified rainfall stations are regulated and processed in relevant unqualified types, the normal operation of the unqualified rainfall stations is guaranteed, the validity of measured data is further guaranteed, and an effective data basis is provided for the total rainfall of the basin obtained after the specific rainfall is achieved subsequently.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention and not to limit the invention. In the drawings:
FIG. 1 is a flow chart of a method for rainfall feature index statistics based on short-duration data according to an embodiment of the present invention;
FIG. 2 is a diagram illustrating an initial division result of a control area according to an embodiment of the present invention;
FIG. 3 is a diagram illustrating the result of the Steady region and the residual region in the embodiment of the present invention.
Detailed Description
The preferred embodiments of the present invention will be described in conjunction with the accompanying drawings, and it will be understood that they are described herein for the purpose of illustration and explanation and not limitation.
The invention provides a rainfall characteristic index statistical method based on short-duration data, as shown in figure 1, comprising the following steps:
step 1, collecting and calculating precipitation data of each rainfall station in a flow domain, and simultaneously determining the control area of each rainfall station;
step 2, calculating the rainfall intensity of each rainfall record according to the beginning and ending time of each rainfall excerpted data of each rainfall station in a preset month and the rainfall;
step 3, screening the rainfall records of the single rainfall station according to different rainfall intensity standards, accumulating the rainfall of all the rainfall records reaching the specific rainfall intensity, and multiplying the rainfall records by the control area of the rainfall station to obtain a first total rainfall amount reaching the specific rainfall intensity corresponding to the single rainfall station;
and 4, overlapping the total rainfall amount reaching the specific rainfall intensity corresponding to each rainfall station in the river basin to obtain a second total rainfall amount reaching the specific rainfall intensity corresponding to the river basin.
In this embodiment, the preset months may refer to months 7 and 8, and the thieson polygon method may be used for calculation in determining the control area of each rainfall station.
In this embodiment, the rainfall characteristic index is to determine the total rainfall at a certain rainfall intensity, and determine the total rainfall of the whole drainage basin by determining the single rainfall of the single rainfall station.
The beneficial effects of the above technical scheme are: the rainfall amount excerpted data with short duration is fully utilized, the rainfall amount and the rainfall duration of each rainfall record are counted one by one, the observation time precision is utilized to the maximum extent, and the comprehensive rainfall characteristic index capable of reflecting the rainfall intensity and the rainfall amount change condition of the drainage basin is constructed. The method has the advantages of simple steps, clear concept, simple and convenient calculation and easy operation, is an effective tool for researching rainfall change characteristics of the drainage basin, and particularly has positive significance for researching rainfall change in the area with the super-seepage as a main area.
The invention provides a rainfall characteristic index statistical method based on short-duration data, and a step 4 of superposing the total rainfall amounts reaching specific rainfall intensity corresponding to all rainfall stations in a river basin to obtain a second total rainfall amount reaching specific rainfall intensity corresponding to the river basin comprises the following steps:
calculating a second total rainfall amount corresponding to the drainage basin reaching a specific rainfall intensity according to the following formula:
Figure BDA0003313983530000091
in the formula: wrA second total rainfall amount corresponding to the drainage basin and reaching a specific rainfall intensity (r mm per minute) or more, wherein the unit is billionth of cubic meters; j is the jth rainfall station with the intensity more than r millimeters per minute in the drainage basin; m is the total number of the rainfall stations reaching the strength of more than r millimeters per minute in the drainage basin; n is the total rainfall record number of the rainfall station corresponding to the intensity of more than r millimeters per minute;
Figure BDA0003313983530000092
recording corresponding rainfall capacity in millimeter for the ith rainfall with the rainfall intensity of more than r millimeter per minute in the preset month flow domain; a is the control area of the corresponding rainfall station, and the unit is square kilometer,
Figure BDA0003313983530000093
and the j-th rainfall station is expressed to reach a first rainfall total amount corresponding to the intensity of more than r millimeters per minute in the preset month.
The beneficial effects of the above technical scheme are: by counting the rainfall amount and the rainfall duration of each rainfall record one by one, the observation time precision is utilized to the maximum extent, and a comprehensive rainfall characteristic index capable of reflecting the rainfall intensity and the rainfall change condition of the drainage basin is constructed.
The invention provides a rainfall characteristic index statistical method based on short-duration data, which comprises the following steps of 1, before collecting and calculating rainfall extract data of each rainfall station in a flow domain, the method further comprises the following steps:
determining the current position of each rainfall station, constructing a historical weather model in a control area corresponding to the control area of each rainfall station according to the weather change condition of a preset month in a historical year, and simulating the rainfall distribution of the control area according to the historical weather model;
constructing regional terrain of a control region corresponding to each rainfall station, and carrying out terrain type division on the regional terrain;
carrying out region matching on the rainfall distribution and the divided terrain results one by one, and determining an independent rainfall measurement region and a dependent rainfall measurement region in each control region;
acquiring a distribution mode according to the regional attribute of the dependent rainfall measurement region, and distributing corresponding auxiliary rainfall stations according to the distribution mode, wherein the auxiliary rainfall stations are used for monitoring and recording the rainfall of the corresponding dependent rainfall measurement region;
the auxiliary area is in the control area, and the number of the auxiliary rainfall stations is at least one.
As shown in fig. 2, for example, the area a represents a control area corresponding to the rainfall station 1, and the historical weather model is obtained (for example, the historical weather model corresponding to the control area is obtained from the weather data between 2000 and 2021), so as to simulate the rainfall distribution of the control area in a preset month, and the rainfall distribution can be estimated according to the estimated rainfall in the history of 7 to 8 months, so as to obtain a distribution composed of different rainfall, and since the different areas have different terrains and may have mountains, plains and the like, an independent rainfall measurement area (i.e., an effective measurement area corresponding to the rainfall station) is obtained by performing area-to-area matching between the rainfall distribution and the corresponding area, for example, the sub-areas a1, a2, a3 and a4 in the area a are areas of normal rainfall, and the remaining areas a5 and a6 are areas of abnormal rainfall, c is an effective measuring area, wherein c is an independent rainfall measuring area, and a1-a6 is a non-independent rainfall measuring area.
The beneficial effects of the above technical scheme are: the control area corresponding to the original control area of the rainfall station is divided according to the weather model to obtain different measurement areas, an effective accurate basis can be provided for follow-up rainfall measurement, the accuracy of follow-up rainfall measurement is guaranteed, a data basis is provided for short-time measurement, and the accurate acquisition of the total rainfall of the drainage basin is improved.
The invention provides a rainfall characteristic index statistical method based on short-duration data, and the method comprises the following steps of 4, after obtaining a second rainfall total amount corresponding to a drainage basin and reaching specific rainfall intensity: and performing correction adjustment on the second total rainfall, which comprises the following steps:
marking a first rainfall station which is initially set;
acquiring a control area corresponding to the first rainfall station and taking the control area as an initial area Aj；
Acquiring monitoring information of each first rainfall station reaching a specific rainfall intensity in a preset month, wherein the monitoring information comprises: a first time point when the specific rainfall is reached, a second time set of continuous rainfall before the specific rainfall is reached, and a third time set of continuous rainfall after the specific rainfall is reached;
acquiring intra-domain weather change information of a sky area of the size corresponding to the initial area of each first rainfall station from a weather observation database corresponding to the preset month, and constructing a weather change model corresponding to a time period according to the intra-domain weather change information, wherein the intra-domain weather change information is related to cloud coverage area, cloud thickness and cloud moving speed;
according to the weather change model, the weather change layers corresponding to the specific rainfall intensity related to the first rainfall station at different time points are divided and output, and are placed in an overlapping mode to obtain a stability maintaining area and a residual area of the first rainfall station, wherein the areas of the stability maintaining area and the residual area form corresponding initial areas;
meanwhile, according to overlapped color results after overlapping, performing region division on the remaining region to obtain a plurality of sub-regions, calling and monitoring the rainfall of the corresponding sub-regions based on the fact that a second rainfall station is arranged in the sub-regions in advance, and obtaining the optimized rainfall corresponding to the initial area by combining the rainfall of the stable region monitored by the corresponding first rainfall station;
based on the optimized rainfall, the second total rainfall is corrected and adjusted to obtain a third total rainfall, the third total rainfall is compared with the second total rainfall, and if the absolute value of the difference value of the third total rainfall and the second total rainfall is within a preset difference value range, rainfall data monitored by the first rainfall station is continuously used for estimating the total rainfall;
otherwise, estimating the total rainfall of the rainfall data monitored by the first rainfall station and the second rainfall station;
preferably, the dimensionally stable region includes an independent rainfall measurement region, and the dependent rainfall measurement region includes a remaining region;
the second rainfall stations are obtained by screening in auxiliary rainfall stations arranged in the corresponding controlled area, and the number of the auxiliary rainfall stations is larger than that of the corresponding second rainfall stations.
In this embodiment, the first rain station is the initial default rain station in each control area.
In this embodiment, the monitoring information of the specific raininess is to obtain a time period under the characteristic index.
In this embodiment, as shown in fig. 3, 1, 2, and 3 may be regarded as different weather change layers, and the cloud states in the weather regions corresponding to the control regions are different, so that the cloud state block diagrams set in different segments in the weather regions are represented by thick and thin filling lines, and since the initial setting of the first rainfall station is determined according to a certain existing reference technique, the stable region and the remaining region of the first rainfall station may be obtained by combining the overlapping placement results, the corresponding d1 region in the diagram may be regarded as the stable region, and d2 is the remaining region.
In this embodiment, the determination of the rainfall stations in the remaining area d1 is based on and selected from rainfall stations preset in non-independent rainfall measurement areas.
The beneficial effects of the above technical scheme are: the weather change model is constructed by acquiring the detection information and the weather information of specific rainfall intensity, the overlapping mode of layers is carried out after the separation, the stable region and the residual region are acquired, the rainfall corresponding to the initial area is optimized by calling the preset second rainfall station, and then the optimization of the total amount of the second rainfall is realized, and the accuracy of the rainfall is ensured to be acquired under the characteristic index.
The invention provides a rainfall characteristic index statistical method based on short-duration data, which is used for acquiring optimized rainfall corresponding to an initial area and comprises the following steps:
Figure BDA0003313983530000121
wherein, YjRepresenting the optimized rainfall of the original corresponding jth initial area; g represents the total number of the sub-areas of the area corresponding to the jth initial area, a second rainfall station is started in each sub-area, and the total number of the second rainfall stations is also G; fg1,iIndicating that the recorded rainfall corresponding to the g1 th second rainfall station when the recorded area corresponding to the j initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute; sg1,iRepresenting the area of the corresponding g1 th sub-area when the recorded area corresponding to the jth initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute; fj,iRecording rainfall measured by a corresponding first rainfall station when the area corresponding to the jth initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute; sj,iRepresenting the area of the corresponding stable region when the recorded corresponding region of the jth initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute; deltag1,iThe effective measurement coefficient of the corresponding rainfall station is represented when the recorded corresponding area of the jth initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute, and the value range is (0.8, 1); deltaj,iThe effective measurement coefficient of the corresponding rainfall station is represented when the recorded corresponding area of the jth initial area reaches the ith rainfall with the rainfall intensity of more than r millimeters per minute, and the value range is (0.8, 1);
wherein when is deltag1,iWhen the content is equal to 0, the content,
Figure BDA0003313983530000122
when deltaj,iWhen the content is equal to 0, the content,
Figure BDA0003313983530000131
in one possible implementation manner, performing a correction adjustment on the second total rainfall amount based on the optimized rainfall amount to obtain a third total rainfall amount includes:
Figure BDA0003313983530000132
wherein, WsIs the third total amount of rainfall present,
Figure BDA0003313983530000133
the ratio of the area corresponding to the jth initial area to the second rainfall station when the ith area reaches the rainfall with the rainfall intensity of more than r millimeters per minute; when in use
Figure BDA0003313983530000134
When the ratio of (a) is greater than or equal to 1/2, a takes the value of 1; when in use
Figure BDA0003313983530000135
When the ratio of (a) is less than 1/2, a is 2.
The beneficial effects of the above technical scheme are: the accurate rainfall of the sub-areas set in the corresponding control area of each first rainfall station is obtained, the optimized rainfall of all the sub-areas in the corresponding control area is determined according to the accurate rainfall, the rainfall of the control area under different weather conditions at the same time point is conveniently and accurately obtained, the total rainfall of the control area is accurately obtained, the second total rainfall is corrected and adjusted finally, the third total rainfall is effectively obtained, the subsequent effective comparison of the rainfall is convenient, reasonable labor division is carried out on the area size of the control area, and the measuring accuracy is further guaranteed.
The invention provides a rainfall characteristic index statistical method based on short-duration data, which comprises the following steps before calculating the rainfall intensity of each rainfall record:
determining the rainfall time period of each rainfall station according to the starting time and the ending time of each rainfall excerpted data of each rainfall station in a preset month, and simultaneously acquiring weather data corresponding to each rainfall time period;
calculating a first matching degree P1 between the rainfall time length of the rainfall time period corresponding to the rainfall station and the estimated rainfall time length corresponding to the weather data according to the following formula;
Figure BDA0003313983530000136
wherein F1 represents the total rainfall time period in the preset month measured by the same rainfall station; p1f1A rainfall period representing the f1 th rainfall period measured at the same rainfall station; p2f1Indicating an estimated rainfall duration determined based on weather data when measuring the rainfall duration of the f1 th rainfall time period of the same rainfall station;
calculating a second matching degree P2 between the rainfall capacity of the rainfall time period corresponding to the rainfall station and the estimated rainfall capacity corresponding to the weather data according to the following formula;
Figure BDA0003313983530000141
wherein F1 represents the total rainfall time period in the preset month measured by the same rainfall station; y is1f1Represents the rainfall of the f1 th rainfall period measured by the same rainfall station; y is2f1Represents an estimated rainfall determined based on weather data when measuring the rainfall at the same rainfall station at the f1 th rainfall period; df1The weather correction parameters based on the weather data are represented, and the value range is (0.1, 0.2);
when the first matching degree P1 is greater than the corresponding first preset matching degree and the second matching degree P2 is greater than the corresponding second preset matching degree, judging that the corresponding rainfall station is qualified, and controlling the qualified rainfall station to continue rainfall measurement;
otherwise, collecting the operation information of the unqualified rainfall station in a preset month, and analyzing the unqualified type of the unqualified rainfall station;
when the unqualified type is irrelevant to equipment upgrading, alarming and reminding are carried out, meanwhile, a standby rainfall station is started to carry out measurement work, meanwhile, a corresponding operation information is reversely analyzed based on a historical estimation model of the unqualified rainfall station, and the estimated rainfall of the unqualified rainfall station is estimated to replace a rainfall record in unqualified operation;
and when the unqualified type is related to equipment upgrading, upgrading the corresponding unqualified type of the corresponding unqualified rainfall station.
In this embodiment, the estimated rainfall duration is obtained according to the corresponding weather data, and the required duration is obtained by inputting the weather data into the pre-trained duration obtaining model.
In this embodiment, the unqualified type refers to a type of the rainfall station, where the operation parameter is related to software without or without messy codes, or a type of hardware failure and structural damage of the rainfall station itself.
In this embodiment, the reverse analysis means, for example, that the historical estimation model is a normal operation model of the estimated rainfall station under the normal operation condition, and the normal operation model includes various normal parameters and normal rainfall corresponding to the various normal parameters, so according to the parameters corresponding to the unqualified type and the unqualified type, the corresponding inaccurate recorded rainfall under the unqualified condition can be reversely determined, and further, the inaccurate recorded rainfall can be used as replacement information to replace the unqualified information, so as to perform recording.
The beneficial effects of the above technical scheme are: whether the corresponding rainfall station has a problem or not is determined through the two aspects of rainfall duration and rainfall, wherein the rainfall station is determined comprehensively through comparison of the calculated matching degree and the preset matching degree, the accuracy of the rainfall recorded in the initial process is guaranteed, normal operation of the unqualified rainfall station is guaranteed through adjustment processing of relevant unqualified types of the unqualified rainfall station, the validity of measured data is further guaranteed, and an effective data basis is provided for the total rainfall of the drainage basin obtained after the specific rainfall is achieved subsequently.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.

Claims (8)

1.一种基于短历时数据的降雨特征指标统计方法，其特征在于，包括：1. a rainfall characteristic index statistical method based on short duration data, is characterized in that, comprises: 步骤1，收集计算流域内各雨量站降水量摘录数据，同时，确定各雨量站的控制面积；Step 1: Collect and calculate the extracted data of the precipitation of each rainfall station in the basin, and at the same time, determine the control area of each rainfall station; 步骤2，根据各雨量站在预设月份的每条降水量摘录数据的开始、结束时间及降雨量，计算每条降雨记录的降雨强度；Step 2: Calculate the rainfall intensity of each rainfall record according to the start, end time and rainfall of each precipitation extract data of each rainfall station in the preset month; 步骤3，按照不同的降雨强度标准对单个雨量站降雨记录进行筛选，将达到特定雨强的所有降雨记录的降雨量累加，并与该雨量站控制面积进行乘积，得到单个雨量站对应的达到特定雨强的第一降雨总量；Step 3: Screen the rainfall records of a single rainfall station according to different rainfall intensity standards, accumulate the rainfall of all rainfall records that reach a specific rainfall intensity, and multiply with the control area of the rainfall station to obtain the corresponding rainfall of a single rainfall station. The first total rainfall of the rain intensity; 步骤4，将流域内各雨量站对应的达到特定雨强的降雨总量叠加，得到对应流域达到特定雨强的第二降雨总量。Step 4, superimpose the total amount of rainfall reaching a specific rain intensity corresponding to each rainfall station in the watershed, to obtain a second total rainfall amount corresponding to the watershed reaching the specific rain intensity. 2.如权利要求1所述的基于短历时数据的降雨特征指标统计方法，其特征在于，步骤4，将流域内各雨量站对应的达到特定雨强的降雨总量叠加，得到对应流域达到特定雨强的第二降雨总量，包括：2. the rainfall characteristic index statistical method based on short duration data as claimed in claim 1, is characterized in that, in step 4, the total amount of rainfall that reaches the specific rain intensity corresponding to each rainfall station in the watershed is superimposed, and obtains that the corresponding watershed reaches the specific rainfall. The second rainfall total of the rain intensity, including: 根据如下公式，计算对应流域达到特定雨强的第二降雨总量：According to the following formula, calculate the second total rainfall that reaches a certain rainfall intensity in the corresponding watershed:
Figure FDA0003313983520000011
Figure FDA0003313983520000011
 式中：Wr为所述流域对应的达到特定雨强(r毫米每分钟)以上的第二降雨总量，单位为亿立方米；j为流域内第j个大于r毫米每分钟以上强度的雨量站；m为流域内达到r毫米每分钟以上强度雨量站的总数；n为达到r毫米每分钟以上强度对应的雨量站的降雨记录总数；
Figure FDA0003313983520000012
为预设月份流域内第i个r毫米每分钟以上雨强的降雨记录相应的降雨量，单位为毫米；A为相应雨量站的控制面积，单位为平方公里，
Figure FDA0003313983520000013
表示第j个雨量站在预设月份达到r毫米每分钟以上强度对应的第一降雨总量。In the formula: W r is the second total amount of rainfall that reaches a specific rainfall intensity (r mm per minute) or more corresponding to the watershed, the unit is 100 million cubic meters; j is the jth in the watershed greater than the intensity of r mm per minute or more. Rain gauge station; m is the total number of rain gauge stations with an intensity of more than r mm per minute in the basin; n is the total number of rainfall records corresponding to the rain gauge station with an intensity of more than r mm per minute;
Figure FDA0003313983520000012
is the corresponding rainfall record for the i-th r mm of rain intensity per minute or more in the basin in the preset month, the unit is mm; A is the control area of the corresponding rainfall station, the unit is square kilometers,
Figure FDA0003313983520000013
Indicates the first total rainfall corresponding to the jth rainfall station reaching an intensity of more than r mm per minute in the preset month. 3.如权利要求1所述的基于短历时数据的降雨特征指标统计方法，其特征在于，步骤1，收集计算流域内各雨量站降水量摘录数据之前，还包括：3. the rainfall characteristic index statistical method based on short duration data as claimed in claim 1, is characterized in that, step 1, before collecting and calculating the precipitation extract data of each rainfall station in the basin, also comprises: 确定每个雨量站的当前位置，并根据历史年份中对应预设月份的天气变化情况，构建每个雨量站所在控制面积对应的控制区域内的历史天气模型，并根据所述历史天气模型，模拟所述控制区域的降雨分布；Determine the current position of each rainfall station, and build a historical weather model in the control area corresponding to the control area where each rainfall station is located according to the weather changes in the corresponding preset months in the historical year, and simulate the historical weather model according to the historical weather model. Rainfall distribution in said control area; 构建每个雨量站对应控制区域的区域地形，并对所述区域地形进行地形类型划分；Construct the regional terrain of the control area corresponding to each rainfall station, and divide the terrain type of the regional terrain; 将所述降雨分布与划分的地形结果进行一一的区域匹配，确定每个控制区域中独立雨量测量区域以及非独立雨量测量区域；Perform regional matching between the rainfall distribution and the divided terrain results, and determine the independent rainfall measurement area and the dependent rainfall measurement area in each control area; 按照所述非独立雨量测量区域的区域属性，获取分配方式，并按照所述分配方式，进行相应辅助雨量站的分配，且所述辅助雨量站，用于监测对应的所述非独立雨量测量区域的降雨量，并记录；According to the regional attribute of the dependent rainfall measurement area, the distribution method is obtained, and according to the distribution method, the corresponding auxiliary rainfall stations are allocated, and the auxiliary rainfall stations are used to monitor the corresponding dependent rainfall measurement area. rainfall, and record; 其中，所述辅助区域在控制区域内，且所述辅助雨量站的个数至少为一。Wherein, the auxiliary area is within the control area, and the number of the auxiliary rainfall stations is at least one. 4.如权利要求3所述的基于短历时数据的降雨特征指标统计方法，其特征在于，步骤4，在得到对应流域达到特定雨强的第二降雨总量之后，还包括：对所述第二降雨总量进行修正调整，其包括：4. The rainfall characteristic index statistical method based on short-duration data as claimed in claim 3, characterized in that, in step 4, after obtaining the second total amount of rainfall that the corresponding watershed reaches a specific rain intensity, it also comprises: 2. The total rainfall is revised and adjusted, which includes: 标记初始设定的第一雨量站；Mark the initial set of first rainfall station; 获取所述第一雨量站对应的控制面积，并作为初始面积Aj；Obtain the control area corresponding to the first rainfall station, and use it as the initial area A j ; 获取每个第一雨量站在预设月份内，获取达到特定雨强的监测信息，所述监测信息包括：达到特定雨强的第一时间点、达到特定雨强前的连续降雨的第二时间集合、达到特定雨强后的连续降雨的第三时间集合；Acquiring each first rainfall station within a preset month, and obtaining monitoring information that reaches a specific rain intensity, the monitoring information includes: the first time point when the specific rain intensity is reached, and the second time of continuous rainfall before reaching the specific rain intensity Collection, the third time collection of continuous rainfall after reaching a specific rain intensity; 从所述预设月份对应的天气观测数据库中，获取每个第一雨量站的初始面积所对应面积大小的天空区域的域内天气变化信息，并按照所述域内天气变化信息构建对应时间段的天气变化模型，所述域内天气变化信息与云覆盖面积、云厚度、云移动速度相关；From the weather observation database corresponding to the preset month, obtain the intra-domain weather change information of the sky area corresponding to the initial area of each first rainfall station, and construct the weather for the corresponding time period according to the intra-domain weather change information a change model, where the weather change information in the domain is related to the cloud coverage area, cloud thickness, and cloud moving speed; 根据所述天气变化模型，拆分并输出构建与所述第一雨量站相关的特定雨强在不同时间点下对应的天气变化层，并将所述天气变化层进行重叠放置，获得第一雨量站的维稳区域以及剩余区域，其中，所述维稳区域以及剩余区域的面积构成对应的初始面积；According to the weather change model, split and output and construct the weather change layers corresponding to the specific rain intensity related to the first rainfall station at different time points, and overlap the weather change layers to obtain the first rainfall The stability maintenance area and the remaining area of the station, wherein the areas of the stability maintenance area and the remaining area constitute the corresponding initial area; 同时，按照重叠放置后的重叠颜色结果，对所述剩余区域进行区域划分，获得若干子区域，并基于预先在所述子区域内设置第二雨量站，来调取监测得到的对应子区域的降雨量，并结合对应第一雨量站所监测的所述维稳区域的降雨量，来获取初始面积对应的优化降雨量；At the same time, according to the overlapping color results after overlapping placement, the remaining area is divided into regions to obtain several sub-regions, and based on pre-setting a second rainfall station in the sub-regions, the corresponding sub-regions obtained by monitoring are retrieved. rainfall, and combined with the rainfall in the stability maintenance area monitored by the first rainfall station to obtain the optimized rainfall corresponding to the initial area; 基于所述优化降雨量，来对所述第二降雨总量进行修正调整，得到第三降雨总量，并与第二降雨总量进行比较，若两者的差值绝对值在预设差值范围内，继续使用第一雨量站监测的降雨数据进行总降雨量的估算；Based on the optimized rainfall, the second total rainfall is corrected and adjusted to obtain a third total rainfall, and compared with the second total rainfall, if the absolute value of the difference between the two is within the preset difference Within the range, continue to use the rainfall data monitored by the first rainfall station to estimate the total rainfall; 否则，将所述第一雨量站与第二雨量站监测的降雨数据进行总降雨量的估算。Otherwise, the total rainfall is estimated based on the rainfall data monitored by the first rainfall station and the second rainfall station. 5.如权利要求4所述的基于短历时数据的降雨特征指标统计方法，其特征在于，获取初始面积对应的优化降雨量，包括：5. the rainfall characteristic index statistical method based on short duration data as claimed in claim 4, is characterized in that, obtains the optimized rainfall corresponding to initial area, comprises:
Figure FDA0003313983520000031
Figure FDA0003313983520000031
 其中，Yj表示原始对应的第j个初始面积的优化降雨量；G表示所述第j个初始面积对应区域的子区域的总个数，且每个子区域中开启有一个第二雨量站，且第二雨量站的总个数也为G；Fg1,i表示记录的所述第j个初始面积对应区域达到第i个r毫米每分钟以上雨强的降雨时，对应第g1个第二雨量站测量的降雨量；Sg1，i表示记录的所述第j个初始面积对应区域达到第i个r毫米每分钟以上雨强的降雨时，对应的第g1个子区域的区域面积；Fj，i记录的所述第j个初始面积对应区域达到第i个r毫米每分钟以上雨强的降雨时，对应的第一雨量站测量的降雨量；Sj，i表示记录的所述第j个初始面积对应区域达到第i个r毫米每分钟以上雨强的降雨时，对应的维稳区域的区域面积；δg1，i表示记录的所述第j个初始面积对应区域达到第i个r毫米每分钟以上雨强的降雨时，对应雨量站的有效测量系数，且取值范围为(0.8,1)；δj,i表示记录的所述第j个初始面积对应区域达到第i个r毫米每分钟以上雨强的降雨时，对应雨量站的有效测量系数，且取值范围为(0.8,1)；Among them, Y j represents the optimized rainfall of the original corresponding j-th initial area; G represents the total number of sub-regions of the area corresponding to the j-th initial area, and each sub-region has a second rainfall station enabled, And the total number of second rainfall stations is also G; F g1,i indicates that when the recorded area corresponding to the jth initial area reaches the ith r mm/min or more rainfall intensity, the corresponding g1 second The rainfall measured by the rainfall station; S g1, i represents the area of the corresponding g1 sub-area when the recorded area corresponding to the j-th initial area reaches the i-th r mm/min rainfall intensity; F j , when the corresponding area of the jth initial area recorded by i reaches the ith r mm/min rainfall intensity, the corresponding rainfall measured by the first rainfall station; S j, i represents the recorded jth When the area corresponding to the initial area reaches the i-th r mm or more rainfall per minute, the area of the corresponding stability maintenance area; δ g1, i indicates that the recorded area corresponding to the j-th initial area reaches the i-th r mm When it rains with a rainfall intensity of more than one minute, the effective measurement coefficient of the corresponding rainfall station, and the value range is (0.8 , 1); When it rains with a rainfall intensity of more than one minute, the effective measurement coefficient of the corresponding rainfall station, and the value range is (0.8, 1); 其中，当δg1,i＝0时，
Figure FDA0003313983520000041
Among them, when δ g1,i =0,
Figure FDA0003313983520000041
 当δj,i＝0时，
Figure FDA0003313983520000042
When δ j,i =0,
Figure FDA0003313983520000042
 6.如权利要求5所述的基于短历时数据的降雨特征指标统计方法，其特征在于，基于所述优化降雨量，来对所述第二降雨总量进行修正调整，得到第三降雨总量，包括：6 . The statistical method for rainfall characteristic indicators based on short-duration data according to claim 5 , wherein, based on the optimized rainfall, the second total rainfall is corrected and adjusted to obtain a third total rainfall. 7 . ,include:
Figure FDA0003313983520000043
Figure FDA0003313983520000043
 其中，Ws为所述第三降雨总量，
Figure FDA0003313983520000044
表示第j个初始面积对应区域达到第i个r毫米每分钟以上雨强的降雨时，与第二雨量站相关的比值；当
Figure FDA0003313983520000045
的比值大于或等于1/2时，a取值为1；当
Figure FDA0003313983520000046
的比值小于1/2时，a取值为2。Wherein, W s is the third total rainfall amount,
Figure FDA0003313983520000044
Represents the ratio related to the second rainfall station when the area corresponding to the j-th initial area reaches the i-th r mm/min rainfall intensity or more; when
Figure FDA0003313983520000045
When the ratio of is greater than or equal to 1/2, the value of a is 1; when
Figure FDA0003313983520000046
When the ratio is less than 1/2, a takes the value 2. 7.如权利要求4所述的基于短历时数据的降雨特征指标统计方法，其特征在于，7. The rainfall characteristic index statistical method based on short duration data as claimed in claim 4, is characterized in that, 所述维稳区域是包括独立雨量测量区域在内的，且所述非独立雨量测量区域是包括剩余区域在内的；The stability maintenance area includes the independent rainfall measurement area, and the non-independent rainfall measurement area includes the remaining area; 所述第二雨量站是对应受控区域中设置的辅助雨量站中筛选得到的，且所述辅助雨量站的数量大于对应的第二雨量站的数量。The second rainfall stations are selected from auxiliary rainfall stations set in the corresponding controlled area, and the number of the auxiliary rainfall stations is greater than the number of the corresponding second rainfall stations. 8.如权利要求1所述的基于短历时数据的降雨特征指标统计方法，其特征在于，计算每条降雨记录的降雨强度之前，还包括：8. the rainfall characteristic index statistical method based on short duration data as claimed in claim 1 is characterized in that, before calculating the rainfall intensity of each rainfall record, also comprises: 根据各雨量站在预设月份的每条降水量摘录数据的开始、结束时间，确定每个雨量站的降雨时间段，同时，获取每个降雨时间段对应的天气数据；Determine the rainfall time period of each rainfall station according to the start and end time of each precipitation extract data of each rainfall station in the preset month, and at the same time, obtain the weather data corresponding to each rainfall time period; 根据如下公式，计算所述雨量站对应的所述降雨时间段的降雨时长与所述天气数据对应的预估降雨时长之间的第一匹配度P1；Calculate the first matching degree P1 between the rainfall duration of the rainfall time period corresponding to the rainfall station and the estimated rainfall duration corresponding to the weather data according to the following formula;
Figure FDA0003313983520000047
Figure FDA0003313983520000047
 其中，F1表示同个雨量站测量的预设月份内的总降雨时间段；P1f1表示同个雨量站测量的第f1个降雨时间段的降雨时长；P2f1表示在测量同个雨量站第f1个降雨时间段的降雨时长时，基于天气数据确定的预估降雨时长；Among them, F1 represents the total rainfall period in the preset month measured by the same rainfall station; P 1f1 represents the rainfall duration of the f1th rainfall period measured by the same rainfall station; P 2f1 represents the f1th rainfall period measured by the same rainfall station. The estimated rainfall duration determined based on the weather data when the rainfall duration of each rainfall time period; 根据如下公式，计算所述雨量站对应的所述降雨时间段的降雨量与基于所述天气数据对应的预估降雨量之间的第二匹配度P2；According to the following formula, calculate the second matching degree P2 between the rainfall amount of the rainfall time period corresponding to the rainfall station and the estimated rainfall amount corresponding to the weather data;
Figure FDA0003313983520000051
Figure FDA0003313983520000051
 其中，F1表示同个雨量站测量的预设月份内的总降雨时间段；Y1f1表示同个雨量站测量的第f1个降雨时间段的降雨量；Y2f1表示在测量同个雨量站第f1个降雨时间段的降雨量时，基于天气数据确定的预估降雨量；Df1表示基于天气数据的天气校正参数，且取值范围为(0.1,0.2)；Among them, F1 represents the total rainfall time period in the preset month measured by the same rainfall station; Y 1f1 represents the rainfall in the f1th rainfall time period measured by the same rainfall station; Y 2f1 represents the f1th rainfall measured by the same rainfall station. The estimated rainfall determined based on the weather data when the rainfall amount in each rainfall time period; D f1 represents the weather correction parameter based on the weather data, and the value range is (0.1, 0.2); 当所述第一匹配度P1大于对应的第一预设匹配度时，且所述第二匹配度P2大于对应的第二预设匹配度时，判定对应的雨量站合格，并控制合格雨量站继续进行雨量测量；When the first matching degree P1 is greater than the corresponding first preset matching degree, and the second matching degree P2 is greater than the corresponding second preset matching degree, it is determined that the corresponding rainfall station is qualified, and the qualified rainfall station is controlled Continue rain measurement; 否则，采集不合格雨量站在预设月份的运行信息,分析所述不合格雨量站的不合格类型；Otherwise, collect the operation information of the unqualified rainfall station in the preset month, and analyze the unqualified type of the unqualified rainfall station; 当所述不合格类型与设备升级无关时，进行报警提醒，同时，启动备用雨量站进行测量工作，同时，基于所述不合格雨量站的历史估计模型对对应的运行信息进行反向分析，预估所述不合格雨量站的估计雨量，来替换处于不合格运行时的降雨记录；When the unqualified type has nothing to do with the equipment upgrade, an alarm reminder is given, and at the same time, the standby rainfall station is started to perform measurement work. estimating the estimated rainfall of said unqualified rain gauge station to replace the rainfall record at the time of unqualified operation; 当所述不合格类型与设备升级有关时，对对应不合格雨量站进行对应不合格类型的升级。When the unqualified type is related to equipment upgrade, the corresponding unqualified rainfall station is upgraded corresponding to the unqualified type.
CN202111225876.2A 2021-10-21 2021-10-21 Rainfall characteristic index statistical method based on short-duration data Pending CN113961623A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111225876.2A CN113961623A (en) 2021-10-21 2021-10-21 Rainfall characteristic index statistical method based on short-duration data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111225876.2A CN113961623A (en) 2021-10-21 2021-10-21 Rainfall characteristic index statistical method based on short-duration data

Publications (1)

Publication Number Publication Date
CN113961623A true CN113961623A (en) 2022-01-21

Family

ID=79465801

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111225876.2A Pending CN113961623A (en) 2021-10-21 2021-10-21 Rainfall characteristic index statistical method based on short-duration data

Country Status (1)

Country Link
CN (1) CN113961623A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114781769A (en) * 2022-06-27 2022-07-22 长江水利委员会水文局 A dynamic and refined early warning method for over-standard floods in a watershed

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170261646A1 (en) * 2016-03-08 2017-09-14 China Institute Of Water Resources And Hydropower Research Self-correcting multi-model numerical rainfall ensemble forecasting method
CN108951531A (en) * 2018-06-22 2018-12-07 中国水利水电科学研究院 The construction method of watershed model time point investigation based on landform scale bar
CN109376940A (en) * 2018-11-02 2019-02-22 中国水利水电科学研究院 Method and device for obtaining temporal and spatial distribution law of rainfall during rainfall
JP2020092367A (en) * 2018-12-07 2020-06-11 日本放送協会 Rain attenuation prediction device and program
CN112785053A (en) * 2021-01-15 2021-05-11 北京市水科学技术研究院 Method and system for forecasting urban drainage basin flood
CN113449435A (en) * 2021-07-22 2021-09-28 北京大学 Multi-objective optimization method for sponge urban engineering design

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170261646A1 (en) * 2016-03-08 2017-09-14 China Institute Of Water Resources And Hydropower Research Self-correcting multi-model numerical rainfall ensemble forecasting method
CN108951531A (en) * 2018-06-22 2018-12-07 中国水利水电科学研究院 The construction method of watershed model time point investigation based on landform scale bar
CN109376940A (en) * 2018-11-02 2019-02-22 中国水利水电科学研究院 Method and device for obtaining temporal and spatial distribution law of rainfall during rainfall
JP2020092367A (en) * 2018-12-07 2020-06-11 日本放送協会 Rain attenuation prediction device and program
CN112785053A (en) * 2021-01-15 2021-05-11 北京市水科学技术研究院 Method and system for forecasting urban drainage basin flood
CN113449435A (en) * 2021-07-22 2021-09-28 北京大学 Multi-objective optimization method for sponge urban engineering design

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JEAN-LUC MARTEL等: "Climate Change and Rainfall Intensity–Duration–Frequency Curves: Overview of Science and Guidelines for Adaptation", JOURNAL OF HYDROLOGIC ENGINEERING, vol. 26, no. 10, 3 August 2021 (2021-08-03) *
周祖昊;贾仰文;王浩;朱厚华;王建华;严登华;: "大尺度流域基于站点的降雨时空展布", 水文, no. 01, 25 February 2006 (2006-02-25) *
幸新涪;周火明;秦维;谭昭娣;卢宇;张乾柱;卢阳;: "基于流量-面积比值的小流域设计洪水计算方法对比研究", 长江科学院院报, no. 01, 15 January 2018 (2018-01-15) *
陈松伟;王鹏;李超群;: "基于短历时降雨数据的降雨输沙经验模型研究", 人民黄河, no. 06, 10 June 2019 (2019-06-10) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114781769A (en) * 2022-06-27 2022-07-22 长江水利委员会水文局 A dynamic and refined early warning method for over-standard floods in a watershed

Similar Documents

Publication Publication Date Title
CN113222296B (en) Flood control scheduling method based on digital twinning
EP1801682B1 (en) Information analysis system for water distribution and pipelines
US8594936B1 (en) Providing current estimates of precipitation accumulations
CN112800379B (en) MODIS remote sensing snow information processing method and device
WO2020253338A1 (en) Traffic big data-based road capacity extraction method
CN119066610A (en) A comprehensive monitoring system of hydrology and water quality based on multi-source data fusion
Sauer Standards for the analysis and processing of surface-water data and information using electronic methods
CN118211721B (en) Building energy consumption prediction method based on artificial intelligence
CN114565279A (en) Highway sanitation vehicle operation assessment method based on mathematical modeling
CN111429028A (en) Power transmission line icing disaster risk assessment method suitable for mountainous terrain
CN106650239A (en) Climate zone-based station rainfall observation error correction method and system
CN113269240B (en) Rainfall station site selection information output method and device, electronic equipment and medium
CN116804569A (en) Railway subgrade collapse monitoring method and system
CN113961623A (en) Rainfall characteristic index statistical method based on short-duration data
CN108877225B (en) Traffic flow index determination method and device
CN118294002A (en) Ecological flow calibration method and system based on three-dimensional hydrodynamic model
CN117152926A (en) Automatic pressure detection early warning system based on data analysis
CN117033935B (en) Prediction method of rainfall characteristic under statistics and monitoring based on Bayesian fusion
Osman et al. Estimating the Clark instantaneous unit hydrograph parameters for selected gauged catchments in the west Coast of Peninsular Malaysia
CN113888043B (en) A full-cycle visual management and analysis system for girder defects of beam bridges
CN110020000A (en) Determine the method and apparatus of abnormal air speed data
CN118551656B (en) River flow velocity prediction method, device and medium
CN118503658B (en) Mountain front water system extraction method based on optical remote sensing data and digital elevation model
CN119509466A (en) Dynamic environment settlement observation method
CN119646431A (en) Ice thickness prediction method and system based on micro-meteorological data

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