Hydrological element non-stationary sequence division method, system, device and storage medium

Abstract

The embodiment of the application provides a hydrological element non-stationary sequence division method, system, device and storage medium, and belongs to the technical field of data processing. The method first acquires a first sequence about annual hydrological element indexes, then constructs a plurality of time-continuous first sample sequences according to the first sequence, and performs data distribution consistency test between two first sample sequences to determine a plurality of first sub-sequences of the hydrological element stationary distribution in the first sequence, acquires the daily hydrological element index value of each date under the interannual cycle corresponding to the first sub-sequence to construct a second sequence, constructs a plurality of time-continuous second sample sequences according to the second sequence, and performs consistency test between the second sample sequences and a target distribution function to determine a plurality of second sub-sequences of the hydrological element stationary distribution in the second sequence. The embodiment of the application realizes annual and intra-annual double-scale division, and can improve the reliability of hydrological element non-stationary sequence division.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of data processing, and particularly relates to a hydrological element non-stationary sequence division method, system, device and storage medium. BACKGROUND

[0002] Affected by the changing environment, the hydrological regime of the basin, the underlying surface conditions and the runoff generation and collection mechanism change, which superimposes the non-stationary and nonlinear characteristics of the time series of hydrological elements such as rainfall and runoff. Dividing the non-stationary time series of hydrological elements into multiple sub-sequences with good consistency is one of the effective ways to adapt to the changing environment and fully apply and mine the measured data information of the model assuming that the sample obeys the "independent and identically distributed" assumption, that is, by dividing the non-stationary time series of hydrological elements, it is helpful for researchers to study the regularity of the hydrological situation. In the related technology, the hydrological element non-stationary time series division method is usually based on the stage division of the annual flood period (pre-flood period, main flood period and post-flood period), and does not consider the influence of long-term climate change on the annual flood period, resulting in low reliability of the hydrological element non-stationary sequence division, thereby affecting the regularity of the hydrological situation. SUMMARY

[0003] The main purpose of the embodiments of the present application is to provide a hydrological element non-stationary sequence division method, system, device and storage medium, which aims to improve the reliability of the hydrological element non-stationary sequence division.

[0004] To achieve the above-mentioned purpose, one aspect of the embodiments of the present application provides a hydrological element non-stationary sequence division method, comprising:

[0005] obtaining a first sequence, wherein the first sequence comprises a plurality of annual hydrological element index values arranged in chronological order;

[0006] constructing a plurality of time-continuous first sample sequences according to the first sequence, and performing data distribution consistency test between two first sample sequences to determine a plurality of first sub-sequences of the first sequence in which the hydrological elements are stably distributed;

[0007] obtaining daily hydrological element index values of each date under the interannual cycle corresponding to the first sub-sequence, and constructing a second sequence according to the daily hydrological element index values of each date;

[0008] constructing a plurality of time-continuous second sample sequences according to the second sequence, and performing consistency test between the second sample sequences and a target distribution function to determine a plurality of second sub-sequences of the second sequence in which the hydrological elements are stably distributed.

[0009] In some embodiments, the constructing a plurality of time-continuous first sample sequences according to the first sequence comprises the following steps:

[0010] determining a plurality of candidate division sequence numbers according to the data sequence numbers of the first sequence;

[0011] performing truncation and division on the first sequence according to each candidate division sequence number to obtain a plurality of sample sequence pairs corresponding to each candidate division sequence number, wherein each sample sequence pair comprises two first sample sequences.

[0012] In some embodiments, performing truncation and division on the first sequence according to each candidate division sequence number to obtain a plurality of sample sequence pairs corresponding to each candidate division sequence number comprises the following steps:

[0013] initializing a start sequence number as the sequence number of the first data in the first sequence and initializing an end sequence number as the candidate division sequence number plus 2;

[0014] performing a first step, wherein the first step comprises performing truncation and division on the first sequence according to the start sequence number, the candidate division sequence number and the end sequence number to obtain a corresponding sample sequence pair;

[0015] performing a second step, wherein the second step comprises updating the end sequence number as the current end sequence number plus 1 and repeating the first step until the end sequence number is updated as the sequence number of the last data in the first sequence and the first step is performed;

[0016] updating the start sequence number as the current start sequence number plus 1 and updating the end sequence number as the candidate division sequence number plus 2, and repeating the first step to the second step until the start sequence number is updated as the candidate division sequence number minus 1 and the first step and the second step are performed.

[0017] In some embodiments, the data distribution consistency test between the two first sample sequences comprises the following steps:

[0018] calculating the data sample difference degree between the two first sample sequences in each sample sequence pair to obtain a consistency test set;

[0019] filtering out sample sequence pairs with a data sample difference degree less than a first preset value from the consistency test set, and obtaining a division sequence number set according to the division point descriptions corresponding to the filtered sample sequence pairs, wherein the division point descriptions comprise a start sequence number, a candidate division sequence number and an end sequence number;

[0020] dividing the first sequence according to the division sequence number set to obtain a plurality of first sub-sequences.

[0021] In some embodiments, the dividing the first sequence according to the set of division indexes to obtain a plurality of first sub-sequences comprises the following steps:

[0022] According to the division point description in each division index in the set of division indexes, determine a plurality of different division schemes for the complete continuous division of the first sequence;

[0023] According to each division scheme, divide the first sequence, and calculate the sum of the data sample difference degrees between each sub-sequence after division;

[0024] The division scheme corresponding to the minimum sum of data sample difference degrees is taken as the target division scheme;

[0025] According to the target division scheme, the first sequence is divided to obtain a plurality of first sub-sequences.

[0026] In some embodiments, the first sub-sequence corresponding to the annual cycle is obtained. The daily hydrological element index value of each date is obtained, and the second sequence is constructed according to the daily hydrological element index value of each date, comprising the following steps:

[0027] Obtain the daily hydrological element index data set of each year in the annual cycle corresponding to the first sub-sequence;

[0028] Extract the daily hydrological element index value of the target date from the daily hydrological element index data set of each year, and take the average value of the extracted daily hydrological element index values from each year as the daily hydrological element index value of the target date;

[0029] According to the daily hydrological element index value of different target dates, a second sequence is constructed.

[0030] In some embodiments, the second sequence is constructed according to the second sequence, and the consistency test between the second sample sequence and the target distribution function is performed to determine the second sub-sequence of the second sequence in which the hydrological element is stationary distributed, comprising the following steps:

[0031] The first division index is initialized as the sequence number of the first data in the second sequence;

[0032] The third step is executed, wherein the third step comprises sequentially dividing the second sequence from the starting point of the first division index and according to different sequence lengths to obtain a plurality of corresponding second sample sequences, calculating the data distribution difference degree between the second sample sequence and the target distribution function, screening out the second sample sequence with a data distribution difference degree less than a second preset value, and determining the target division point of the second sequence according to the end point of the second sample sequence with the maximum sequence length screened out.

[0033] updating the first division number as the target division point, and repeating the third step until the first division number is updated as the sequence number of the last data in the second sequence;

[0034] dividing the second sequence according to the plurality of target division points to obtain a plurality of second sub-sequences.

[0035] To achieve the above object, another aspect of the embodiment of the present application provides a hydrological element non-stationary sequence division system, comprising:

[0036] A first module is configured to acquire a first sequence, wherein the first sequence comprises a plurality of annual hydrological element index values arranged in chronological order;

[0037] A second module is configured to construct a plurality of time-continuous first sample sequences according to the first sequence, and perform data distribution consistency test between two first sample sequences to determine a plurality of first sub-sequences of the first sequence in which the hydrological element is in a stationary distribution.

[0038] A third module is configured to acquire daily hydrological element index values of each date under an inter-annual cycle corresponding to the first sub-sequence, and construct a second sequence according to the daily hydrological element index values of each date.

[0039] A fourth module is configured to construct a plurality of time-continuous second sample sequences according to the second sequence, and perform consistency test between the second sample sequences and a target distribution function to determine a plurality of second sub-sequences of the second sequence in which the hydrological element is in a stationary distribution.

[0040] To achieve the above object, another aspect of the embodiment of the present application provides an electronic device, which comprises a memory, a processor, a program stored in the memory and executable on the processor, and a data bus for realizing connection communication between the processor and the memory, and the program is executed by the processor to realize the method of the above embodiment.

[0041] To achieve the above object, another aspect of the embodiment of the present application provides a storage medium, which is a computer readable storage medium for computer readable storage, and the storage medium stores one or more programs, and the one or more programs can be executed by one or more processors to realize the method of the above embodiment.

[0042] The hydrological element non-stationary sequence division method, system, device and storage medium provided by the application first obtain a plurality of annual hydrological element index values arranged in chronological order, that is, a first sequence, then construct a plurality of time-continuous first sample sequences according to the first sequence, and perform data distribution consistency test between two first sample sequences, so as to determine a plurality of first sub-sequences of the hydrological element stationary distribution in the first sequence, and realize the stationary division on the interannual scale. Then, the daily hydrological element index values of each date under the interannual cycle corresponding to the first sub-sequence are obtained, and a second sequence is constructed according to the daily hydrological element index values of each date, a plurality of time-continuous second sample sequences are constructed according to the second sequence, and consistency test is performed between the second sample sequences and the target distribution function, so as to determine a plurality of second sub-sequences of the hydrological element stationary distribution in the second sequence. The application divides the non-stationary sequence in the year on the basis of the stationary division on the interannual scale, considers the influence of long-term climate change on the annual flood season, and improves the reliability of the hydrological element non-stationary sequence division. BRIEF DESCRIPTION OF DRAWINGS

[0043] Figure 1 is a flowchart of the hydrological element non-stationary sequence division method provided by the embodiment of the application;

[0044] Figure 2 is a first sequence diagram provided by the embodiment of the application about the formation of the annual runoff of a hydrological station;

[0045] Figure 3 is an interannual characteristic index sequence provided by the embodiment of the application with different starting points of years a calculation result diagram;

[0046] Figure 4 is a second sequence diagram provided by the embodiment of the application about the daily average runoff of a hydrological station formed under different interannual cycles;

[0047] Figure 5 is a hardware structure diagram of an electronic device provided by the embodiment of the application. DETAILED DESCRIPTION

[0048] In order to make the purpose, technical scheme and advantages of the application more clear, the application will be further described in detail below in combination with the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the application, and are not used to limit the application.

[0049] It is to be understood that even though a function module is described in the system, and the logical order is shown in the flowchart, in some cases, the steps shown or described can be performed in an order different from the order shown in the module in the system or the flowchart. The terms "first", "second", and the like in the specification and claims and the above-described drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence.

[0050] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used in the specification herein is for the purpose of describing the embodiments of the present application only and is not intended to be limiting of the present application.

[0051] First, the meanings of several terms involved in the present application are explained:

[0052] K-S test (Kolmogorov-Smirnov test) is a non-parametric statistical test method, which is used to compare sample data with reference distribution (single sample K-S test) or compare the distribution of two samples (two sample K-S test). Its core idea is to judge whether the sample comes from a certain distribution or two samples come from the same distribution by comparing the maximum difference between cumulative distribution functions (CDF). The test statistic of K-S test is D value, which represents the maximum vertical distance between two cumulative distribution functions. According to the D value and the sample size, the p value is calculated. Assuming the null hypothesis is that two samples come from the same distribution, if the p value is less than the significance level (such as 0.05), the null hypothesis is rejected; if the p value is greater than the significance level, the null hypothesis cannot be rejected.

[0053] The embodiments of the present application provide a hydrological element non-stationary sequence division method, system, device and storage medium, which aims to improve the reliability of hydrological element non-stationary sequence division.

[0054] The hydrological element non-stationary sequence division method, system, device and storage medium provided by the embodiments of the present application are specifically explained by the following embodiments. First, the hydrological element non-stationary sequence division method in the embodiments of the present application is described.

[0055] The hydrological element non-stationary sequence division method provided by the embodiments of the present application relates to the technical field of data processing. The hydrological element non-stationary sequence division method provided by the embodiments of the present application can be applied to a terminal, can be applied to a server end, and can also be software running in a terminal or a server end. In some embodiments, the terminal can be a smart phone, a tablet computer, a notebook computer, a desktop computer, or the like; the server end can be configured as a stand-alone physical server, can be configured as a server cluster or a distributed system formed by multiple physical servers, can also be configured as a cloud server providing basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, CDNs, and big data and artificial intelligence platforms; and the software can be an application that implements the hydrological element non-stationary sequence division method, but is not limited to the above forms.

[0056] The present application can be used in many general or special computer system environments or configurations. For example: personal computers, server computers, handheld devices or portable devices, tablet devices, multi-processor systems, microprocessor-based systems, set-top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, distributed computing environments including any of the above systems or devices, and the like. The present application can be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, and the like that perform specific tasks or implement specific abstract data types. The present application can also be practiced in a distributed computing environment in which tasks are performed by remote processing devices connected by a communication network. In a distributed computing environment, program modules can be located in local and remote computer storage media, including storage devices.

[0057] Figure 1 The method in the optional flowchart of the hydrological element non-stationary sequence division method provided by the embodiments of the present application, Figure 1 The method in the optional flowchart of the hydrological element non-stationary sequence division method provided by the embodiments of the present application,

[0058] Step S101, obtaining a first sequence, wherein the first sequence includes a plurality of annual hydrological element index values arranged in chronological order;

[0059] Step S102, constructing a plurality of time-continuous first sample sequences according to the first sequence, and performing data distribution consistency test between two first sample sequences to determine a plurality of first sub-sequences of the hydrological element stationary distribution in the first sequence;

[0060] Step S103, obtaining the daily hydrological element index value of each date under the inter-annual cycle corresponding to the first sub-sequence, and constructing a second sequence according to the daily hydrological element index value of each date;

[0061] Step S104, constructing a plurality of time-continuous second sample sequences according to the second sequence, and performing consistency test between the second sample sequences and the target distribution function to determine a plurality of second sub-sequences of the second sequence in which the hydrological element is in a stationary distribution.

[0062] The steps S101 to S104 shown in the embodiments of the present application first obtain a plurality of annual hydrological element index values arranged in chronological order, i.e. a first sequence, then construct a plurality of time-continuous first sample sequences according to the first sequence, and perform data distribution consistency test between two first sample sequences to determine a plurality of first sub-sequences of the first sequence in which the hydrological element is in a stationary distribution, thereby achieving stationary division on the inter-annual scale. Then, the daily hydrological element index value of each date under the inter-annual cycle corresponding to the first sub-sequence is obtained, and a second sequence is constructed according to the daily hydrological element index value of each date, a plurality of time-continuous second sample sequences are constructed according to the second sequence, and consistency test is performed between the second sample sequences and the target distribution function, thereby determining a plurality of second sub-sequences of the second sequence in which the hydrological element is in a stationary distribution. The present application divides the non-stationary sequence in the year on the basis of the stationary division on the inter-annual scale, considers the influence of long-term climate change on the intra-annual flood season, and improves the reliability of the hydrological element non-stationary sequence division.

[0063] In step S101 of some embodiments, the first sequence includes a plurality of annual hydrological element index values arranged in chronological order. The type of hydrological element index can be set according to actual analysis requirements, for example, the type of hydrological element index can be rainfall index or runoff index, etc., and the embodiments of the present application do not make specific limitations.

[0064] Further, the purpose of the hydrological element non-stationary sequence division is to divide the inter-annual cycle and the intra-annual cycle. The embodiments of the present application can select data of an index type convenient for analysis as the first sequence to be analyzed.

[0065] Illustratively, the preliminary measured runoff and rainfall data are used as data, the sequence length of the measured runoff data is , the sequence length of the measured rainfall data is , the unit is year, and the selection method is as follows:

[0066] When ≥ preset data amount (for example, 30), the runoff data is used as subsequent analysis data, i.e. the annual runoff of each year forms a first sequence.

[0067] When <When the preset data volume is used, if ≥Preset data volume, using rainfall data as the subsequent analysis data, that is, forming the first sequence based on the annual rainfall of each year; if <The preset data volume allows for the simultaneous use of runoff and rainfall data as subsequent analysis data, thus forming a series based on annual rainfall for each year. and the annual runoff sequence The two sequences are used as the first sequence to divide the interannual and intra-annual cycles, respectively, and finally the cycle division results based on the two sequences are combined.

[0068] In step S102 of some embodiments, the first sequence can be pre-divided into various first sample sequences. The annual hydrological element index values ​​contained in the first sample sequences are continuous in time, and different first sample sequences differ in length, time start point, or time end point. By performing a data distribution consistency test on all pairs of first sample sequences, two first sample sequences with significantly different data distributions can be identified, indicating that the data distribution within these two first sample sequences may be stationary. Based on each pair of two first sample sequences with significantly different data distributions, the non-stationary first sequence of hydrological elements can be continuously divided, resulting in multiple first sub-sequences with stationary hydrological element distributions within the sequence, thereby achieving the division of interannual cycles.

[0069] In this embodiment, the method for verifying the consistency of the data distribution of the two first sample sequences can be the two-sample KS test, the Anderson-Darling test, the Epps-Singleton test, or the maximum mean difference method, etc.

[0070] In step S103 of some embodiments, each first subsequence corresponds to an interannual cycle, and the ranges of the interannual cycles do not overlap. Daily hydrological element index values ​​for each date within the interannual cycle are obtained, and a second sequence is constructed based on these daily hydrological element index values. For example, assuming an interannual cycle is 1997-2002, the daily hydrological element index values ​​for each date in any year within 1997-2002 can be obtained, thus obtaining a second sequence with dates as the time unit. In another example, assuming an interannual cycle is 1997-2002, the daily hydrological element index values ​​for each date corresponding to these six years within 1997-2002 can be obtained. For each date (month and day), the average daily hydrological element index value for that date across the six years can be calculated, and the second sequence is constructed using the average values ​​of all dates.

[0071] In step S104 of some embodiments, the second sequence can be pre-divided into various second sample sequences, the diurnal hydrological element index values contained in the second sample sequences are continuous in time, and different second sample sequences are different in length, time starting point or time ending point. By respectively testing the data distribution consistency of each second sample sequence with the target distribution function, the second sample sequence that is obviously consistent with the data distribution of the target distribution function can be determined, indicating that the second sample sequence is consistent with the data distribution of the target distribution function. According to each second sample sequence that is obviously consistent with the data distribution of the target distribution function, the second sequence can be continuously divided to obtain multiple second sub-sequences with stable distribution of hydrological elements in the sequence, thereby realizing the division of the intra-annual cycle.

[0072] In the present embodiment, the target distribution function can be selected according to the regular variation rule of the intra-annual cycle of the hydrological element index. The target distribution function can adopt a Poisson distribution function. The Poisson distribution curve is a graphical representation of the Poisson distribution probability density function (PMF) or cumulative distribution function (CDF). The Poisson distribution is a discrete probability distribution used to describe the probability distribution of the number of events occurring within a fixed time or space. It is suitable for cases where the probability of event occurrence is low and independent.

[0073] In the present embodiment, the method for testing the data distribution consistency of the second sample sequence with the Poisson distribution function can adopt a two-sample K-S test method, an Anderson-Darling test method, an Epps-Singleton test method or a maximum mean difference method, etc.

[0074] According to some embodiments of the present application, in step S102, the step of constructing multiple time-continuous first sample sequences according to the first sequence can include but is not limited to the following steps:

[0075] Step S201, determining multiple candidate division sequence numbers according to the data sequence numbers of the first sequence;

[0076] Step S202, according to each candidate division sequence number, intercepting and dividing the first sequence to obtain multiple sample sequence pairs corresponding to each candidate division sequence number, wherein each sample sequence pair includes two first sample sequences.

[0077] In the embodiment, the data sequence number of the first sequence is used to determine the corresponding candidate division sequence number. For example, assuming that the first sequence is {40, 34, 23, 31, 24}, the candidate division sequence numbers are 0, 1, 2, 3, and 4 respectively. For each candidate division sequence number, at least one data before the candidate division sequence number and including the candidate division sequence number in the first sequence is taken to form a first sample sequence, and at least one data after the candidate division sequence number in the first sequence is taken to form a first sample sequence. For example, for the candidate division sequence number "2", one data before the sequence number 2 and the sequence number 2 in the first sequence is taken to obtain the first sample sequence {23}, two data before the sequence number 2 and the sequence number 2 in the first sequence is taken to obtain the first sample sequence {34, 23}, three data before the sequence number 2 and the sequence number 2 in the first sequence is taken to obtain the first sample sequence {40, 34, 23}, one data after the sequence number 2 in the first sequence is taken to obtain the first sample sequence {31}, and so on, thereby obtaining multiple sample sequence pairs corresponding to each candidate division sequence number.

[0078] According to some embodiments of the present application, the step of obtaining multiple sample sequence pairs corresponding to each candidate division sequence number by intercepting and dividing the first sequence according to each candidate division sequence number in step S202 can include but is not limited to the following steps:

[0079] In step S301, the starting sequence number is initialized as the sequence number of the first data in the first sequence, and the ending sequence number is initialized as the candidate division sequence number plus 2.

[0080] In step S302, a first step is performed, wherein the first step includes intercepting and dividing the first sequence according to the starting sequence number, the candidate division sequence number, and the ending sequence number to obtain the corresponding sample sequence pair.

[0081] In step S303, a second step is performed, wherein the second step includes updating the ending sequence number as the current ending sequence number plus 1, and repeatedly performing the first step until the ending sequence number is updated as the sequence number of the last data in the first sequence and the first step is performed.

[0082] In step S304, the starting sequence number is updated as the current starting sequence number plus 1, and the ending sequence number is updated as the candidate division sequence number plus 2, and the first step and the second step are repeatedly performed until the starting sequence number is updated as the candidate division sequence number minus 1 and the first step and the second step are performed.

[0083] In the embodiment, the number of the first data in the first sequence is 0. For each candidate division sequence, i.e., the year to be analyzed , the first sample sequence is extracted as follows:

[0084] ​S11, initializing a start sequence number x=0 (i.e. the sequence number of the first data in the first sequence), and initializing an end sequence number y=n+2;

[0085] S12, extracting data with the first sequence number from x to n, and extracting data with the first sequence number from n+1 to y, to obtain two corresponding first sample sequences, which are a sample sequence pair;

[0086] S13, setting y=y+1, and repeating step S12 until step S12 with the condition of y=N is executed, and then executing step S14; in the above steps S11 to S13, the sequence from x to n in the first sequence is extracted, and the sequence from n+1 to y in the first sequence is extracted, respectively, and the sample sequence pair is obtained correspondingly;

[0087] S14, setting x=x+1 and y=n+2, and then repeating the above steps S12 to S13 until step S14 with the condition of x=n-1 is executed, and then ending the current process, and traversing to obtain multiple sample sequence pairs of candidate division sequence numbers; in step S14, the distance from x to n-1 is shortened by 1 as a step, i.e. the sequence from x to n-1 in the first sequence is extracted, and the sequence from n+1 to y in the first sequence is extracted, respectively, and the sample sequence pair is obtained correspondingly.

[0088] In the above traversal process, the first sequence is divided into two samples before and after , and the corresponding annual hydrological element index value is attributed to the former sample.

[0089] According to some embodiments of the present application, in step S102, the step of performing data distribution consistency test between the two first sample sequences to determine multiple first sub-sequences of the hydrological element stable distribution in the first sequence can include but is not limited to the following steps:

[0090] Step S401, calculating the data sample difference degree between the two first sample sequences in each sample sequence pair to obtain a consistency test set;

[0091] Step S402, screening out the sample sequence pair with the data sample difference degree less than the first preset value from the consistency test set, and obtaining a division sequence number set according to the division point description corresponding to the screened sample sequence pair, wherein the division point description includes a start sequence number, a candidate division sequence number and an end sequence number;

[0092] ​​​​​​​​​​Step S403, dividing the first sequence according to the division sequence number set to obtain a plurality of first sub-sequences.

[0093] In the embodiment, the two-sample K-S test method can be used to perform consistency test on the data distribution of the two first sample sequences in the sample sequence pair. The sample consistency can be measured by the data sample difference degree, that is, the progressive significance (two-tailed, denoted as ) between the two samples can be calculated. The consistency test set is formed by calculating the value of each sample sequence pair.

[0094] From the consistency test set, sample sequence pairs with less than a first preset value (for example, ) are screened out. The sample sequence pair with indicates that the difference between the corresponding two samples is large, and they may belong to different interannual cycles. According to the division point description (the division point description includes the start sequence number, the candidate division sequence number and the end sequence number) corresponding to all the screened sample sequence pairs, the division sequence number set is obtained. The continuous division scheme from to in the first division can be determined by the division point description in the division sequence number set, so as to obtain M first sub-sequences .

[0095] According to some embodiments of the present application, in step S403, the step of dividing the first sequence according to the division sequence number set to obtain a plurality of first sub-sequences can include but is not limited to the following steps:

[0096] Step S501, determining a plurality of different division schemes for the complete and continuous division of the first sequence according to the division point description in the division sequence number set;

[0097] Step S502, dividing the first sequence according to each division scheme and calculating the sum of the data sample difference degrees between the divided sub-sequences;

[0098] Step S503, taking the division scheme corresponding to the sum of the data sample difference degrees as the target division scheme;

[0099] Step S504, dividing the first sequence according to the target division scheme to obtain a plurality of first sub-sequences.

[0100] In the embodiment, the possible multiple different continuous division schemes can be described by dividing each division point in the division number set. For example, assuming that the time period of the first sequence is from 1962 to 2024, and by the process of the above embodiment, all division point descriptions in the division number set can correspond to division time periods of 1962-1994, 1995-1996, 1997-2000, 1988-1996, 2001-2024, 1962-1987, 1997-2003, and 2004-2024. According to the principle of complete and continuous division of the first sequence, the above division time periods can be combined to obtain division scheme 1: 1962-1994, 1995-1996, 1997-2000, and 2001-2024, and division scheme 2: 1962-1987, 1988-1996, 1997-2003, and 2004-2024.

[0101] According to each division scheme, the first sequence is divided, and the sum of the data sample difference degrees between each sub-sequence after division is calculated, that is, the sum of the data sample difference degrees between the first sample sequences corresponding to each division time period in the division scheme. Then, the sum of the data sample difference degrees between each is calculated. The division scheme corresponding to the sum of the data sample difference degrees is taken as the target division scheme, and the first sequence is divided by using the target division scheme to obtain multiple first sub-sequences. In the embodiment of the present application, the division scheme corresponding to the sum of the data sample difference degrees is taken as the target division scheme, which can maximize the difference degrees between the first sub-sequences and improve the rationality of the interannual division.

[0102] Further, the division scheme with the largest number of division time periods and the smallest sum of the data sample difference degrees can be selected as the target division scheme, which can improve the rationality of the interannual division and the precision of the interannual period division.

[0103] According to some embodiments of the present application, in step S103, the daily hydrological element index values of each date in the interannual period corresponding to the first sub-sequence are obtained, and the second sequence is constructed according to the daily hydrological element index values of each date. The step can include but is not limited to the following steps:

[0104] In step S601, the daily hydrological element index data set of each year in the interannual period corresponding to the first sub-sequence is obtained.

[0105] In step S602, the daily hydrological element index value of the target date is extracted from the daily hydrological element index data set of each year, and the average value of the extracted daily hydrological element index values from each year is taken as the daily hydrological element index value of the target date.

[0106] ​​​Step S603, constructing a second sequence according to the daily hydrological element index values of different target dates.

[0107] In this embodiment, the daily hydrological element index values of each date in a year are included in the second sequence. Taking January 1st as an example, the calculation process of the corresponding daily hydrological element index value (such as daily average rainfall) in the second sequence of the interannual cycle from 1999 to 2002 is as follows: obtaining the rainfall of January 1st in 1999, the rainfall of January 1st in 2000, the rainfall of January 1st in 2001, and the rainfall of January 1st in 2002, and calculating the average of the above four rainfall values to obtain the daily average rainfall of January 1st in the second sequence.

[0108] According to some embodiments of the present application, in step S104, the step of constructing a plurality of time-continuous second sample sequences from the second sequence and performing consistency test between the second sample sequences and the target distribution function to determine a plurality of second sub-sequences of the second sequence in which the hydrological element is in a stationary distribution can include but is not limited to the following steps:

[0109] Step S701, initializing the first division sequence number as the sequence number of the first data in the second sequence;

[0110] Step S702, performing the third step, wherein the third step includes dividing the second sequence in sequence according to different sequence lengths starting from the first division sequence number to obtain a plurality of corresponding second sample sequences, calculating the data distribution difference degree between the second sample sequence and the target distribution function, screening out the second sample sequence with a data distribution difference degree less than a second preset value, and determining the target division point of the second sequence according to the end point of the second sample sequence with the maximum sequence length screened out;

[0111] Step S703, updating the first division sequence number to the target division point and repeating the third step until the first division sequence number is updated to the sequence number of the last data in the second sequence;

[0112] Step S704, dividing the second sequence according to a plurality of target division points to obtain a plurality of second sub-sequences.

[0113] In this embodiment, the starting date sequence number in the second sequence is 0, and the date sequence number of the day to be analyzed is The second sample data is divided, and the data distribution difference degree between each second sample sequence and the Poisson distribution (i.e. the target distribution function) is calculated in sequence using the single-parameter K-S test, i.e. the asymptotic significance (two-tailed, denoted as ), the calculation process is as follows:

[0114] The first division sequence number (initialized as the starting date sequence number) in the second sequence is extracted in sequence to Given multiple second sample sequences, calculate the sum of the values ​​of each second sample sequence and the Poisson distribution. value;

[0115] Filter out The value is greater than the second preset value (e.g.) The second sample sequence corresponding to) Remember the largest one The corresponding serial number is ;

[0116] Update the first partition number to Extract sequentially To each , , ..., Multiple second sample sequences Calculate the relationship between each second sample sequence and the Poisson distribution. Value, select one of them The value is greater than the second preset value (e.g.) The termination date corresponding to the second sample sequence of ) is denoted as the sequence number corresponding to the longest termination date. ;

[0117] Update the first partition number to Extract sequentially To each , , ..., Multiple second sample sequences Calculate the relationship between each second sample sequence and the Poisson distribution. Value, select one of them The value is greater than the second preset value (e.g.) The termination date corresponding to the second sample sequence of ) is denoted as the sequence number corresponding to the longest termination date. ;

[0118] This process continues until the first partition index is updated to the index of the last data in the second sequence, resulting in a dynamically partitioned set within the year. The set is dynamically divided into multiple target partitioning points within the year, and each target partitioning point is the sequence number corresponding to the maximum termination date obtained from each step of the above process. , , .... Divide points based on multiple targets. The second sequence is divided into multiple second subsequences, thus corresponding to multiple intra-annual cycles, i.e. The second sequence is divided into two samples, with the boundary set before and after. The corresponding daily hydrological element index values ​​are assigned to the previous sample.

[0119] According to some embodiments of the present application, only the runoff or rainfall index type is analyzed, and the interannual cycle is divided according to the first sequence of one of the annual runoff or annual rainfall data , and the intra-annual cycle is divided according to the second sequence of the daily runoff or daily rainfall data. .

[0120] If runoff data and rainfall index types are analyzed at the same time, the final annual cycle can be determined by combining the division results of the interannual cycle corresponding to the first sequence of the two types of data, and the final intra-annual cycle can be determined by combining the division results of the intra-annual cycle corresponding to the second sequence of the two types of data. For example, the interannual cycle of runoff , the interannual cycle of rainfall ; merging and deduplicating the two gives ; deleting the runoff data division node that appears in the rainfall division node, gives , which is the interannual cycle division node.

[0121] Next, taking the measured runoff data of a certain river station from 1950 to 2021 as an example, the calculation process of the method according to the embodiments of the present application is described.

[0122] First, the data is selected, and the length of the measured runoff data of the hydrological station from 1950 to 2021 is =72>30, and the runoff data is used as the data to carry out dynamic division.

[0123] Second, the interannual characteristic index (i.e., the annual hydrological element index value) is calculated to obtain the first sequence, i.e., the first sequence of the annual average runoff of the hydrological station , and the first sequence is shown in Figure 2 .

[0124] Third, interannual scale consistency test:

[0125] Using two-sample K-S test to traverse the year to be analyzed , the of the interannual characteristic index sequence (i.e., the first sample sequence) starting from 1950, 1998, and 2002 is calculated , and the calculation results are shown in Figure 3 .

[0126] Fourth, interannual cycle division:

[0127] After screening, the scheme that satisfies and can realize continuous division from the starting point to the end point of the interannual characteristic index sequence is shown in Table 1.

[0128] Table 1 Interannual cycle division scheme

[0129]

[0130] As shown in Table 1, Scheme 1 is the interannual cycle division scheme that meets the requirements, and the 1950–2021 runoff sequence is divided into... =1950-1998 =1999—2002, =2003—2021 three first subsequences.

[0131] Fifth, calculation of characteristic indicators for the year: calculate separately. , and The annual characteristic index values ​​(i.e., daily hydrological element index values) are used to obtain the corresponding second sequence. ,like Figure 4 As shown.

[0132] Sixth, dynamic division within the year:

[0133] Traversal using a single-parameter KS test , and Date to be analyzed calculate The results of the dynamic classification within the year after screening are shown in Table 2.

[0134] Table 2. Intra-year Period Division Scheme

[0135]

[0136] Seventh, dynamic division of interannual and intraannual scales:

[0137] The runoff sequence of the hydrological station from 1950 to 2021 contains three interannual cycles. The subsequence from 1950 to 1998 has 9 intra-annual periods, the subsequence from 1999 to 2002 has 11 intra-annual periods, and the subsequence from 2003 to 2021 has 8 intra-annual periods. In practical applications, the interannual cycle of inflow can be determined based on the characteristic parameters of the inflow flood forecast, and then measured data can be used and extracted according to the corresponding intra-annual division scheme.

[0138] This application also proposes a system for classifying non-stationary sequences of hydrological elements, including:

[0139] The first module is used to obtain a first sequence, wherein the first sequence includes multiple annual hydrological element index values ​​arranged in chronological order.

[0140] The second module is used to construct multiple time-continuous first sample sequences based on the first sequence, and to perform a data distribution consistency test between two first sample sequences in order to determine multiple first subsequences in the first sequence with a stable distribution of hydrological elements.

[0141] The third module is configured to obtain the daily hydrological element index values of each date under the inter-annual cycle corresponding to the first sub-sequence, and construct a second sequence according to the daily hydrological element index values of each date;

[0142] The fourth module is configured to construct a plurality of second sample sequences which are continuous in time according to the second sequence, and perform consistency test between the second sample sequences and a target distribution function to determine a plurality of second sub-sequences of the second sequence in which the hydrological element is in a stationary distribution.

[0143] It can be understood that the contents in the above hydrological element non-stationary sequence division method embodiments are applicable to the present system embodiment, the present system embodiment specifically implements the same functions as the above hydrological element non-stationary sequence division method embodiments, and achieves the same beneficial effects as the above hydrological element non-stationary sequence division method embodiments.

[0144] The present application also provides an electronic device, which comprises a memory, a processor, a program stored in the memory and executable on the processor, and a data bus for realizing connection and communication between the processor and the memory, and the program is executed by the processor to realize the above hydrological element non-stationary sequence division method. The electronic device can be any intelligent terminal including a tablet computer, a vehicle-mounted computer, etc.

[0145] Please refer to Figure 5 , Figure 5 The hardware structure of the electronic device of another embodiment is illustrated, which comprises:

[0146] The processor 501 can be implemented in the form of a general-purpose CPU (Central Processing Unit), a microprocessor, an ASIC (Application Specific Integrated Circuit), or one or more integrated circuits, etc., and is used to execute related programs to realize the technical solutions provided by the present application;

[0147] The memory 502 can be implemented in the form of a ROM (ReadOnly Memory), a static storage device, a dynamic storage device, or a RAM (Random Access Memory), etc. The memory 502 can store an operating system and other application programs, and when the technical solutions provided by the present application are implemented by software or firmware, the related program codes are saved in the memory 502 and are called and executed by the processor 501 to realize the hydrological element non-stationary sequence division method of the present application;

[0148] The input / output interface 503 is configured to realize information input and output.

[0149] The communication interface 504 is configured to realize communication interaction between the device and other devices, and the communication can be realized in a wired manner (for example, a USB, a network cable, or the like) or in a wireless manner (for example, a mobile network, WIFI, Bluetooth, or the like).

[0150] The bus 505 is configured to transmit information between various components (for example, the processor 501, the memory 502, the input / output interface 503, and the communication interface 504) of the device.

[0151] The processor 501, the memory 502, the input / output interface 503, and the communication interface 504 are connected to each other in the device through the bus 505.

[0152] The embodiment of the present application further provides a storage medium, the storage medium is a computer readable storage medium, and the storage medium is used for computer readable storage. The storage medium stores one or more programs, and the one or more programs can be executed by one or more processors to implement the above-mentioned hydrological element non-stationary sequence division method.

[0153] The memory is a non-transient computer readable storage medium, and can be used to store non-transient software programs and non-transient computer executable programs. In addition, the memory can include a high-speed random access memory, and can further include a non-transient memory, for example, at least one magnetic disk storage device, a flash memory device, or other non-transient solid-state memory device. In some embodiments, the memory can optionally include a memory remotely arranged relative to the processor, and the remote memory can be connected to the processor through a network. Examples of the network include but are not limited to the Internet, an intranet, a local area network, a mobile communication network, and a combination thereof.

[0154] The embodiments described in the embodiments of the present application are used to more clearly illustrate the technical solutions of the embodiments of the present application, and do not constitute a limitation on the technical solutions provided by the embodiments of the present application. Those skilled in the art can know that, with the evolution of technology and the appearance of new application scenarios, the technical solutions provided by the embodiments of the present application are also applicable to similar technical problems.

[0155] Those skilled in the art can understand that the technical solutions shown in the figure do not constitute a limitation on the embodiments of the present application, and can include more or fewer steps than the figure, or combine certain steps, or different steps.

[0156] The system embodiments described above are merely illustrative, wherein the units described as separate components can or can not be physically separated, i.e., can be located in one place, or can be distributed to multiple network units. Part or all of the modules can be selected according to actual needs to achieve the purpose of the embodiment.

[0157] Those skilled in the art can understand that all or some steps in the method disclosed above, and the functional modules / units in the system, the device can be implemented as software, firmware, hardware and appropriate combinations thereof.

[0158] The terms "first", "second", "third", "fourth" and the like in the description of the application and in the claims of the foregoing drawings, if any, are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the use of the terms so

[0159] It should be understood that in this application, "at least one" means one or more, and "multiple" means two or more. "And / or" is used to describe the relationship between the associated objects, which means that there can be three relationships, for example, "A and / or B" can mean that there are three cases: only A, only B, and A and B exist at the same time, where A and B can be singular or plural. The character " / " generally represents that the associated objects before and after are in an "or" relationship. "At least one of the following" or similar expressions means any combination of these items, including any combination of single or multiple items. For example, at least one of a, b or c, can mean: a, b, c, "a and b", "a and c", "b and c", or "a and b and c", where a, b, and c can be single or multiple.

[0160] In several embodiments provided in the present application, it should be understood that the disclosed system and method can be implemented in other manners. For example, the system embodiments described above are merely schematic. For example, the division of the units is only a logical function division. There can be another division manner for the actual implementation, for example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the displayed or discussed mutual couplings or direct couplings or communication connections between different units, can be indirect couplings or communication connections through some interfaces, and electrical, mechanical or other forms.

[0161] The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they can be located in one place, or can be distributed on multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the embodiment.

[0162] In addition, the functional units in each embodiment of the present application can be integrated in one processing unit, or each unit can be physically present separately, or two or more units can be integrated in one unit. The integrated unit can be realized in the form of hardware or in the form of a software functional unit.

[0163] If the integrated unit is realized in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer readable storage medium. Based on this understanding, the technical solutions of the present application essentially or the part that contributes to the prior art or the whole or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, and includes multiple instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods in the various embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, and various program storage media.

[0164] The preferred embodiments of the present application are described above with reference to the accompanying drawings, but this does not limit the scope of the present application. Any modifications, equivalent replacements and improvements made by those skilled in the art without departing from the scope and spirit of the present application should be within the scope of the present application.

Claims

1. A method for partitioning non-stationary sequences of hydrological elements, characterized in that, Includes the following steps: Obtain a first sequence, wherein the first sequence includes multiple annual hydrological element index values ​​arranged in chronological order; Multiple time-continuous first sample sequences are constructed based on the first sequence, and a data distribution consistency test is performed between two first sample sequences to determine multiple first subsequences in which hydrological elements are stably distributed in the first sequence. Obtain the daily hydrological element index values ​​for each date under the interannual cycle corresponding to the first subsequence, and construct a second sequence based on the daily hydrological element index values ​​for each date; Based on the second sequence, multiple time-continuous second sample sequences are constructed, and a consistency test is performed between the second sample sequences and the target distribution function to determine multiple second subsequences in which hydrological elements are stationary in the second sequence; The step of constructing multiple time-continuous second sample sequences based on the second sequence, and performing a consistency test between the second sample sequences and the target distribution function to determine multiple second subsequences with stationary distribution of hydrological elements in the second sequence, includes the following steps: Initialize the first partition index to the index of the first data in the second sequence; The third step includes dividing the second sequence sequentially according to different sequence lengths, starting from the first division number, to obtain multiple corresponding second sample sequences. The third step also includes calculating the degree of data distribution difference between the second sample sequences and the target distribution function; selecting second sample sequences with a data distribution difference degree less than a second preset value; and determining the target division point of the second sequence based on the endpoint of the second sample sequence with the longest selected sequence. Update the first partition number to the target partition point, and repeat the third step until the first partition number is updated to the number of the last data in the second sequence; The second sequence is divided according to multiple target dividing points to obtain multiple second subsequences.

2. The method for partitioning non-stationary sequences of hydrological elements according to claim 1, characterized in that, The step of constructing multiple time-continuous first sample sequences based on the first sequence includes the following steps: Multiple candidate partition numbers are determined based on the data sequence number of the first sequence; Based on each candidate partition number, the first sequence is truncated and divided to obtain multiple sample sequence pairs corresponding to each candidate partition number, wherein each sample sequence pair includes two first sample sequences.

3. The method for dividing non-stationary sequences of hydrological elements according to claim 2, characterized in that, The step of truncating and dividing the first sequence according to each candidate partition number to obtain multiple sample sequence pairs corresponding to each candidate partition number includes the following steps: The starting sequence number is initialized to the sequence number of the first data in the first sequence, and the ending sequence number is initialized to the candidate partition number plus 2. Perform the first step, wherein the first step includes truncating and dividing the first sequence according to the start number, the candidate partition number and the end number to obtain corresponding sample sequence pairs; Perform the second step, wherein the second step includes updating the termination number to the current termination number plus 1, and repeating the first step until the termination number is updated to the number of the last data in the first sequence and the first step is performed. Update the starting number to the current starting number plus 1, and update the ending number to the candidate partition number plus 2. Repeat the first step to the second step until the starting number is updated to the candidate partition number minus 1 and the first step and the second step are executed.

4. The method for dividing non-stationary sequences of hydrological elements according to claim 3, characterized in that, The step of performing a data distribution consistency test between two first sample sequences to determine multiple first subsequences with a stable distribution of hydrological elements in the first sequence includes the following steps: Calculate the degree of data sample difference between the two first sample sequences in each of the sample sequence pairs to obtain a consistency test set; From the consistency test set, sample sequence pairs with a difference in data sample degree less than a first preset value are selected. Based on the corresponding split point descriptions of the selected sample sequence pairs, a set of split number is obtained, wherein the split point descriptions include a start number, a candidate split number, and an end number. The first sequence is divided according to the set of partition numbers to obtain multiple first subsequences.

5. The method for partitioning non-stationary sequences of hydrological elements according to claim 4, characterized in that, The step of dividing the first sequence according to the set of partition numbers to obtain multiple first subsequences includes the following steps: Based on the description of each partition point in the partition sequence set, determine multiple different partitioning schemes that can completely and continuously partition the first sequence. The first sequence is divided according to each partitioning scheme, and the sum of the differences in data samples between each subsequence after partitioning is calculated; The partitioning scheme that minimizes the sum of the differences in the data samples is taken as the target partitioning scheme; The first sequence is divided according to the target partitioning scheme to obtain multiple first subsequences.

6. The method for partitioning non-stationary sequences of hydrological elements according to claim 1, characterized in that, The step of obtaining the daily hydrological element index values ​​for each date under the interannual cycle corresponding to the first subsequence, and constructing the second sequence based on the daily hydrological element index values ​​for each date, includes the following steps: Obtain the daily hydrological element index dataset for each year in the interannual cycle corresponding to the first subsequence; Extract the daily hydrological element index values ​​for the target date from the daily hydrological element index datasets of each year, and use the average value of the extracted daily hydrological element index values ​​from each year as the daily hydrological element index value for the target date. A second sequence is constructed based on the daily hydrological element index values ​​for different target dates.

7. A system for classifying non-stationary sequences of hydrological elements, characterized in that, include: The first module is used to obtain a first sequence, wherein the first sequence includes multiple annual hydrological element index values ​​arranged in chronological order. The second module is used to construct multiple time-continuous first sample sequences based on the first sequence, and to perform a data distribution consistency test between two first sample sequences to determine multiple first subsequences in which hydrological elements are stably distributed in the first sequence. The third module is used to obtain the daily hydrological element index values ​​of each date in the interannual cycle corresponding to the first subsequence, and to construct a second sequence based on the daily hydrological element index values ​​of each date. The fourth module is used to construct multiple time-continuous second sample sequences based on the second sequence, and to perform a consistency test between the second sample sequences and the target distribution function to determine multiple second subsequences in which hydrological elements are stationary in the second sequence; The fourth module is specifically used to perform the following steps: Initialize the first partition index to the index of the first data in the second sequence; The third step includes dividing the second sequence sequentially according to different sequence lengths, starting from the first division number, to obtain multiple corresponding second sample sequences. The third step also includes calculating the degree of data distribution difference between the second sample sequences and the target distribution function; selecting second sample sequences with a data distribution difference degree less than a second preset value; and determining the target division point of the second sequence based on the endpoint of the second sample sequence with the longest selected sequence. Update the first partition number to the target partition point, and repeat the third step until the first partition number is updated to the number of the last data in the second sequence; The second sequence is divided according to multiple target dividing points to obtain multiple second subsequences.

8. An electronic device, characterized in that, The electronic device includes a memory, a processor, a program stored in the memory and executable on the processor, and a data bus for enabling communication between the processor and the memory, wherein the program, when executed by the processor, implements the steps of the method as described in any one of claims 1 to 6.

9. A storage medium, said storage medium being a computer-readable storage medium for computer-readable storage, characterized in that, The storage medium stores one or more programs, which can be executed by one or more processors to implement the steps of the method according to any one of claims 1 to 6.

Images