A river section flow field data acquisition and analysis method and system

Abstract

The application discloses a river section flow field data collection and analysis method and system, which collects initial data of a river section flow field in a preset time period through a flow field collection system, and processes the initial data to sequentially obtain effective data of the river section flow field, interpolation data of the river section flow field and storage data of the river section flow field; the effective data of the river section flow field is analyzed through a flow field analysis system to obtain average flow velocities at different characteristic lines in a horizontal X direction and average flow velocities at different characteristic lines in a vertical Y direction, and the interpolation data of the river section flow field is analyzed to obtain spatial distribution change characteristics of the river section flow field; and the above data is dynamically displayed. The application comprehensively collects overall characteristic data of the river section flow field, automatically analyzes horizontal, vertical and spatial distribution change characteristics of the river section flow field, calculates time-space change characteristic time series of the river section flow field, dynamically displays result data in combination with a time axis, and quickly, comprehensively and visually presents time-space change characteristics of the river section flow field.

Description

TECHNICAL FIELD

[0001] The present application relates to hydrology and water resources and data processing technology, and particularly relates to a river section flow field data acquisition and analysis method and system. BACKGROUND

[0002] River section flow field is an important hydrological parameter, and it is crucial to dynamically monitor the spatio-temporal variation characteristics of river section flow field. However, the existing river section flow velocity acquisition equipment is mostly aimed at a single point, and cannot comprehensively reflect the overall characteristics of the section flow field, and the spatio-temporal variation characteristics of the river section flow field cannot be visually presented. SUMMARY

[0003] The present application aims to solve the problems existing in the prior art, and provides a river section flow field data acquisition and analysis method and system. The present application realizes the full-process integration of flow field data acquisition, processing, analysis and display through a flow field acquisition system, a data processing system, a flow field analysis system and a flow field display system, and quickly, comprehensively and visually presents the spatio-temporal variation characteristics of the river section flow field, thereby providing a scientific basis for water resource management and water ecological protection.

[0004] TECHNICAL SOLUTION The present application provides a river section flow field data acquisition and analysis method, which comprises the following steps:

[0005] Step 1, acquiring initial data [VO] of the river section flow field in a preset time period through a flow field acquisition system,

[0006] ;

[0007] ;

[0008] ;

[0009] ;

[0010] m is the total number of acquisition time series in the preset acquisition time period; VO m is the initial data of the river section flow field of the mth time series; VO (r, s) is the initial data of the river flow velocity at the coordinate (r, s) monitoring point; r is the total number of monitoring points in the river X direction; s is the total number of monitoring points in the river Y direction; W is the width of the river section; and D is the depth of the river section.

[0011] Step 2, processing the initial data [VO] of the river section flow field obtained in step 1 through a data processing system to obtain effective data [VE] of the river section flow field, interpolation data [VI] of the river section flow field and storage data [VS] of the river section flow field in sequence;

[0012]

[0013] (VI m , VE m , VI m ) is the flow field data of the mth time series of river cross section;

[0014]

[0015]

[0016] EFF is the 3σ effective analysis function; VE m is the effective data of the mth time series of river cross section flow field; VE (r, s) is the effective data of river flow velocity at the monitoring point of coordinates (r, s);

[0017]

[0018]

[0019]

[0020]

[0021] KRI is the Kriging interpolation function; VI m is the interpolation data of the mth time series of river cross section flow field; VI (j, k) is the interpolation data of river flow velocity at the monitoring point of coordinates (j, k); j is the total number of interpolation points in the X direction of the river; k is the total number of interpolation points in the Y direction of the river; ΔX is the interpolation resolution in the X direction of the river; ΔY is the interpolation resolution in the Y direction of the river;

[0022] Step 3, analyze the effective data of the river cross section flow field [VE] through the flow field analysis system, get the average flow velocity [AX] at different characteristic lines in the horizontal X direction and the average flow velocity [AY] at different characteristic lines in the vertical Y direction, and analyze the interpolation data of the river cross section flow field [VI] to get the spatial distribution change characteristics C of the river cross section flow field;

[0023] ;

[0024] ;

[0025] ;

[0026] AX m is the horizontal variation characteristic of the mth time series of river cross section flow field; AX X=r is the average flow velocity of the rth characteristic line in the horizontal direction;

[0027] ;

[0028] ;

[0029] ;

[0030] AY m is the vertical variation characteristic of the flow field of the river section for the mth time series; AY Y=s is the average flow velocity of the vertical s th characteristic line;

[0031] ;

[0032] β is the change slope of the spatial distribution time series of the flow field of the river section; Z is the significance of the change trend of the spatial distribution time series of the flow field of the river section;

[0033] Step 4, dynamically displaying the data obtained in step 3, including: dynamically presenting the horizontal variation characteristic broken line graph time series [F AX ] of the flow field of the river section in the preset collection period, the vertical variation characteristic broken line graph time series [F AY ], the spatial distribution heat map time series [H VI ] and the spatial distribution variation characteristic heat map H C of the flow field of the river section.

[0034] Further, the flow field collection system in step 1 is provided with a collection module, which collects the horizontal flow velocity and the vertical flow velocity of each coordinate point of the river section in real time through a sensing probe, and generates the initial data [VO] of the flow field of the river section in the preset collection period.

[0035] Further, the data processing system in step 2 includes an effective processing module, an interpolation processing module and a data storage module, the effective processing module performs effectiveness judgment on the initial data [VO] of the flow field of the river section to obtain the effective data [VE] of the flow field of the river section, the interpolation processing module performs spatial interpolation on the obtained effective data [VE] of the flow field of the river section to generate the interpolation data [VI] of the flow field of the river section, and the data storage module stores the above data; the specific method is:

[0036] Step 2.1, based on the initial data [VO] of the flow field of the river section, the effective processing module uses 3σ effective analysis method to judge the effectiveness of the data, eliminates invalid values and device null values, and generates the effective data [VE] of the flow field of the river section;

[0037] Step 2.2, based on the effective data [VE] of the flow field of the river section, the interpolation processing module uses Kriging interpolation method to perform spatial interpolation processing on the data, and generates the interpolation data [VI] of the flow field of the river section;

[0038] Step 2.3, the data storage module stores the initial data [VO], the valid data [VE] and the interpolation data [VI] in JSON format; and generates river section flow field storage data [VS];

[0039] ;

[0040] (VO m , VE m , VI m ) is the river section flow field storage data of the mth time series.

[0041] Further, the flow field analysis system comprises river section flow field transverse variation characteristic analysis, vertical variation characteristic analysis and spatial distribution variation characteristic analysis. The transverse variation characteristic analysis reflects the variation law of the river section flow velocity in the transverse direction, the vertical variation characteristic analysis reflects the variation law of the river section flow velocity in the vertical direction, and the spatial distribution variation characteristic analysis reflects the variation trend of the river section flow field;

[0042] In step 3, after obtaining the average flow velocity [AX] at different characteristic lines in the transverse X direction and the average flow velocity [AY] at different characteristic lines in the vertical Y direction, the spatial distribution variation characteristic of the river section flow field is analyzed. The specific method is: based on the river section flow field interpolation data [VI], the Theil-Sen Median trend analysis is used to explore the variation trend of the spatial distribution time series of the river section flow field, the M-K test is coupled to analyze the significance of the variation trend, and the river section flow field spatial distribution variation characteristic C is generated.

[0043] ;

[0044] β is the variation slope of the spatial distribution time series of the river section flow field; and Z is the significance of the variation trend of the spatial distribution time series of the river section flow field.

[0045] ;

[0046] VI m is the interpolation data of the mth time series of the river section flow field; VI n is the interpolation data of the nth time series of the river section flow field; n < m;

[0047] When β > 0, it indicates that the spatial distribution time series of the river section flow field presents an accelerating trend; and when β < 0, it indicates that the spatial distribution time series of the river section flow field presents a decelerating trend.

[0048] ;

[0049] ;

[0050] ;

[0051] ;

[0052] S is the test statistic; VAR is the variance function; SIGN is the sign function; VI h is the interpolated data of the flow field of the river section of the hth time series; VI i is the interpolated data of the flow field of the river section of the ith time series;

[0053] If |Z| > Z 1−α / 2 , it indicates that there is a significant trend in the spatial distribution of the flow field of the river section over time.

[0054] Taking the confidence level α as 0.05, Z 1−α / 2 = 1.96, the spatial distribution change characteristic C of the flow field of the river section is divided into four types: significant slowing down (β < 0, |Z| > 1.96), slight slowing down (β < 0, |Z| < 1.96), slight speeding up (β > 0, |Z| < 1.96), and significant speeding up (β > 0, |Z| > 1.96).

[0055] Further, the detailed method of the step 4 flow field dynamic display is:

[0056] a) Based on the lateral variation characteristic [AX] of the flow field of the river section, the lateral variation characteristic line graph time series [F AX ] of the flow field of the river section is generated by the river section flow field time display module;

[0057] ;

[0058] LINE is a line graph function; F AX,m is the lateral variation characteristic line graph of the flow field of the river section of the mth time series;

[0059] b) Based on the vertical variation characteristic [AY] of the flow field of the river section, the vertical variation characteristic line graph time series [F AY ] of the flow field of the river section is generated by the river section flow field time display module;

[0060] ;

[0061] LINE is a line graph function; F AY,m is the vertical variation characteristic line graph of the flow field of the river section of the mth time series;

[0062] c) generating the spatial distribution heat map time series [H of the river cross-section flow field based on the river cross-section flow field interpolation data [VI] through the river cross-section flow field spatial display module VI ] using a gradient color band from red to green to present the river flow rate from slow to fast;

[0063] ;

[0064] HEAT is the heat map function; H VI,m is the spatial distribution heat map of the mth time series river cross-section flow field;

[0065] d) generating the time display results [RT] based on the horizontal variation characteristic line graph time series [F AX ], the vertical variation characteristic line graph time series [F AY ] and the spatial distribution heat map time series [H VI ] of the river cross-section flow field, and dynamically presenting the time variation characteristics of the river cross-section flow field through the time axis;

[0066] ;

[0067] ;

[0068] RT m is the time display results of the mth time series river cross-section flow field.

[0069] e) generating the spatial distribution variation characteristic heat map H C of the river cross-section flow field based on the spatial distribution variation characteristic C of the river cross-section flow field through the river cross-section flow field spatial display module, and using red, orange, blue and green to respectively present the significant slowing down, slight slowing down, slight speeding up and significant speeding up of the river flow rate;

[0070] ;

[0071] Finally, based on the spatial distribution variation characteristic heat map H C of the river cross-section flow field, the spatial display results RS are generated; The application also discloses a river cross section flow field data collection and analysis method system, which comprises a flow collection system, a data processing system, a flow field analysis system and a flow field display system; the river cross section flow field data collected by the flow collection system is processed by the data processing system, the flow field analysis system and the flow field display system in sequence and finally displayed dynamically; wherein the flow collection system comprises a control module, a mechanical module and a collection module, the mechanical module is in the shape of "H", two supports are vertically arranged on the two banks of the river, the top ends of the two supports are horizontally arranged with guide rails, the collection module is hung on the river surface through a cable and a pulley and is moved horizontally along the river surface through the guide rails, the depth in the vertical direction is adjusted through the cable, and thus the collection module can collect the initial data [VO] of the river cross section flow field in a preset collection period; the control module controls the motion instruction of the mechanical module.

[0072] Further, the control module is implanted with a default instruction set [Rule] through automatic programming, the instruction set [Rule] comprises a start instruction ON, a pause instruction OF and a reset instruction RE, each instruction is provided with a corresponding instruction code, the instruction code is composed of an identification code VV, a function code FF, a parameter code PP and an end code EE, and the parameter code PP is of indefinite length; the expression is as follows:

[0073] Rule = [Rule1, Rule2, …, Rule t ];

[0074] Rule t = [VV t ,FF t ,PP t ,EE t ];

[0075] The start instruction ON realizes the start function of mechanical movement, and the parameter code comprises the vector path information of mechanical movement;

[0076] PP = DTxxDXxxDYxx…DXxxDYxx…

[0077] In the formula, DX is the movement step length identifier in the X direction, the data segment length is 2 characters; DY is the movement step length identifier in the Y direction, the data segment length is 2 characters; DT is the single step movement time interval identifier, the data segment length is 2 characters, and the unit is second; wherein the DX and DY identifiers can appear multiple times, and the data segment is encoded and converted in hexadecimal;

[0078] The pause instruction OF realizes the pause function of mechanical movement, the parameter code PP = DTxxxxxx, and xxxxxx is the total length of the pause, wherein x is a data code symbol, x = 1~F, and the unit is second;

[0079] The reset command RE returns the probe to its default position. The parameter code PP = DTxxDXDY can be configured to follow the order of X and Y directions.

[0080] Beneficial effects: Compared with the prior art, the present invention has the following advantages:

[0081] (1) This invention provides an integrated system for collecting, processing, analyzing and displaying the flow field of a river cross section, thereby realizing the automated execution of the entire process of flow field data collection, processing, analysis and display.

[0082] (2) This invention comprehensively collects the overall characteristic data of the flow field of the river section, and the analysis results are more scientific.

[0083] (3) This invention automatically analyzes the lateral, vertical and spatial distribution variation characteristics of the flow field of the river section, calculates the time series of the spatiotemporal variation characteristics of the flow field of the river section, and dynamically displays the results data in combination with the time axis, so as to quickly, comprehensively and visually present the spatiotemporal variation characteristics of the flow field of the river section. Attached Figure Description

[0084] Figure 1 This is an overall flowchart of the present invention;

[0085] Figure 2 This is a schematic diagram of the flow field acquisition system in the embodiment;

[0086] Figure 3 This is a schematic diagram of the mechanical module structure of the flow field acquisition system in the embodiment;

[0087] Figure 4 This is a schematic diagram of the acquisition module structure of the flow field acquisition system in the embodiment;

[0088] Figure 5 This is a schematic diagram of the lateral variation characteristics in the embodiment;

[0089] Figure 6 This is a schematic diagram of the vertical variation characteristics in the embodiment. Detailed Implementation

[0090] The technical solution of the present invention will be described in detail below, but the scope of protection of the present invention is not limited to the embodiments described.

[0091] like Figure 1 As shown, the present invention provides a method for acquiring and analyzing cross-sectional flow field data of a river, comprising the following steps:

[0092] Step 1: Collect initial flow field data [VO] of the river cross-section within a preset time period using a flow field acquisition system.

[0093] ;

[0094] ;

[0095] ;

[0096] ;

[0097] m is the total number of collection time series in a preset collection period; VO m is the initial data of the river section flow field of the mth time series; VO (r, s) is the initial data of the river flow rate at the monitoring point of coordinates (r, s); r is the total number of monitoring points in the X direction of the river; s is the total number of monitoring points in the Y direction of the river; W is the width of the river section; and D is the depth of the river section;

[0098] Step 2: processing the initial data [VO] of the river section flow field obtained in step 1 by the data processing system to obtain, in sequence, effective data [VE] of the river section flow field, interpolation data [VI] of the river section flow field, and storage data [VS] of the river section flow field;

[0099] ;

[0100] (VO m , VE m , VI m ) is the storage data of the river section flow field of the mth time series;

[0101] ;

[0102] ;

[0103] EFF is a 3σ effective analysis function; VE m is the effective data of the river section flow field of the mth time series; VE (r, s) is the effective data of the river flow rate at the monitoring point of coordinates (r, s);

[0104]

[0105] ;

[0106] ;

[0107] ;

[0108] KRI is a Kriging interpolation function; VI m is the interpolation data of the river section flow field of the mth time series; VI (j, k)The river flow rate interpolation data at the monitoring point of coordinates (j, k); j is the total number of interpolation points in the X direction of the river; k is the total number of interpolation points in the Y direction of the river; ΔX is the interpolation resolution in the X direction of the river; and ΔY is the interpolation resolution in the Y direction of the river;

[0109] Step 3, analyzing the effective data [VE] of the river section flow field by the flow field analysis system to obtain the average flow rate [AX] at different characteristic lines in the transverse X direction and the average flow rate [AY] at different characteristic lines in the vertical Y direction, and analyzing the interpolation data [VI] of the river section flow field to obtain the spatial distribution change characteristic C of the river section flow field;

[0110] ;

[0111] ;

[0112] ;

[0113] AX m is the transverse change characteristic of the river section flow field of the mth time sequence; AX X=r is the average flow rate of the rth characteristic line in the transverse direction;

[0114] ;

[0115] ;

[0116] ;

[0117] AY m is the vertical change characteristic of the river section flow field of the mth time sequence; AY Y=s is the average flow rate of the st characteristic line in the vertical direction;

[0118] ;

[0119] β is the change slope of the spatial distribution time sequence of the river section flow field; and Z is the significance of the change trend of the spatial distribution time sequence of the river section flow field;

[0120] Step 4, dynamically displaying the data obtained in step 3, and the specific content includes: dynamically presenting the transverse change characteristic broken line graph time sequence [F AX ] of the river section flow field in a preset collection period, the vertical change characteristic broken line graph time sequence [F AY ], the spatial distribution heat map time sequence [H VI ], and the spatial distribution change characteristic heat map H C of the river section flow field.

[0121] Further, the flow field acquisition system in step 1 is provided with an acquisition module, which acquires the transverse flow velocity and longitudinal flow velocity of each coordinate point of the river section through a sensing probe in real time to generate the initial data [VO] of the river section flow field in a preset acquisition period.

[0122] Further, the data processing system in step 2 includes an effective processing module, an interpolation processing module and a data storage module, the effective processing module performs effectiveness judgment on the initial data [VO] of the river section flow field to obtain effective data [VE] of the river section flow field, the interpolation processing module performs spatial interpolation on the obtained effective data [VE] of the river section flow field to generate interpolation data [VI] of the river section flow field, and the data storage module stores the above data; the specific method is as follows:

[0123] Step 2.1, based on the initial data [VO] of the river section flow field, the effective processing module performs effectiveness judgment on the data by using 3σ effective analysis method to eliminate invalid values and device null values, and generates effective data [VE] of the river section flow field;

[0124] Step 2.2, based on the effective data [VE] of the river section flow field, the interpolation processing module performs spatial interpolation processing on the data by using Kriging interpolation method to generate interpolation data [VI] of the river section flow field;

[0125] Step 2.3, the data storage module stores the initial data [VO], effective data [VE] and interpolation data [VI] in JSON format; and generates storage data [VS] of the river section flow field;

[0126] ;

[0127] (VO m , VE m , VI m ) is the storage data of the river section flow field of the mth time sequence.

[0128] Further, the flow field analysis system includes river section flow field transverse variation characteristic analysis, vertical variation characteristic analysis and spatial distribution variation characteristic analysis. The transverse variation characteristic analysis reflects the variation law of the river section flow velocity in the transverse direction, the vertical variation characteristic analysis reflects the variation law of the river section flow velocity in the vertical direction, and the spatial distribution variation characteristic analysis reflects the variation trend of the river section flow field.

[0129] In step 3, after obtaining the average flow rate [AX] at different characteristic lines in the transverse X direction and the average flow rate [AY] at different characteristic lines in the vertical Y direction, the spatial distribution variation characteristics of the river section flow field are analyzed. The specific method is as follows: based on the interpolation data [VI] of the river section flow field, the Theil-Sen Median trend analysis is used to explore the variation trend of the spatial distribution time series of the river section flow field, the M-K test is coupled to analyze the significance of the variation trend, and the variation characteristics C of the spatial distribution of the river section flow field are generated.

[0130] ;

[0131] β is the change slope of the spatial distribution time series of the river section flow field; Z is the significance of the variation trend of the spatial distribution time series of the river section flow field.

[0132] ;

[0133] VI m is the interpolation data of the river section flow field of the mth time series; VI n is the interpolation data of the river section flow field of the nth time series; n < m. When β > 0, it indicates that the spatial distribution time series of the river section flow field presents an accelerating trend; when β < 0, it indicates that the spatial distribution time series of the river section flow field presents a slowing trend.

[0134] ;

[0135] ;

[0136] ;

[0137] ;

[0138] S is the test statistic; VAR is the variance function; SIGN is the sign function; VI h is the interpolation data of the river section flow field of the hth time series; VI i is the interpolation data of the river section flow field of the ith time series.

[0139] If |Z| > Z 1−α / 2 , it indicates that there is a significant variation trend in the spatial distribution time series of the river section flow field.

[0140] Take the confidence level α as 0.05, Z 1−α / 2= 1.96, the spatial distribution variation characteristics C of the river cross-section flow field is divided into four types: significant slowing down (β < 0, |Z| > 1.96), slight slowing down (β < 0, |Z| < 1.96), slight speeding up (β > 0, |Z| < 1.96) and significant speeding up (β > 0, |Z| > 1.96).

[0141] Further, the detailed method of the step 4 flow field dynamic display is:

[0142] a) Based on the lateral variation characteristic [AX] of the river cross-section flow field, the lateral variation characteristic line graph time sequence [F AX ] of the river cross-section flow field is generated by the river cross-section flow field time display module;

[0143] ;

[0144] LINE is a line graph function; F AX,m is the lateral variation characteristic line graph of the mth time sequence river cross-section flow field;

[0145] b) Based on the vertical variation characteristic [AY] of the river cross-section flow field, the vertical variation characteristic line graph time sequence [F AY ] of the river cross-section flow field is generated by the river cross-section flow field time display module;

[0146] ;

[0147] LINE is a line graph function; F AY,m is the vertical variation characteristic line graph of the mth time sequence river cross-section flow field;

[0148] c) Based on the interpolation data [VI] of the river cross-section flow field, the spatial distribution heat map time sequence [H VI ] of the river cross-section flow field is generated by the river cross-section flow field time display module, and a gradual color band from red to green is used to present the river flow speed from slow to fast;

[0149] ;

[0150] HEAT is a heat map function; H VI,m is the spatial distribution heat map of the mth time sequence river cross-section flow field;

[0151] d) Based on the lateral variation characteristic line graph time sequence [F AX ], the vertical variation characteristic line graph time sequence [F AY ] and the spatial distribution heat map time sequence [H VI, generate time display results [RT], and dynamically present the time variation characteristics of the river section flow field through the time axis;

[0152] ;

[0153] ;

[0154] RT m is the time display result of the mth time series river section flow field.

[0155] e) Based on the spatial distribution variation characteristics C of the river section flow field, the spatial distribution variation characteristics heat map H of the river section flow field is generated through the river section flow field spatial display module C , red, orange, blue, and green are used to represent significant slowing down, slight slowing down, slight acceleration, and significant acceleration of river flow velocity, respectively;

[0156] ;

[0157] Finally, based on the spatial distribution variation characteristics heat map H of the river section flow field C , the spatial display result RS is generated. .

[0158] The system of the above-mentioned river section flow field data collection and analysis method includes a flow collection system, a data processing system, a flow field analysis system, and a flow field display system. The river section flow field data collected by the flow field collection system is processed by the data processing system, the flow field analysis system, and the flow field display system in turn, and finally dynamically displayed.

[0159] As shown in Figures 2 to 4 , the flow field collection system includes a control module, a mechanical module, and a collection module. The mechanical module is in the shape of "H". Two H-shaped supports 1 are vertically erected on both sides of the river. The top ends of the two supports 1 are horizontally laid with guide rails 3. The collection module is suspended on the river surface through cables 4 and pulleys 2, and is moved horizontally along the river surface through the guide rails 3. The depth in the vertical direction is adjusted through the cables 4, so as to realize the collection of the initial data [VO] of the river section flow field in the preset collection period. The control module controls the motion instructions of the mechanical module.

[0160] The control module implants a default instruction set [Rule] through automatic programming. The instruction set [Rule] includes start instruction ON, pause instruction OF, and reset instruction RE. Each instruction is provided with a corresponding instruction code. The instruction code is composed of identification code VV, function code FF, parameter code PP, and end code EE. The parameter code PP is of indefinite length. The expression is as follows:

[0161] Rule = [Rule1, Rule2, …, Rulet

[0162] Rule t = [VV t ,FF t ,PP t ,EE t ];

[0163] The start instruction ON realizes the mechanical motion start function, and the parameter code includes the vector path information of the mechanical motion.

[0164] PP = DTxxDXxxDYxx…DXxxDYxx…

[0165] In the above formula, DX is a movement step length identifier in the X direction, and the data segment length thereof is 2 characters; DY is a movement step length identifier in the Y direction, and the data segment length thereof is 2 characters; and DT is a single step moving time interval identifier, and the data segment length thereof is 2 characters, with the unit being seconds; wherein the DX and DY identifiers can appear multiple times, and the data segment is converted according to the hexadecimal encoding;

[0166] The pause instruction OF realizes the mechanical motion pause function, and the parameter code PP = DTxxxxxx, wherein xxxxxx is the total time length of the pause, x is a data code symbol, x = 1 ~ F, and the unit is seconds;

[0167] The reset instruction RE realizes the probe returning to the default position, and the parameter code PP = DTxxDXDY, which can be configured in the order of the X direction and the Y direction.

[0168] Table 1 Control instruction description table

[0169]

[0170] Embodiment

[0171] To verify the technical effect of the present application, the embodiment is applied to a typical plain river channel in a certain area in the south, the width of the river channel is 20 m, and the depth is 5 m. The specific collection and analysis method is as follows:

[0172] Step S1 collects data through the flow field collection system;

[0173] The total number of collection time sequences in the preset collection period is m = 10, the movement step length in the X direction is DX = 2, the movement step length in the Y direction is DY = 1, the total number of monitoring points in the X direction of the river channel is r = 10, the total number of monitoring points in the Y direction of the river channel is s = 5, and the initial data [VO] of the river channel section flow field is collected through the flow field collection system.

[0174] ​​

[0175] Take the first acquisition time sequence as an example (m = 1):

[0176] ;

[0177] S2 processes the collected data through the data processing system;

[0178] (1) Effective processing module

[0179] Based on the initial data [VO] of the river section flow field, 3σ effective analysis is used to judge the effectiveness of the data, and invalid values and device null values are removed to generate effective data [VE] of the river section flow field.

[0180]

[0181] Take the first acquisition time sequence as an example (m = 1):

[0182] ;

[0183] (2) Interpolation processing module

[0184] Based on the effective data [VE] of the river section flow field, the Kriging interpolation method is used for spatial interpolation processing of the data to generate interpolation data [VI] of the river section flow field. The interpolation resolution of the river X direction is ΔX = 0.5, the interpolation resolution of the river Y direction is ΔY = 0.5, the total number of interpolation points of the river X direction is j = 40, and the total number of interpolation points of the river Y direction is k = 10.

[0185] ;

[0186] Take the first acquisition time sequence as an example (m = 1):

[0187] ;

[0188] (3) Data storage module

[0189] The data storage module stores the original data, effective data and interpolation data in JSON format to generate storage data [VS] of the river section flow field.

[0190] ;

[0191] Step S3 analyzes the processed data through the flow field analysis system;

[0192] (1) River section flow field transverse variation characteristics

[0193] Based on the effective data [VE] of the river section flow field, the average flow velocity [AX] at different characteristic lines in the transverse X direction is calculated.

[0194] ;

[0195] Take the first acquisition time sequence as an example (m = 1):

[0196] ;

[0197] (2) Vertical variation characteristics of river section flow field

[0198] Based on the effective data [VE] of the river section flow field, the average flow velocity [AY] at different characteristic lines in the vertical Y direction is calculated.

[0199] ;

[0200] Take the first acquisition time sequence as an example (m = 1):

[0201] ;

[0202] (3) Spatial distribution variation characteristics of river section flow field

[0203] Based on the interpolation data [VI] of the river section flow field, Theil-Sen Median trend analysis is used to explore the variation trend of the spatial distribution time sequence of the river section flow field, coupled with M-K test analysis of the significance of the variation trend, to generate the spatial distribution variation characteristics C of the river section flow field.

[0204] ;

[0205] Step S4 dynamically displays the obtained analyzed river section data through a flow field display system;

[0206] (1) River section flow field time display module

[0207] Based on the lateral variation characteristics [AX] of the river section flow field, the lateral variation characteristic line graph time sequence [F AX ] of the river section flow field is generated.

[0208] ;

[0209] Take the first acquisition time sequence as an example (m = 1), the lateral variation characteristic line graph F AX,1 , as shown in Figure 5 ;

[0210] Based on the vertical variation characteristics [AY] of the river section flow field, the vertical variation characteristic line graph time sequence [F AY ] of the river section flow field is generated.

[0211] ;

[0212] Taking the first collected time series as an example (m = 1), the vertical variation characteristic broken line graph F AY,1 ,like Figure 6 As shown;

[0213] Based on the interpolated flow field data of the river cross section [VI], a time series of spatial distribution heatmaps of the flow field of the river cross section [H] is generated. VI The gradient color band from red to green is used to represent the river's flow speed from slow to fast.

[0214] ;

[0215] Time series of broken line graphs based on the lateral variation characteristics of the flow field in the river cross section [F] AX Vertical variation characteristic broken line graph time series [F] AY ] and spatial distribution heatmap time series [H VI The [RT] method generates time-based display results, dynamically presenting the temporal variation characteristics of the river cross-section flow field through a timeline.

[0216] ;

[0217] Taking the first collected time series as an example (m = 1):

[0218] ;

[0219] (2) Spatial display module of river cross-section flow field

[0220] Based on the spatial distribution variation characteristics C of the river cross-section flow field, a heat map H of the spatial distribution variation characteristics of the river cross-section flow field is generated. C The colors red, orange, blue, and green respectively represent a significant decrease in river flow velocity, a slight decrease, a slight increase, and a significant increase.

[0221] ;

[0222] Heatmap H based on the spatial distribution variation characteristics of the river cross-section flow field C Generate spatial display results RS, .

Claims

1. A method for collecting and analyzing river cross-section flow field data, characterized in that, The method comprises the following steps: Step 1, collecting initial data [VO] of the river section flow field in a preset period of time through a flow field collection system, ； ； ； ； m is the total number of collection time series in a preset collection period; VO m is the initial data of the river section flow field of the mth time series; VO (r, s) is the initial data of the river flow velocity at the monitoring point of coordinates (r, s); r is the total number of monitoring points in the X direction of the river; s is the total number of monitoring points in the Y direction of the river; W is the width of the river section; and D is the depth of the river section. Step 2, processing the initial data [VO] of the river section flow field obtained in step 1 through a data processing system to obtain effective data [VE] of the river section flow field, interpolation data [VI] of the river section flow field and storage data [VS] of the river section flow field in sequence; ； (VO m , VE m , VI m ) is the data of the river section flow field of the mth time sequence; ； ； EFF is the 3σ effective analysis function; VE m is the effective data of the river section flow field of the mth time series; VE (r, s) is the effective data of the river flow velocity at the monitoring point with coordinates (r, s). ； ； ； ； KRI is a Kriging interpolation function; VI m is the interpolation data of the river section flow field of the mth time series; VI (j, k) is the interpolation data of the river flow velocity at the coordinate (j, k) monitoring point; j is the total number of interpolation points in the X direction of the river; k is the total number of interpolation points in the Y direction of the river; ΔX is the interpolation resolution in the X direction of the river; and ΔY is the interpolation resolution in the Y direction of the river. Step 3, analyzing the effective data [VE] of the river section flow field through a flow field analysis system to obtain average flow velocities [AX] at different characteristic lines in the horizontal X direction and average flow velocities [AY] at different characteristic lines in the vertical Y direction, and analyzing the interpolation data [VI] of the river section flow field to obtain a spatial distribution change characteristic C of the river section flow field; ； ； ； AX m is the lateral variation characteristic of the flow field of the river section for the mth time series; AX X=r is the average flow velocity of the rth characteristic line in the lateral direction. ； ； ； AY m is the vertical variation characteristic of the flow field of the river section for the mth time series;AY Y=s is the average flow velocity of the vertical sth characteristic line. ； β is a change slope of the spatial distribution time sequence of the river section flow field; and Z is a significance of the change trend of the spatial distribution time sequence of the river section flow field; In step 3, after the average flow velocities [AX] at different characteristic lines in the horizontal X direction and the average flow velocities [AY] at different characteristic lines in the vertical Y direction are obtained, the spatial distribution change characteristic of the river section flow field is analyzed, and the specific method is as follows: based on the interpolation data [VI] of the river section flow field, a Theil-Sen Median trend analysis is used to explore the change trend of the spatial distribution time sequence of the river section flow field, and a M-K test is coupled to analyze the significance of the change trend, so as to generate the spatial distribution change characteristic C of the river section flow field; ； β is a change slope of the spatial distribution time sequence of the river section flow field; and Z is a significance of the change trend of the spatial distribution time sequence of the river section flow field; ； VI m is the interpolated data of the flow field of the river section for the mth time series;VI n is the interpolated data of the flow field of the river section for the nth time series; n < m; When β > 0, it indicates that the spatial distribution time sequence of the river section flow field presents an accelerating trend; When β < 0, it indicates that the spatial distribution time sequence of the river section flow field presents a slowing trend; ； ； ； ； S is the test statistic; VAR is the variance function; SIGN is the sign function; VI h is the interpolated data of the flow field of the river section for the hth time series; VI i is the interpolated data of the flow field of the river section for the ith time series; If |Z| > Z 1−α / 2 , it indicates that there is a significant trend in the spatial distribution of the time series of the river cross-section flow field. Step 4, the data obtained in step 3 is dynamically displayed, and the specific content includes: dynamically presenting the horizontal variation characteristic broken line graph time sequence [F AX ] of the river section flow field in the preset collection period, the vertical variation characteristic broken line graph time sequence [F AY ], the spatial distribution heat map time sequence [H VI ], and the river section flow field spatial distribution variation characteristic heat map H C .

2. The river cross-section flow field data acquisition and analysis method according to claim 1, characterized in that, The flow field collection system in step 1 is provided with a collection module, the collection module collects the horizontal flow velocity and the vertical flow velocity of each coordinate point of the river section through a sensing probe in real time, and generates the initial data [VO] of the river section flow field in a preset collection period.

3. The river cross-section flow field data acquisition and analysis method according to claim 1, characterized in that, The data processing system in step 2 comprises an effective processing module, an interpolation processing module and a data storage module, the effective processing module judges the effectiveness of the initial data [VO] of the river section flow field to obtain the effective data [VE] of the river section flow field, the interpolation processing module performs spatial interpolation on the obtained effective data [VE] of the river section flow field to generate the interpolation data [VI] of the river section flow field, and the data storage module stores the above data; and the specific method is as follows: Step 2.1, based on the initial data [VO] of the river section flow field, the effective processing module uses a 3σ effective analysis method to judge the effectiveness of the data, removes invalid values and device null values, and generates the effective data [VE] of the river section flow field; Step 2.2, based on the effective data [VE] of the river section flow field, the interpolation processing module uses a Kriging interpolation method to perform spatial interpolation processing on the data, and generates the interpolation data [VI] of the river section flow field; Step 2.3, the data storage module stores the initial data [VO], valid data [VE] and interpolation data [VI] in JSON format; and generates river section flow field storage data [VS]; ； (VO m , VE m , VI m ) are data storage for the flow field of the river section of the mth time sequence.

4. The river cross-section flow field data acquisition and analysis method according to claim 1, characterized in that The detailed method of the step 4 flow field dynamic display is: a) generating a time series of the lateral variation characteristic polyline graph of the river section flow field [F based on the lateral variation characteristic [AX] of the river section flow field through the river section flow field time display module AX ] ； LINE is a broken line graph function; F AX,m is the broken line graph of the transverse variation characteristics of the mth time series river cross-section flow field; b) generating a vertical variation characteristic broken line graph time sequence [F of the river section flow field based on the vertical variation characteristic [AY] of the river section flow field, by a river section flow field time display module AY ] ； LINE is a broken line graph function; F AY,m is the vertical variation characteristic broken line graph of the mth time series river cross-section flow field; c) generating a spatial distribution thermograph time series [H] of the river section flow field by a river section flow field time display module based on the river section flow field interpolation data [VI] VI ] ； HEAT is the heat map function; H VI,m is the heat map of the spatial distribution of the flow field of the mth time series river cross-section; d) a time series of the broken line graph of the lateral variation characteristics of the river cross-section flow field [F AX ], a time series of the broken line graph of the vertical variation characteristics of the river cross-section flow field [F AY ], and a time series of the spatial distribution thermodynamic map of the river cross-section flow field [H VI ], to generate a time display result [RT] and dynamically present the time variation characteristics of the river cross-section flow field through a time axis. ； ； RT m is the time display result of the mth time series river cross-section flow field; e) generating a heat map H of the spatial distribution variation characteristics of the flow field of the river section based on the spatial distribution variation characteristics C of the flow field of the river section, through a spatial distribution variation characteristics heat map generation module of the flow field of the river section C ; ； Finally, based on the spatial distribution of the river cross-section flow field change characteristics heat map H C , generate spatial display results RS; .

5. A system for implementing the method for collecting and analyzing data of flow field of river cross section according to any one of claims 1 to 4, characterized in that, The system comprises a flow acquisition system, a data processing system, a flow field analysis system and a flow field display system; The river section flow field data acquired by the flow acquisition system is processed by the data processing system, the flow field analysis system and the flow field display system in sequence, and is finally dynamically displayed; The flow acquisition system comprises a control module, a mechanical module and an acquisition module, the mechanical module is in the shape of H as a whole, two supports are vertically arranged on the two banks of the river, and a guide rail is horizontally arranged at the top of the two supports, the acquisition module is suspended on the river surface through a cable and a pulley, and is moved along the guide rail in the horizontal direction of the river surface, the depth in the vertical direction is adjusted through the cable, so that the acquisition module can acquire the initial data [VO] of the river section flow field in a preset acquisition period; and the control module controls the movement instruction of the mechanical module.

6. The system for river cross-section flow field data acquisition and analysis method according to claim 5, wherein The control module is programmed to implant a default instruction set [Rule], the instruction set [Rule] comprises a start instruction ON, a pause instruction OF and a reset instruction RE, each instruction is provided with a corresponding instruction code, and the instruction code is composed of an identification code VV, a function code FF, a parameter code PP and an end code EE, wherein the parameter code PP is of indefinite length; and the expression is as follows: Rule = [Rule1, Rule2, …, Rule t ] ; Rule t = [VV t ,FF t ,PP t ,EE t ] The start instruction ON realizes the start function of the mechanical movement, and the parameter code comprises the vector path information of the mechanical movement; PP = DTxxDXxxDYxx…DXxxDYxx… In the formula, DX is an identifier of the movement step length in the X direction, DY is an identifier of the movement step length in the Y direction, and DT is an identifier of the movement time interval of a single step; wherein the identifiers DX and DY can appear multiple times; The pause instruction OF realizes the pause function of the mechanical movement, and the parameter code PP = DTxxxxxx, wherein xxxxxx is the total length of the pause, x is a data code symbol, and x = 1~F; The reset instruction RE realizes the function of returning the probe to the default position, and the parameter code PP = DTxxDXDY, which configures the sequence of the movement in the X direction and the Y direction.

Images