Water level estimation program, water level estimation system, and water level estimation method

Abstract

To provide a water level condition estimation system and an estimation method capable of recognizing a wide range of water level conditions in real time.SOLUTION: A water level condition estimation method according to one aspect of the present invention has the steps of: acquiring the estimation conditions; recording image data; recording water level data based on the image data; creating water level condition estimation data by an extended Kalman filter FEM based on the estimation conditions and the water level data; and displaying the water level condition estimation data. The water level condition estimation system according to another aspect of the present invention includes a detection unit, an image data recording unit, a water level data recording unit, a water level status estimation data acquisition unit for acquiring water level condition estimation data by the extended Kalman filter FEM, an estimation condition acquisition unit, and a display unit for executing the method.SELECTED DRAWING: Figure 4

Description

[Technical field]

[0001] The present invention relates to a water level estimation program, a water level estimation system, and a water level estimation method. [Background technology]

[0002] For example, Patent Document 1 discloses a state prediction device capable of accurately predicting water level conditions of known tsunamis and the like in real time, which includes a prediction unit that predicts a state vector at the next time point using two-dimensional shallow water equations that represent the propagation of a tsunami, for a state vector composed of the tsunami flow rate and water level at multiple points set two-dimensionally in an area including the radar coverage area, a smoothing unit that uses a Kalman filter to smooth the state vector predicted by the prediction unit collectively within the coverage area using a Kalman gain, the state vector predicted by the prediction unit, and an observation vector composed of sea surface flow velocity observation values ​​in multiple cells spanning multiple range directions and multiple beam directions within the coverage area, and a setting unit that sets an initial value used to predict the state vector in the prediction unit.

[0003] However, not only is it important to predict the water level conditions at the time of tsunami occurrence, but it is also important from the perspective of disaster prevention and mitigation to always be aware of the water level conditions at each point in a specific range, i.e., the water level situation. As a technology for estimating the water level situation over a wide area without actually observing the water level conditions at each point, the inventors of the present invention have developed a tidal current prediction technology for rivers and bays that uses a shallow water area flow condition estimation analysis technology based on a Kalman filter finite element method (FEM) that takes into account the movement of the waterfront, as described in Non-Patent Document 1. [Prior art documents] [Patent documents]

[0004] [Patent Document 1] Patent No. 6641532 [Non-patent literature]

[0005] [Non-Patent Document 1] Proceedings of the 32nd Annual Meeting of the Japanese Society for Computational Mechanics (2019) Vol. 32, No. 013 Summary of the Invention [Problem to be solved by the invention]

[0006] When grasping water level conditions over a wide area, it is desirable to use actual observation values ​​from each point, but from the perspective of system scale, there is a need for flow regime estimation analysis technology that can match all actual observation values ​​as closely as possible in real time while using as few observation values ​​as possible.

[0007] While the inventors were studying the development of the above-mentioned analysis technology, they performed a flow regime estimation analysis using water level data obtained using conventional observation technology compared to data obtained by a water level gauge installed in an actual bay, but discovered a problem in that the resulting water levels differed from the actually measured water levels, even when a numerical analysis method that included the effect of high numerical viscosity was used.

[0008] In view of the above problems, the present invention aims to provide a water level estimation program, a water level estimation system, and a water level estimation method that can grasp water level conditions over a wide area in real time with greater accuracy. [Means for solving the problem]

[0009] The inventors of the present invention conducted intensive research into the above-mentioned problems and discovered that by performing estimation analysis using an extended Kalman filter (FEM) that takes into account the movement of the waterfront using observed water levels, it is possible to grasp the water level conditions over a wide area in real time, which led to the completion of the present invention.

[0010] That is, a water level situation estimation program according to one aspect of the present invention for solving the above problem is provided on a computer. Observation points located upstream of the water flow The method includes the steps of recording image data, recording water level data based on the image data, and creating water level status estimation data using the extended Kalman filter finite element method based on the estimation conditions and the water level data. In addition, the extended Kalman filter finite element method uses equation (1) as the system equation, and calculates the state transition matrix after linear approximation at time step n [Fn ] and the prediction error covariance matrix [P (+) n ], the estimation error covariance matrix [P (-) n+1 ] and the estimation error covariance matrix [P (-) n+1 ], the Kalman gain matrix [K 1 n+1 ] and the estimation error covariance matrix [P (-) n+1 ] and the Kalman gain matrix [K 1 n+1 ], the prediction error covariance matrix [P (+) n+1 Calculating Further, the method includes a step of providing minute water depths to adjacent junctions in elements including nodes of waterfront boundaries consisting of wet and dry regions; Furthermore, the state transition matrix [A] and the estimated value vector {φ^ (-) n }, the estimated vector {φ^ (-) n+1 } and calculate the estimated vector {φ^ (-) n+1 } and water level data and the Kalman gain matrix [K 1 n+1 ], the optimal estimate vector {φ^ (+) n+1 Preferably, the method further comprises the step of calculating .times. ...

[0011] In this aspect, the step of recording the water level data includes binarizing the image data and performing noise reduction processing. Remove It is preferable to carry out at least one of the above processes.

[0012] In addition, in this respect, the water level situation estimation data is The micro depth was given to the elements including the nodes of the water boundary consisting of the wet area and the dry area. It is preferable that the extended Kalman filter is created by the finite element method.

[0013] In addition, in this aspect, it is preferable that the water level condition estimation data includes water area estimation data and water depth estimation data.

[0014] In addition, in this aspect, the step of recording the image data includes: The state before flooding or inundation occurs It is preferable that the step of recording the image data and the image data at the time of estimation and recording the water level data includes recording the water level data based on the image data at the time of estimation.

[0015] In addition, in this aspect, it is preferable to have a step of displaying the water level situation estimation data.

[0016] In addition, a water level situation estimation system according to another aspect of the present invention includes: Upstream of the flow of a body of water It comprises a plurality of detection units installed at different observation points, an image data recording unit which records image data based on the output from the detection units, a water level data recording unit which records water level data based on the image data, and a water level condition estimation data acquisition unit which acquires water level condition estimation data using the extended Kalman filter finite element method based on the estimation conditions and the water level data. In addition, in the water level estimation data acquisition section, the system equation in the extended Kalman filter finite element method is used to obtain the state transition matrix after linear approximation at time step n [F n ] and the prediction error covariance matrix [P (+) n ], the estimation error covariance matrix [P (-) n+1 ] and calculate the estimated error covariance matrix [P (-) n+1 ], the Kalman gain matrix [K 1 n+1 ] and calculate the estimated error covariance matrix [P (-) n+1 ] and the Kalman gain matrix [K 1 n+1 ], the prediction error covariance matrix [P (+) n+1 ], Furthermore, for elements that include nodes on the waterfront boundary consisting of wet and dry areas, a small water depth is given to the adjacent junctions, Furthermore, the state transition matrix [A] and the estimated value vector {φ^ (-) n }, the estimated vector {φ^ (-) n+1 } and calculate the estimated vector {φ^ (-) n+1 } and water level data and the Kalman gain matrix [K 1 n+1 ], the optimal estimate vector {φ^ (+) n+1 It is preferable to calculate

[0017] Furthermore, a water level estimation method according to another aspect of the present invention program teeth, Water level estimation method using a computer execution to make It is something. Effect of the Invention

[0018] As described above, the present invention can provide a water level estimation program, a water level estimation system, and a water level estimation method that can grasp water level conditions over a wide area in real time with higher accuracy. [Brief description of the drawings]

[0019] [Figure 1] FIG. 2 is a schematic diagram of functional blocks of the water level estimation system. [Diagram 2] This is an example of a connection including the hardware of a water level estimation system. [Diagram 3] 13 is a connection example including hardware relating to another example of a water level situation estimation system. [Figure 4] 13 is a flowchart showing a water level status estimation procedure. [Diagram 5] 13 is a flowchart showing a water level data recording procedure. [Figure 6] 13 is a flowchart showing a procedure for acquiring water level estimation data. [Figure 7] 1 is an example of image data for acquiring water level data. [Figure 8] 13 is an example of the water level data obtained. [Figure 9] This is an example of comparing water level estimation data (estimated values) with observed values. [Figure 10] This is an example of comparing estimated water level data with observed values. [Figure 11] 1 is an example showing distribution of flow velocity and water level fluctuation amount in an embodiment. [Figure 12] This is an example of a condition for estimating and analyzing water level conditions using water level data obtained from image analysis. [Figure 13] 13 is an example showing a comparison between a measured water level and an estimated water level in the embodiment. [Figure 14] 13 is an example showing a comparison between a measured water level and an estimated water level in the embodiment. BEST MODE FOR CARRYING OUT THEINVENTION

[0020] Hereinafter, the embodiments of the present invention will be described in detail with reference to the drawings. However, the present invention can be embodied in many different forms, and is not limited to the specific examples mentioned in the following embodiments and examples.

[0021] (Water level estimation system) FIG. 1 shows an outline of the functional blocks of a water level estimation system (hereinafter referred to as "this system") S according to this embodiment. As shown in this figure, this system S includes a detection unit 1, an image data recording unit 2 that records image data based on the output from the detection unit 1, a water level data recording unit 3 that records water level data based on the image data, a water level estimation data acquisition unit 4 that acquires water level estimation data by an extended Kalman filter finite element method (FEM), an estimation condition acquisition unit 5 that acquires analysis conditions such as mesh data and observation point data that are the initial values ​​for the calculation of the extended Kalman filter FEM, and estimation variables such as flow velocity and water level fluctuation, and a display unit 6 that displays the water level estimation data. Furthermore, this system S includes a recording medium 7 that records a water level estimation program that executes the extended Kalman filter FEM and the various data described above, and a volatile recording medium 8 that performs various calculations such as executing the extended Kalman filter FEM.

[0022] The effects of this system S will be described in detail later, but this system makes it possible to grasp the flood damage situation in real time by estimating the water level in a wide area other than the observation point when flooding occurs on land such as residential areas due to floods, typhoons, tsunamis, etc., and to grasp the damage situation around the river or bay area in real time by estimating the water level in a wide area other than the observation point when the water level in a river or bay rises due to a typhoon, etc. In particular, since water level data can be obtained based on image data obtained by the output of the detection unit, it is easy to obtain various information on the water level situation at the observation point, and by using this information, the system can estimate the water level situation with high accuracy.

[0023] Figure 2 shows an example of connections including the hardware of this system S. This system S can be roughly divided into, but is not limited to, two groups of components. One is a detection unit as a mechanical device that is actually installed at an observation point, and the other is an image data recording unit 2, a water level data recording unit 3, and a water level situation estimation data acquisition unit 4 that are virtually installed in an information processing device that processes various data based on the output of the detection unit. This system S may also be equipped with an additional group of components, a display unit 6 as a mechanical device that displays flood situation estimation data to grasp the damage situation at various points in real time.

[0024] 3 shows an example of connections including hardware of the present system S' according to one variation of the present system S according to the present embodiment. The present system S' is characterized in that each element is connected via a network N, and in particular, processing is possible even if the information processing device and the detection unit are in different positions or at different locations.

[0025] A plurality of detection units 1 are installed at observation points different from the points where the information processing devices are installed, and are connected to each element such as the information processing devices through a network N.

[0026] The image data recording unit 2, water level data recording unit 3, and water level condition estimation data acquisition unit 4, which are virtually installed within the information processing device, may be placed within a single information processing device, or may be connected to multiple information processing devices via a network N and distributed within one of them.

[0027] The display unit 6 can be provided in an information processing device, such as a portable information terminal such as a smartphone, and a plurality of display units can be provided in the information processing device directly or indirectly via the network N.

[0028] (Detection unit 1) As described above, the detection unit 1 in the present system S can operate at the same observation point as the information processing device, but it is preferable to install a plurality of detection units at observation points different from the point where the information processing device is installed. By providing at least one detection unit at each of the plurality of observation points, it is possible to perform output according to the surrounding situation. A specific example of the detection unit 1 is a so-called video camera, and it is preferable that the camera can output the surrounding situation as an electric signal, preferably image data.

[0029] Moreover, it is preferable that the detection unit 1 is connected to an electric communication line, i.e., a so-called network N, so that the image data can be output to an external device directly or indirectly via another device. This has the advantage that it is possible to easily synchronize a plurality of detection units, and to accurately obtain a wide range of information. The detection unit 1 may be any device other than a video camera, such as a so-called smartphone with a camera, a tablet terminal with a camera, or a high-speed camera, as long as it can obtain the situation at the observation point as image data as described above and output it directly or indirectly to an external device.

[0030] The more detection units 1 installed at multiple different observation points, the more accurately it will be possible to estimate water level conditions over a wider area. However, it is preferable to install multiple units along river or bay coastlines where levee breaches and overtopping are expected, i.e., along areas expected to be the origins of flooding, so that the water level condition estimation data described below will match the actually observed data as closely as possible, improving the estimation accuracy.

[0031] In this system, when a video camera is used for the detection unit 1, it is preferable that the camera is capable of high-speed processing. For example, if a camera capable of taking 200 images per second is used, the time increment Δt in the extended Kalman filter FEM can be set small, which is advantageous in terms of improving the estimation accuracy of the water level and flow velocity.

[0032] In this system, an information processing device for processing image data may be attached to the detection unit 1. By recording the image data acquired by the detection unit 1 and further recording the water level data, the processing load of the centrally managed information processing device can be reduced.

[0033] In this system, the image data recording unit 2, water level data recording unit 3, and water level condition estimation data acquisition unit 4 are virtually installed in an information processing device, so-called a computer. More specifically, the calculation program, estimation program, and various data are recorded in a recording medium 7 such as a hard disk of the computer, and are read into a volatile recording medium 8 such as a memory as necessary and executed to form virtual units.

[0034] In this system, the image data recording unit 2, water level data recording unit 3, and water level status estimation data acquisition unit 4 may be configured to be installed inside the same computer, but they may also be configured to be connected via an electric communication line, i.e., a so-called network N, or other means that allow data to be exchanged and arranged separately and distributed. Therefore, the example in Fig. 3 shows, as an example, a configuration in which a computer having an image data recording unit 2 installed inside is directly connected near the detection unit 1, a configuration in which a computer having an image data recording unit 2 and a water level data recording unit 3 installed inside are directly connected, and further, a configuration in which another computer is installed at a location (remote location) away from the detection unit (observation point) via the network N and the image data recording unit 2, water level data recording unit 3, and water level status estimation data acquisition unit 4 are installed inside it, and shows examples in which each unit is executed as necessary.

[0035] The computer in this system is not limited as long as it has the above functions, but it is preferable that it includes, for example, a central processing unit (CPU), which is a component of a general computer, a non-volatile recording medium such as a hard disk or flash memory, a volatile recording medium such as memory, a bus connecting these, input devices such as a keyboard or mouse, a display device such as a monitor, etc., but is not limited to these.

[0036] The computer may be a so-called notebook computer or desktop computer, but if high-speed calculation is required, it may be a so-called workstation or supercomputer, or may be a mobile information terminal, specifically a so-called smartphone or tablet terminal, which has become increasingly popular in recent years, in consideration of portability. A typical mobile information terminal is a terminal cover in which the above-mentioned CPU, display device, and other components are stored and integrated, and a sensor is placed on the display device to make it a so-called touch panel, which is very easy to handle and is preferable. In the case of a mobile information terminal, a program for executing the present method is recorded and displayed as a so-called app in the mobile information terminal, and the present method can be easily executed by starting the app. Furthermore, it can also serve as the display unit 6 described later.

[0037] (Water level estimation program) As described above, the image data recording unit 2, water level data recording unit 3, and water level status estimation data acquisition unit 4 are established as virtual units by storing a program in a recording medium 7 such as a computer's hard disk, reading it into a volatile recording medium 8 such as a memory, and executing it. That is, in this system, a water level status estimation program (hereinafter referred to as "this program") is stored in a computer, and by executing this, the image data recording unit 2, water level data recording unit 3, and water level status estimation data acquisition unit 4 are virtually provided in the computer, making it possible to perform each process.

[0038] Focusing on this program, a more detailed explanation will be given with reference to Fig. 4. This program causes a computer to execute the steps of (S1) acquiring calculation conditions for the extended Kalman filter FEM by inputting analysis conditions such as mesh data and observation point data, and estimation conditions including at least one of estimation variables such as flow velocity and water level fluctuation, (S2) recording image data, (S3) recording water level data based on the image data, (S4) creating water level status estimation data by the extended Kalman filter FEM based on the estimation conditions and water level data, and (S5) displaying the water level status estimation data. Note that, in executing this program, the above step (S1) is performed by the estimation condition acquisition unit 5, step (S2) is performed by the image data recording unit 2, step (S3) is performed by the water level data recording unit 3, and step (S4) is performed by the water level status estimation data acquisition unit 4.

[0039] (Water level estimation method) As is clear from the above description, when this program is executed on an information processing device, i.e., a computer, it can realize a water level estimation method (hereinafter also referred to as "this method") having the steps of: (S1) acquiring analysis conditions such as mesh data and observation point data, and estimation conditions including estimation variables such as flow velocity and water level fluctuation, (S2) recording image data, (S3) recording water level data based on the image data, (S4) creating water level estimation data using an extended Kalman filter FEM based on the estimation conditions and the water level data, and (S5) displaying the water level estimation data.

[0040] Next, the processing of each recording unit and acquiring unit will be specifically explained in correspondence with the steps of the present method.

[0041] (S1: Obtain estimation conditions) First, this method includes a step (S1) of acquiring calculation conditions for the extended Kalman filter FEM by inputting analysis conditions such as mesh data and observation point data, and estimation conditions including estimation variables such as flow velocity and water level fluctuation, and by executing this step, an estimation condition acquisition unit 5 is virtually provided. These data are recorded in a recording medium 7.

[0042] Examples of analysis conditions, which are estimation conditions, include mesh data (total number of nodes, total number of elements, coordinate values ​​of each node on the xy plane), observation point data (coordinate points of observation points and observed water level data), time increment Δt, number of time steps, gravitational acceleration, reference water depth (data in the z direction), parameters set to prevent numerical calculations from diverging (e.g., ramping parameters), initial values ​​of estimation variables, initial values ​​of the estimation error covariance matrix, system error covariance matrix, observation error covariance matrix, boundary conditions related to flow velocity and water level given at the boundaries of the area to be analyzed, etc. It is important to use at least one of these, and preferably multiple.

[0043] The estimated variables are also one of the estimation conditions required in the extended Kalman filter FEM, and include information on the flow velocity and water level fluctuation, for example. The estimated variables are also sometimes called observed variables.

[0044] The estimated conditions can be inputted and recorded in a computer using an input device such as a keyboard, a mouse, or a scanner. Information may also be read out from an optical recording medium such as a CD (Compact Disc) or a DVD (Digital Versatile Disc), a magneto-optical recording medium such as an MO (Magneto-Optical disk), a tape medium, a magnetic recording medium, or a semiconductor memory. The estimated conditions may also be acquired from an external storage medium such as a USB memory, and these devices may be acquired directly via an interface that allows the estimated conditions to be inputted to an information processing device, such as a USB (Universal Serial Bus) or Bluetooth (registered trademark), or may be acquired via a network or a communication line.

[0045] Among the analysis conditions that serve as the initial values ​​for the extended Kalman filter FEM, the mesh data specifically includes information indicating the x, y position on the map for a certain area on the map (observation target area) where the water level situation is to be estimated, as well as information z of the depth from the reference plane at a certain x, y point. The mesh can be set by dividing the observation target area using a finite element mesh. When dividing the mesh, it is preferable to set it for the observation point where the detection unit 1 is installed and the area including the river or bay to be observed. Note that the area including the river or bay is sometimes simply called the water area.

[0046] Among the analysis conditions that are the initial values ​​of the extended Kalman filter FEM, the observation point data is information related to the observation point where the detector 1 is installed, specifically, information such as the measured water level and the measured flow velocity in the x and y directions. Although calculations using the extended Kalman filter FEM are possible using only the measured water level or the measured flow velocity in the x and y directions, previous studies have shown that calculations can be performed with higher accuracy when using "measurement values ​​of both (measurement values ​​of flow velocity in the x and y directions and measured water level)" compared to when using "measurement values ​​of only the measured water level". In the estimation analysis of the flow field using the extended Kalman filter FEM, it has been confirmed that the influence of the observation values ​​is transmitted to the downstream side, so it is preferable that the observation point data includes water level information on the upstream side of the water body in the observation target area.

[0047] (S2: Image data recording) This method also includes a step of recording image data (S2). By executing this step, an image data recording unit 2 is virtually provided.

[0048] In this step, it is preferable to record calm time image data, which is image data during calm times, and estimated time image data, which is image data during times other than calm times. "Calm time image data" is data including information on the state before flooding or inundation occurs, and can be used to estimate the water level in a calm state at a time when flooding or inundation damage has not occurred. Also, "estimated time image data" is data including information other than calm times at at least one observation point of the multiple detection units, specifically, data including information at the observation point from immediately after flooding or inundation occurs, and by using it to estimate the water level from the time when flooding or inundation occurs, it becomes possible to grasp the damage situation in real time.

[0049] (S3: Water level data recording) This method also includes a step (S3) of recording water level data based on image data, and by executing this step, a water level data recording unit 3 is virtually provided.

[0050] In addition, in the step of recording the water level data, it is preferable to perform at least one of binarization (S31) and noise processing (S32) on the image data recorded in the image data recording step (S2). Note that, in the step of recording the water level data (S3), if both the binarization (S31) and noise processing (S32) are performed, it is preferable to perform the noise processing (S32) after the binarization (S31).

[0051] In addition, the image data for recording the water level data used to create the water level estimation data described later can be used to grasp the water level situation in real time by using the estimation time image data. In this case, it is preferable to add synchronized time information to the image data in the multiple detection units. This makes it easy to set the time when flooding or inundation is detected at at least one observation point of the multiple detection units as the water level estimation reference time. Examples of a method for adding time information include a method of shooting a synchronized signal, such as a synchronized blinking light or an accurate radio clock, at a fixed position in a part of the camera's field of view that is not the water surface. This makes it possible to record the time information as a part of the image. If the signal position in the captured image does not change, the time can be obtained by the same process in each camera.

[0052] (S31: Binarization processing) Here, the binarization process (S31) will be described. Identification of the water surface position by the binarization process (S31) involves taking the difference between two consecutive image data. Specifically, in an image including the water surface obtained by a detection unit with a fixed angle of view, the difference between each pixel information that changes over time can be taken to extract the water surface position and remove the background. Then, a brightness threshold is determined for the difference information of the two extracted image data, and binarization is performed. Note that, for binarization, it is preferable to convert the image including the water surface or the difference information into a grayscale in advance, since this reduces the number of variables for determining the brightness threshold. In addition, when using grayscaled difference information, for example, the brightness threshold is preferably determined so as to distinguish between image noise caused by a slight brightness difference in the background and a slight brightness difference caused by changes in the water surface. The brightness threshold may be determined manually or automatically in an information processing device by machine learning, and any means may be used. Then, binarization is possible by setting values ​​below the threshold as black = 0 and values ​​above the threshold as white = 1. As described above, using the brightness after grayscaling for binarization is advantageous in terms of reducing the load on the computer, but color information may be used individually, and the parameters used for binarization are not limited.

[0053] (S32: Noise processing) In addition, in the noise processing (S32), when using image data after binarization processing, it is preferable to repeat a process (dilation process) in which a certain pixel in the image data is focused on and the focused pixel is made white if there is at least one white pixel among the eight pixels adjacent to the focused pixel, and a process (erosion process) in which the focused pixel is made black if there is at least one black pixel adjacent to the focused pixel. In addition, in this noise processing, it is preferable to perform a process (closing process) in which the erosion process is performed multiple times and then the expansion process is performed the same number of times. For example, the erosion process is performed twice and then the expansion process is performed twice (erosion → erosion → expansion → expansion). On the other hand, it is also preferable to perform a process (opening process) in which the expansion process is performed multiple times and then the erosion process is performed the same number of times (for example, expansion → expansion → erosion → erosion). If the opening process is performed on the image data that has been binarized and then the closing process is performed, there is an advantage that white dot noise can be removed. Note that, in the case of using color information, noise can also be removed by repeating the expansion process and the erosion process.

[0054] (S33: Get water level) Water level data at the observation point can be obtained by comparing the image after noise processing with an image of an object where the water level can be confirmed before binarization, such as a water level board with a water level scale engraved on it. Either calm time image data or estimated time image data can be used as the image where the object where the water level can be confirmed is used, but it is preferable to use calm time image data in order to obtain accurate water level data without being affected by binarization or noise processing. This step is the water level acquisition step (S33).

[0055] (S4: Obtaining water level estimation data) In addition, this program has a step (S4) of acquiring water level condition estimation data by an extended Kalman filter FEM based on analysis conditions such as mesh data and observation point data, estimation conditions including estimation variables such as flow velocity and water level fluctuation, and water level data, and a water level condition estimation data acquisition unit 4 that executes this step is virtually provided.

[0056] In particular, this method has the advantage that it can obtain highly accurate water level estimation data over a wide range of water areas by estimating the water level using an extended Kalman filter FEM that takes into account the movement of the water shoreline. Here, "movement of the water shoreline" refers to creating a mesh of the entire observation area in advance, dividing the area into a wet area with water and a dry area without water, and providing minute water depths to elements that include nodes on the water shoreline, which makes it possible to analyze shallow water flows that take into account the movement of the water shoreline.

[0057] In addition, in this method, the "water level condition estimation data" preferably includes water area estimation data and water depth estimation data. The water area estimation data is data that indicates a planar extent, and more specifically, it is preferable that the data includes "information on planar coordinates x, y". On the other hand, the water depth estimation data is data that indicates the water depth at a certain x, y point, and includes "information on the depth z direction" at that x, y point.

[0058] In addition, in the step (S4) of creating or acquiring water level estimation data, the water area estimation data and the water depth estimation data can be obtained by performing a step (S41) of calculating the Kalman gain matrix, followed by a step (S42) of estimating the water area and a step (S43) of estimating the water depth.

[0059] (S41, S42: Kalman gain matrix calculation, water area estimation) Incidentally, the Kalman gain matrix calculation (S41) in the extended Kalman filter FEM is preferably performed so as to minimize the trace norm of the estimation error covariance matrix. The calculation and the calculation order for minimizing the trace norm of the estimation error covariance matrix are not limited, but can be performed using, for example, the specific procedure disclosed in Transactions of the Japan Society of Mechanical Engineers, Vol. 82, No. 835, p. 15, 2016.

[0060] Water area estimation (S42) can be output by calculating the system equation and then correcting it using the Kalman gain matrix shown in the following equation (4) and the observed value shown in the following equation (2). Specifically, the system equation describing the water area estimation can be expressed as the following equation (1), and the observation equation describing the observed value can be expressed as the following equation (2). Equation (1) below is expressed by an equation in which the shallow water long wave equation (motion equation and continuity equation including advection term (nonlinear term)) is discretized in the spatial direction based on the finite element method and also in the time direction using the finite difference method.

[0061]

number

[0062]

number

[0063] Here, {φ} is the state vector representing the flow velocity and water level, n is the time step, [A] is the state transition matrix (the reference water depth h and the gravitational acceleration g are also included in the matrix [A]), [Γ] is the driving matrix, {q} is the system noise vector, {z} is the observation value vector, [H] is the observation matrix, and {r} is the observation noise vector.

[0064] In addition, the estimation error covariance matrix before assimilation is expressed by the following equation (3), the estimation error covariance matrix after assimilation is expressed by the following equation (4), the Kalman gain matrix is ​​expressed by the following equation (5), the calculation of the estimated value is expressed by the following equation (6), and the calculation of the optimal estimated value is expressed by the following equation (7).

[0065]

number

[0066]

number

[0067]

number

[0068]

number

[0069]

number

[0070] Here, [P (-) ] is the error covariance matrix before assimilation, [F] is the state transition matrix after linear approximation, [Q] is the system error covariance matrix, [K1] is the Kalman gain matrix, [R] is the observation error covariance matrix, [P (+) ] is the error covariance matrix after assimilation (prediction error covariance matrix), [P (-) ] is the error covariance matrix before assimilation (estimated error covariance matrix), T is the transpose symbol, {φ^ (-)} is the estimated value vector, {φ^ (+)} indicates the optimal estimate vector. Here, ^ (hat) indicates the estimate. The water line processing is performed before equation (6), and in the element (mesh) that contains the node that represents the water line, minute water depth is given to the adjacent node. By calculating the left side of equation (6) (estimate vector at n+1 step) using the estimate vector at n step corrected by the minute water depth on the right side of equation (6), it becomes possible to calculate the water level situation taking into account the movement of the water line.

[0071] The procedure for estimating and analyzing the water level situation using the extended Kalman filter FEM, taking into account the movement of the waterfront line, can be, for example, the following steps 1 to 9.

[0072] (Water level analysis procedure) 1. Input the calculation conditions such as the analysis mesh, the time step Δt for the water level estimation analysis of the water level situation, the time interval ΔT for the water level measurement, the reference water depth, and the initial conditions (flow velocity, water level fluctuation). It is preferable that Δt and ΔT are the same value in the water level estimation analysis of the water level situation, but if Δt and ΔT are different, adjust so that the measured water level value is given for every time step Δt. 2. Water level measurement {z n+1}. 3. The estimation error covariance matrix [P (-) n+1 Calculate (Equation (3)) 4. The Kalman gain matrix [K1 n+1 Calculate (Equation (4)) 5. Prediction error covariance matrix [P (+) n+1 Calculate (Equation (5)) 6. In elements (mesh) that contain nodes on the waterfront, minute water depths are assigned to adjacent nodes. 7. Estimated flow velocity and water level fluctuation at a certain time step n+1 {φ^ (-) n+1} is calculated (Equation (6)). 8. The measured water level at time step n+1 is {z n+1}, the optimal estimate of flow velocity and water level fluctuation {φ^ (+) n+1} is calculated (Equation (7)). 9. Update the time step (add a time step Δt to the current time) and go back to step 2 above.

[0073] The specific calculation sequence of the Kalman filter FEM that does not take into account the movement of the waterfront can be performed, for example, using the specific procedure disclosed in Transactions of the Japan Society of Mechanical Engineers, Vol. 82, No. 835, 2016, p. 15.

[0074] (S43: Water depth estimation) On the other hand, water depth estimation (S43) can be performed using the same process as water area estimation. After calculating the above equation (1), the result can be output by correcting it using the Kalman gain matrix shown in the above equation (4) and the observed value shown in the above equation (2).

[0075] (S5: Display of estimated water level data) This method also includes a step (S5) of displaying the water level status estimation data. This refers to displaying the water level status estimation data on a display unit such as a monitor or display that a computer necessarily has, or a display of a smartphone.

[0076] The display unit 6 may be directly connected via an interface that enables the output of water level estimation data from an information processing device having the water level estimation data acquisition unit 4, or may be connected via a network or communication line. EXAMPLES

[0077] Here, the water level estimation method according to the above embodiment was processed using actual data to confirm its effectiveness. A specific description will be given below.

[0078] (Evaluation experiment 1: Water level data recording section) A gate was installed at one end of a waterway (10 m long, 0.4 m wide) to hold the water mass. Detection units capable of measuring the water level were installed 0.25 m, 1 m, and 2 m from the gate. The detection unit consisted of a water level board installed in the waterway and a high-speed camera (taking 200 images per second).

[0079] The gate was opened momentarily to release the water mass into the waterway. The time when the gate was opened was set to zero seconds, and measurements were started simultaneously for all detection parts. After the measurement started, the image data was recorded in a computer with an image data recording unit. Two consecutive images (example shown in Figure 7) were read out from the recorded image data, and images of the same range near the water level plate were cut out from both images. Next, the cut-out images were each grayscaled, and both grayscaled images were binarized using the same brightness threshold. After noise processing, the water level engraved on the water level plate was compared with the water surface position clarified by noise processing to obtain the water level.

[0080] Figure 8 shows the results of recording the change in water level over time at each point. Even in a situation where the gate was opened momentarily and a large amount of water splashed onto the water surface, it was confirmed that the peak water level became lower and the time it took to reach the peak became slower as the gate moved away from the water.

[0081] (Evaluation experiment 2-1: Water level estimation data acquisition section) The accuracy of estimating flow fields in shallow water areas using the Kalman filter finite element method was verified by changing the settings of observation variables in the water area and the number and locations of observation points.

[0082] The verification used water level information provided by tide gauge stations to estimate the water level situation in Tokyo Bay, and then compared it with water level information provided by tide gauge stations that were not used in the water level estimation. The water level data used was from three locations: Tokyo (2313), Chiba (3305), and Yokohama New Port (3308) (numbers in parentheses are all tide gauge station registration numbers). The water level data used for comparison after estimating the water level situation was from the tide gauge station in Yokosuka (3301). In this study, one type of mesh data was used, and one set of observation point data was also adopted.

[0083] The water depth estimation data was obtained using chart data, and the calculation parameters were as shown in Table 1.

[0084] [Table 1]

[0085] Figure 9 shows the results of the flow estimation analysis, i.e., the estimation of water level conditions. The estimated water level conditions data (estimated value) in Yokosuka (3301), estimated using water level data from three locations, Tokyo (2313), Chiba (3305), and Yokohama New Port (3308), was confirmed to be in good agreement with the water level data (observed value) observed in Yokosuka (3301).

[0086] In addition, instead of using the water level data from the three locations mentioned above, a flow estimation analysis was performed using a numerical analysis method that includes the effect of high numerical viscosity, using water level data obtained using conventional observation technology. The results are shown in Figure 10, but the estimated water level data in Yokosuka did not match the results.

[0087] As a result of estimating the water level situation, we were able to estimate the distribution of flow velocity and water level fluctuations within Tokyo Bay (Figure 11).

[0088] This example shows that using the x- and y-direction flow velocities and water level fluctuations on a two-dimensional plane as observed variables produces more accurate flow field estimation results than using only the water level fluctuations as observed variables. It also shows that it is more important to place the observation points upstream of the flow than downstream.

[0089] (Evaluation experiment 2-2: Water level estimation data acquisition section) Using the water level data obtained from image analysis, we conducted an estimation analysis of the water level situation for the dam break problem, taking into account the moving boundary of the waterfront (Figure 12).

[0090] [Table 2]

[0091] The results are shown in Figures 13 and 14. Figure 13 shows the results for (x=1.1m, y=0.2m), and Figure 14 shows the results for (x=1.8m, y=0.2m). From these results, it was confirmed that the water level estimation results were close to the water level data obtained from image analysis.

[0092] As described above, this example confirmed the usefulness of the present method. [Industrial Applicability]

[0093] The present invention has industrial applicability as a water level estimation program, a water level estimation system, and a water level estimation method. Flood areas can be estimated with high accuracy and in real time, and by preparing the water level estimation system developed according to the present invention, evacuation can be appropriately guided in the event of river flooding or flood disasters caused by, for example, torrential rain. [Explanation of symbols]

[0094] S, S' - This system 1. Detection section 2. Image data recording section 3. Water level data recording section 4. Water level estimation data acquisition section 5... Estimated condition acquisition section 6...Display section

Claims

1. A method for recording image data from an observation point located upstream of a flow in a water body; recording water level data based on the image data; A step of creating water level state estimation data by an extended Kalman filter finite element method based on the estimation conditions and the water level data, The extended Kalman filter finite element method uses the following equation (1) as a system equation: A step of calculating an estimated error covariance matrix [P (−) n+1 ] at the time step n+1 following the time step n based on the state transition matrix [F n ] after linear approximation at the time step n and the prediction error covariance matrix [P (+) n ]; calculating a Kalman gain matrix [K 1 n+1 ] at the next time step n+1 based on the estimated error covariance matrix [P (−) n+1 ]; calculating a prediction error covariance matrix [P (+) n+1 ] at the next time step n+1 based on the estimation error covariance matrix [P (−) n+1 ] and the Kalman gain matrix [K 1 n+1 ]; Further, the method includes a step of providing minute water depths to adjacent junctions in elements including nodes of waterfront boundaries consisting of wet and dry regions; Furthermore, an estimate vector {φ^ (-) n+1} at the next time step n+1 is calculated based on the state transition matrix [A] and the estimate vector {φ^ (-) n} at the time step n, and an optimal estimate vector {φ^ (+) n+1} at the next time step n+1 is calculated based on the estimate vector {φ^ (-) n+1}, the water level data, and the Kalman gain matrix [K 1 n+1 ]. Method for estimating water level conditions. [0010] where {φ} is the state vector, n is the time step, [A] is the state transition matrix, [Γ] is the driving matrix, and {q} is the system noise vector.

2. 2. The method for estimating water level conditions according to claim 1, wherein the step of recording the water level data includes performing at least one of a binarization process and a noise removal process on the image data.

3. The water level estimation method according to claim 1 , wherein the water level estimation data includes water area estimation data and water depth estimation data.

4. The step of recording the image data includes recording image data of a state before the inundation or flooding occurs and image data at the time of estimation, 2. The method for estimating water level conditions according to claim 1, wherein the step of recording the water level data includes recording the water level data based on the image data at the time of estimation.

5. The water level condition estimating method according to claim 1 , further comprising a step of displaying the water level condition estimation data.

6. A plurality of detection units installed at different observation points on the upstream side of the flow of a water body; an image data recording unit that records image data based on an output from the detection unit; a water level data recording unit that records water level data based on the image data; a water level state estimation data acquisition unit that acquires water level state estimation data by an extended Kalman filter finite element method based on the estimation conditions and the water level data; In the water level state estimation data acquisition unit, the following equation (1) is used as a system equation in the extended Kalman filter finite element method, Calculate an estimated error covariance matrix [P (−) n+1 ] at the time step n+1 following the time step n based on the state transition matrix [F n ] after linear approximation and the prediction error covariance matrix [P (+) n ] at the time step n; Calculating a Kalman gain matrix [K 1 n+1 ] at the next time step n+1 based on the estimation error covariance matrix [P (−) n+1 ]; calculating a prediction error covariance matrix [P (+) n+1 ] at the next time step n+1 based on the estimation error covariance matrix [P (−) n+1 ] and the Kalman gain matrix [K 1 n+1 ]; Furthermore, in elements that include nodes on the waterfront boundary consisting of wet and dry areas, minute water depths are given to adjacent junctions, Furthermore, an estimate vector {φ^ (-) n+1} at the next time step n+1 is calculated based on the state transition matrix [A] and the estimate vector {φ^ (-) n} at the time step n, and an optimal estimate vector {φ^ (+) n+1} at the next time step n+1 is calculated based on the estimate vector {φ^ (-) n+1}, the water level data, and the Kalman gain matrix [K 1 n+1 ]. Water level estimation system. [0010] where {φ} is the state vector, n is the time step, [A] is the state transition matrix, [Γ] is the driving matrix, and {q} is the system noise vector.

7. A water level estimation program that causes a computer to execute the water level estimation method according to any one of claims 1 to 5.

Images