JP2019139455A - Runoff analyzer and runoff analysis method - Google Patents

Runoff analyzer and runoff analysis method Download PDF

Info

Publication number
JP2019139455A
JP2019139455A JP2018021427A JP2018021427A JP2019139455A JP 2019139455 A JP2019139455 A JP 2019139455A JP 2018021427 A JP2018021427 A JP 2018021427A JP 2018021427 A JP2018021427 A JP 2018021427A JP 2019139455 A JP2019139455 A JP 2019139455A
Authority
JP
Japan
Prior art keywords
virtual
mesh
storage amount
water level
runoff
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2018021427A
Other languages
Japanese (ja)
Other versions
JP7013270B2 (en
Inventor
理 山中
Satoru Yamanaka
理 山中
由紀夫 平岡
Yukio Hiraoka
由紀夫 平岡
弘明 古米
Hiroaki Komai
弘明 古米
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Toshiba Infrastructure Systems and Solutions Corp
Original Assignee
Toshiba Corp
Toshiba Infrastructure Systems and Solutions Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Corp, Toshiba Infrastructure Systems and Solutions Corp filed Critical Toshiba Corp
Priority to JP2018021427A priority Critical patent/JP7013270B2/en
Publication of JP2019139455A publication Critical patent/JP2019139455A/en
Application granted granted Critical
Publication of JP7013270B2 publication Critical patent/JP7013270B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Sewage (AREA)
  • Alarm Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

To provide a runoff analyzer and a runoff analysis method with which it is possible to obtain a runoff analysis result that is close to reality while using a simple analysis model.SOLUTION: A runoff analyzer 1 comprises a rainfall distribution input unit 101 and a runoff analysis unit 103. The rainfall distribution input unit accepts, as its input, rainfall distribution information for each mesh of an area being diagnosed. The runoff analysis unit calculates a retained quantity in each mesh on the basis of the rainfall distribution information. The runoff analysis unit calculates a hypothetical water level in a mesh where the retained quantity is less than the maximum retained quantity of a hypothetical conduit line, calculates a hypothetical water level in a mesh where the retained quantity is greater than or equal to the maximum retained quantity of the hypothetical conduit line and less than an upper-limit retained quantity that is the sum of the maximum retained quantity of the hypothetical conduit line and the maximum retained quantity of a hypothetical manhole, calculates a hypothetical water level in a mesh where the retained quantity is greater than or equal to the upper-limit retained quantity, and calculates various quantities needed for the evaluation of flooding risk in each mesh using the hypothetical water levels calculated in accordance with the retained quantities.SELECTED DRAWING: Figure 1

Description

本発明の実施形態は、流出解析装置及び流出解析方法に関する。   Embodiments described herein relate generally to an outflow analysis device and an outflow analysis method.

近年、局地的かつ短時間に降る大雨(以下「局所豪雨」という。)が多発しており、メディアはこの局所豪雨のことを、いつどこで発生するか分からないという意味で「ゲリラ豪雨」という言葉で表現し、今ではこの言葉も広く世間一般に認知されてきた。局所豪雨による典型的な被害として、都市の内部で水が溢れる内水氾濫が頻発するようになってきた。これまで行政は、洪水を未然に防ぐため、築堤、河道掘削、護岸整備やダム建設など、主に大規模河川の増水や決壊による洪水を想定した対策を講じてきた。河川の氾濫は外水氾濫と呼ばれ、従来は外水氾濫に対する対策が重点的に行われてきたが、今後は、内水氾濫も考慮した対策が重要になってくると考えられる。実際、洪水被害(外水氾濫、内水氾濫)を被害額でみると、内水氾濫の被害額は全国の被害総額の約半分を占め、東京都では都の被害総額の90%を超えている。このように、堤防の整備が比較的進んだ都市部では、内水氾濫が新たな課題となっている。   In recent years, there has been a lot of localized and short-term heavy rainfall (hereinafter referred to as “local heavy rain”), and the media refers to this local heavy rain as “guerrilla heavy rain” in the sense that it does not know when and where it occurs. Expressed in words, these words are now widely recognized by the general public. As a typical damage caused by local heavy rain, inundation floods occur frequently. In the past, the government has taken measures to prevent flooding, mainly in anticipation of flooding due to flooding or breakage of large rivers, such as embankments, river channel excavation, revetment construction and dam construction. River inundation is called outside water inundation. Conventionally, countermeasures against outside water inundation have been focused on, but in the future, it will be important to take into account inland water inundation. In fact, when we look at flood damage (outside water inundation, inland water inundation) in terms of damage amount, the damage amount of inland flooding accounts for about half of the total damage amount in the country, exceeding 90% of the total damage amount in Tokyo Yes. In this way, inundation is a new issue in urban areas where the development of dykes is relatively advanced.

このような都市全体での浸水リスクを評価する手法として分布型流出解析が広く用いられている。分布型流出解析は、ある地域の土地の利用形態や標高などの地形情報、下水管路の敷設状況などの土木情報に基づき、水文学的なモデルと水理学的なモデルとを適宜併用して構築される流出解析モデルを用いて降雨の流れ(以下「流出」という。)を追跡する流出解析方法である。具体的には、解析対象の地域が分割された複数のメッシュについて、上記の流出解析モデルと各メッシュの雨量データとを適用することによってメッシュごとの流出状況を解析することができる。しかしながら、実際には、取得可能な雨量データがこのような流出解析モデルに必ずしも適した態様で得られるとは限らず、また、土木情報についても必ずしも必要な全ての情報が得られるとも限らない。そのため、取得可能な雨量データや土木情報を流出解析モデルに適応させるために労力を要していた。このような課題を解決するために、流出解析モデルを用いる際の手間を低減する手法がいくつか提案されている。   Distributed runoff analysis is widely used as a method for assessing such inundation risks throughout the city. Distributed runoff analysis is based on the use of land in a certain region, topographical information such as elevation, and civil engineering information such as the status of sewage pipe construction, and a combination of hydrological and hydraulic models as appropriate. This is a runoff analysis method that tracks the flow of rainfall (hereinafter referred to as “runoff”) using the constructed runoff analysis model. Specifically, the runoff situation for each mesh can be analyzed by applying the runoff analysis model and the rainfall data of each mesh to a plurality of meshes obtained by dividing the analysis target area. In practice, however, acquirable rainfall data is not always obtained in a mode suitable for such a runoff analysis model, and all necessary information is not necessarily obtained for civil engineering information. For this reason, labor has been required to adapt the available rainfall data and civil engineering information to the runoff analysis model. In order to solve such a problem, several methods for reducing labor when using the outflow analysis model have been proposed.

その一方で、土木情報を用いずに流出解析を行う手法も提案されている。例えば、各メッシュに仮想的なタンク(以下「仮想タンク」という。)を割り当て、仮想タンク間の水の収支のみで流出解析を行う解析モデル(以下「仮想タンクモデル」という。)が提案されている。しかしながら、このような仮想タンクモデルでは、土木情報を用いていないため、動水勾配や径深、流積を計算するための水位を適切に定義することができず、必ずしも現実の流出に即した解析が行えない可能性があった。   On the other hand, a method for performing runoff analysis without using civil engineering information has also been proposed. For example, an analysis model (hereinafter referred to as “virtual tank model”) in which a virtual tank (hereinafter referred to as “virtual tank”) is assigned to each mesh and the outflow analysis is performed only by the balance of water between the virtual tanks has been proposed. Yes. However, in such a virtual tank model, since civil engineering information is not used, it is not possible to appropriately define the water level for calculating the hydrodynamic gradient, the depth of the depth, and the flow product. There was a possibility that analysis could not be performed.

特開2009−8651号公報JP 2009-8651 A 特開2005−128838号公報JP 2005-128838 A 特許第4682178号公報Japanese Patent No. 4682178 特許第4185910号公報Japanese Patent No. 4185910 特許第4082686号公報Japanese Patent No. 4082686

本発明が解決しようとする課題は、簡易な解析モデルを用いつつ、より現実に即した流出解析結果を得ることのできる流出解析装置及び流出解析方法を提供することである。   The problem to be solved by the present invention is to provide an outflow analysis device and an outflow analysis method capable of obtaining a more realistic outflow analysis result while using a simple analysis model.

実施形態の流出解析装置は、降雨分布入力部と、流出解析部と、を持つ。降雨分布入力部は、複数のメッシュに分割された診断対象地域の降雨量を示す情報であって前記メッシュごとの降雨量を示す降雨分布情報を入力する。流出解析部は、前記降雨分布情報に基づいて各メッシュ間における流量の収支バランスを計算することにより各メッシュの貯留量を算出する。前記流出解析部は、算出された各メッシュの貯留量と、前記メッシュごとに設定された仮想管路の最大貯留量及び前記仮想管路と地表面とを接続する仮想マンホールの最大貯留量とに基づいて、前記貯留量が前記仮想管路の最大貯留量未満であるメッシュの仮想水位を以下の式(1)で算出し、前記貯留量が前記仮想管路の最大貯留量以上であり、かつ前記仮想管路の最大貯留量と前記仮想マンホールの最大貯留量との和である上限貯留量未満であるメッシュの仮想水位を以下の式(2)で算出し、前記貯留量が上限貯留量以上であるメッシュの仮想水位を以下の式(3)で算出し、前記貯留量に応じて算出された前記仮想水位を用いて各メッシュの浸水リスクの評価に必要な諸量を算出する。   The runoff analysis device of the embodiment has a rainfall distribution input unit and a runoff analysis unit. The rainfall distribution input unit inputs the rainfall distribution information indicating the rainfall amount for each mesh that is the rainfall amount in the diagnosis target area divided into a plurality of meshes. The runoff analysis unit calculates a storage amount of each mesh by calculating a flow balance between the meshes based on the rainfall distribution information. The outflow analysis unit calculates the storage amount of each mesh, the maximum storage amount of the virtual pipe set for each mesh, and the maximum storage amount of the virtual manhole connecting the virtual pipe and the ground surface. Based on the virtual water level of the mesh that the storage amount is less than the maximum storage amount of the virtual pipeline is calculated by the following equation (1), the storage amount is greater than or equal to the maximum storage amount of the virtual pipeline, and The virtual water level of the mesh that is less than the upper limit storage amount that is the sum of the maximum storage amount of the virtual pipeline and the maximum storage amount of the virtual manhole is calculated by the following equation (2), and the storage amount is equal to or higher than the upper limit storage amount The virtual water level of the mesh is calculated by the following equation (3), and various amounts necessary for evaluating the inundation risk of each mesh are calculated using the virtual water level calculated according to the storage amount.

実施形態の流出解析装置の構成の具体例を示す図。The figure which shows the specific example of a structure of the outflow analysis apparatus of embodiment. 実施形態の流出解析装置1が用いる流出解析モデルの概念を説明する図。The figure explaining the concept of the outflow analysis model which the outflow analysis apparatus 1 of embodiment uses. 仮想水位と、流積及び径深との関係を示す図。The figure which shows the relationship between a virtual water level, a pour, and a diameter depth. 仮想水位と、動水勾配及び管路勾配との関係を示す図。The figure which shows the relationship between a virtual water level, a dynamic water gradient, and a pipeline gradient.

以下、実施形態の流出解析装置及び流出解析方法を、図面を参照して説明する。   Hereinafter, the outflow analysis device and the outflow analysis method of the embodiment will be described with reference to the drawings.

図1は、実施形態の流出解析装置の構成の具体例を示す図である。流出解析装置1は、
バスで接続されたCPU(Central Processing Unit)やメモリや補助記憶装置などを備え、流出解析装置プログラムを実行する。流出解析装置1は、流出解析装置プログラムの実行によって降雨分布入力部101、流出係数算出部102、流出解析部103、浸水リスク評価部104、地図情報記憶部105、メッシュ面積設定部106、下水管理台帳データ記憶部107、仮想下水管径算出部108、仮想下水管勾配算出部109、仮想下水管粗度算出部110、マンホール総面積算出部111、容量算出部112及び診断結果表示部113を備える装置として機能する。なお、流出解析装置1の各機能の全て又は一部は、ASIC(Application Specific Integrated Circuit)やPLD(Programmable Logic Device)やFPGA(Field Programmable Gate Array)等のハードウェアを用いて実現されてもよい。流出解析装置プログラムは、コンピュータ読み取り可能な記録媒体に記録されてもよい。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ROM、CD−ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置である。流出解析装置プログラムは、電気通信回線を介して送信されてもよい。
FIG. 1 is a diagram illustrating a specific example of the configuration of the outflow analysis apparatus according to the embodiment. The spill analyzer 1
A CPU (Central Processing Unit) connected via a bus, a memory, an auxiliary storage device, and the like are provided, and an outflow analysis program is executed. The runoff analysis device 1 has a rainfall distribution input unit 101, runoff coefficient calculation unit 102, runoff analysis unit 103, inundation risk evaluation unit 104, map information storage unit 105, mesh area setting unit 106, sewage management by executing the runoff analysis device program. A ledger data storage unit 107, a virtual sewage pipe diameter calculation unit 108, a virtual sewage pipe gradient calculation unit 109, a virtual sewage pipe roughness calculation unit 110, a manhole total area calculation unit 111, a capacity calculation unit 112, and a diagnosis result display unit 113 are provided. Functions as a device. All or some of the functions of the outflow analysis device 1 may be realized using hardware such as an application specific integrated circuit (ASIC), a programmable logic device (PLD), or a field programmable gate array (FPGA). . The outflow analysis device program may be recorded on a computer-readable recording medium. The computer-readable recording medium is, for example, a portable medium such as a flexible disk, a magneto-optical disk, a ROM, a CD-ROM, or a storage device such as a hard disk built in the computer system. The outflow analysis device program may be transmitted via a telecommunication line.

降雨分布入力部101は、自装置に降雨分布情報を入力する。降雨分布情報は、雨量レーダ等から取得される、対象地域の降雨量の分布を示す情報である。具体的には、降雨分布情報は、メッシュ状に分割された対象地域について各メッシュの位置及び降雨量を示す情報の集合として得られる。   The rain distribution input unit 101 inputs the rain distribution information to its own device. The rainfall distribution information is information indicating the distribution of rainfall in the target area, which is acquired from a rainfall radar or the like. Specifically, the rainfall distribution information is obtained as a set of information indicating the position of each mesh and the rainfall amount for the target area divided into meshes.

例えば、降雨分布入力部101は、国土交通省が運用するXバンドMPレーダなどの雨量レーダによって取得された降雨分布情報を入力する。XバンドMPレーダによって取得される降雨分布情報は、250m四方のメッシュ状のデータとして取得される。なお、降雨分布情報は、XバンドMPレーダ以外によって取得されたものであってもよい。例えば、降雨分布情報は、アメッシュデータ(Xバンドレーダの一つ)であってもよいし、実用化が進められている垂直方向の降雨情報を計測することができるフェーズドアレイレーダによって取得されるデータであってもよい。流出解析の視点では、降雨分布情報は、Xバンドレーダのように細かいメッシュサイズで取得されることが望ましいが、1km四方のCバンドレーダによって取得されてもよい。また、降雨分布情報は、降雨レーダで計測されたレーダの反射強度を、レーダ方程式などを用いて換算することによって取得されるが、対象地域に設置された地上雨量計によって取得されたデータを補正したものであってもよい。   For example, the rainfall distribution input unit 101 inputs rainfall distribution information acquired by a rainfall radar such as an X-band MP radar operated by the Ministry of Land, Infrastructure, Transport and Tourism. The rainfall distribution information acquired by the X-band MP radar is acquired as 250m square mesh data. Note that the rainfall distribution information may be acquired by other than the X-band MP radar. For example, the rainfall distribution information may be Amesh data (one of X-band radars), or data acquired by a phased array radar that can measure rainfall information in the vertical direction that is being put to practical use. It may be. From the viewpoint of runoff analysis, it is desirable that the rainfall distribution information be acquired with a fine mesh size like an X-band radar, but may be acquired with a 1-km square C-band radar. Rainfall distribution information is obtained by converting the radar reflection intensity measured by the rain radar using a radar equation, etc., but corrects the data obtained by the ground rain gauge installed in the target area. It may be what you did.

流出係数算出部102は、降雨分布入力部101が入力した降雨分布情報に基づいて、各メッシュの流出係数を算出する。流出係数は、降雨量に対して地表を流下する雨水(有効降雨量)の割合を表す数値である。   The runoff coefficient calculation unit 102 calculates the runoff coefficient of each mesh based on the rain distribution information input by the rain distribution input unit 101. The runoff coefficient is a numerical value representing the ratio of rainwater (effective rainfall) that flows down the ground surface to the rainfall.

例えば、流出係数算出部102では、人工衛星によって取得される電子的な地図情報データであるGIS(Geographic Information System)データを用いて、降雨分布情報が取得されるメッシュごとの流出係数を算出する。例えば、流出係数算出部102は、以下のような方法でメッシュごとの流出係数を算出することができる。まず、流出係数算出部102は、GISデータが示す色や色の変化、形状などを用いて各土地の利用形態を推定する。例えば、土地の利用形態とGISデータが示す色(具体的には数値として定義される)との対応関係が予め定められても良い。例えば、緑色の部分は山、灰色の部分は道路、所定サイズ以下の四角い部分は屋根というように定義される。流出係数算出部102は、このような対応関係に基づいて、各メッシュにおける土地の利用形態の割合を算出する。流出係数算出部102は、土地の利用形態と流出係数との対応関係を予め記憶する。流出係数算出部102は、土地の利用形態に対応付けられた流出係数を、上記の割合を重みとしてメッシュごとに加算した値をそのメッシュの流出係数として算出する。   For example, the runoff coefficient calculation unit 102 uses the GIS (Geographic Information System) data, which is electronic map information data acquired by an artificial satellite, to calculate the runoff coefficient for each mesh from which rainfall distribution information is acquired. For example, the outflow coefficient calculation unit 102 can calculate the outflow coefficient for each mesh by the following method. First, the runoff coefficient calculation unit 102 estimates the use form of each land using the color, the change in color, the shape, and the like indicated by the GIS data. For example, the correspondence relationship between the land use form and the color (specifically defined as a numerical value) indicated by the GIS data may be determined in advance. For example, a green part is defined as a mountain, a gray part as a road, and a square part having a predetermined size or less as a roof. The runoff coefficient calculation unit 102 calculates the proportion of land use forms in each mesh based on such correspondence. The runoff coefficient calculation unit 102 stores in advance a correspondence relationship between land use patterns and runoff coefficients. The runoff coefficient calculation unit 102 calculates, as the runoff coefficient of the mesh, a value obtained by adding the runoff coefficient associated with the land use form for each mesh with the above ratio as a weight.

例えば、GISデータに基づいて、山の割合が0.6、道路の割合が0.3、屋根の割合が0.1としてあるメッシュの土地の利用形態の割合が算出されたとする。そして、山の流出係数が0.6、道路の流出係数が0.8、屋根の流出係数が1として予め定義されているとする。この場合、流出係数算出部102は、そのメッシュの流出係数は0.7(=0.6×0.6+0.3×0.8+0.1×1)となる。なお、土地の利用形態に関する情報が詳細に取得できない場合には、対象地域を山林部、田畑部、都市部などのおおまかな利用形態に分類し、これらのおおまかな利用形態ごとに流出係数を予め定義しておいてもよい。   For example, based on the GIS data, it is assumed that the proportion of the usage pattern of the mesh land with the mountain proportion being 0.6, the road proportion being 0.3, and the roof proportion being 0.1 is calculated. The mountain runoff coefficient is 0.6, the road runoff coefficient is 0.8, and the roof runoff coefficient is 1. In this case, the outflow coefficient calculation unit 102 has an outflow coefficient of 0.7 (= 0.6 × 0.6 + 0.3 × 0.8 + 0.1 × 1). If information on land use forms cannot be obtained in detail, the target area is classified into rough use forms such as forests, fields, and urban areas, and the runoff coefficient is preliminarily determined for each of these rough use forms. You may define it.

流出解析部103は、流出係数算出部102によって算出された各メッシュの流出係数を降雨分布情報が示すメッシュごとの降雨量に乗算することによって、各メッシュの有効降雨量を示す有効雨量情報を取得する。流出解析部103は、有効雨量情報を入力として流出解析を行うことにより、各メッシュの貯留量の変化を計算する。   The runoff analysis unit 103 obtains effective rainfall information indicating the effective rainfall amount of each mesh by multiplying the runoff coefficient of each mesh calculated by the runoff coefficient calculation unit 102 by the rainfall amount for each mesh indicated by the rainfall distribution information. To do. The runoff analysis unit 103 calculates the change in the storage amount of each mesh by performing the runoff analysis with the effective rainfall information as an input.

具体的には、流出解析部103は、各メッシュに後述の仮想下水管を設定した流出解析モデルを用いて流出解析を行う。この点において、本実施形態の流出解析モデルは、各メッシュに仮想タンクを設定した従来の仮想タンクモデル(例えば、特開2015−004245号公報参照。)と異なる。   Specifically, the outflow analysis unit 103 performs an outflow analysis using an outflow analysis model in which a virtual sewer pipe described later is set for each mesh. In this respect, the outflow analysis model of the present embodiment is different from a conventional virtual tank model (see, for example, JP-A-2015-004245) in which a virtual tank is set for each mesh.

浸水リスク評価部104は、流出解析部103によって解析された各メッシュの流出量に基づいて、メッシュごとの浸水有無や浸水の度合い(以下「浸水リスク」という。)を判定する。   Based on the outflow amount of each mesh analyzed by the outflow analysis unit 103, the inundation risk evaluation unit 104 determines whether or not the mesh is inundated and the degree of inundation (hereinafter referred to as “inundation risk”).

地図情報記憶部105は、磁気ハードディスク装置や半導体記憶装置などの記憶装置を用いて構成される。地図情報記憶部105は、流出係数算出部102による流出係数の算出に必要となる対象地域の地図情報を記憶する。例えば、GISデータであってもよいし、他の方法で取得された地図情報であってもよい。   The map information storage unit 105 is configured using a storage device such as a magnetic hard disk device or a semiconductor storage device. The map information storage unit 105 stores map information of the target area that is necessary for the calculation of the runoff coefficient by the runoff coefficient calculation unit 102. For example, it may be GIS data or map information obtained by another method.

メッシュ面積設定部106は、流出係数算出部102による流出係数の算出に必要となる各メッシュの面積を設定する。   The mesh area setting unit 106 sets the area of each mesh necessary for the calculation of the outflow coefficient by the outflow coefficient calculation unit 102.

下水道管理台帳データ記憶部107は、磁気ハードディスク装置や半導体記憶装置などの記憶装置を用いて構成される。下水道管理台帳データ記憶部107は、流出解析部103や浸水リスク評価部104の処理に必要となる下水道の物理緒元を示す情報(以下「下水道管理台帳データ」という。)を記憶する。   The sewer management ledger data storage unit 107 is configured using a storage device such as a magnetic hard disk device or a semiconductor storage device. The sewer management ledger data storage unit 107 stores information (hereinafter referred to as “sewer management ledger data”) indicating the physical specifications of the sewer necessary for the processing of the outflow analysis unit 103 and the inundation risk evaluation unit 104.

仮想下水管径算出部108は、下水道管理台帳データに基づいて、流出解析部103による流出解析に用いられる仮想下水管の管路径を算出する。   The virtual sewage pipe diameter calculation unit 108 calculates the pipe diameter of the virtual sewage pipe used for the outflow analysis by the outflow analysis unit 103 based on the sewer management ledger data.

仮想下水管勾配算出部109は、下水道管理台帳データに基づいて、流出解析部103による流出解析に用いられる仮想下水管の勾配を算出する。   The virtual sewage pipe gradient calculation unit 109 calculates the gradient of the virtual sewage pipe used for the outflow analysis by the outflow analysis unit 103 based on the sewer management ledger data.

仮想下水管粗度算出部110は、下水道管理台帳データに基づいて、流出解析部103による流出解析に用いられる仮想下水管の粗度を算出する。   The virtual sewage pipe roughness calculation unit 110 calculates the roughness of the virtual sewage pipe used for the outflow analysis by the outflow analysis unit 103 based on the sewer management ledger data.

マンホール総面積算出部111は、流出解析部103による流出解析に用いられるマンホール総面積を算出する。マンホール総面積は、対象地域に存在するマンホールのメッシュごとの総面積を表す。   The manhole total area calculation unit 111 calculates the total manhole area used for the outflow analysis by the outflow analysis unit 103. The total manhole area represents the total area of each manhole mesh existing in the target area.

容量算出部112は、流出解析部103や浸水リスク評価部104で用いられる最大貯留量及び上限貯留量を算出する。最大貯留量は、対象地域に存在する管路のメッシュごとの総容量である。上限貯留量は、対象地域に存在するマンホールのメッシュごとの総容量と、上記の最大貯留量との和である。   The capacity calculation unit 112 calculates the maximum storage amount and the upper limit storage amount used by the outflow analysis unit 103 and the inundation risk evaluation unit 104. The maximum storage amount is the total capacity for each mesh of pipelines existing in the target area. The upper limit storage amount is the sum of the total capacity of each manhole mesh existing in the target area and the maximum storage amount.

診断結果表示部113は、上記各機能部の処理結果に関する情報を、メッシュごとに表示する。   The diagnosis result display unit 113 displays information on the processing result of each functional unit for each mesh.

図2は、実施形態の流出解析装置1が用いる流出解析モデルの概念を説明する図である。図2は、隣接する2つのメッシュ(メッシュi及びメッシュj)が、それぞれのメッシュに存在するマンホール(マンホールi及びマンホールj)を介して下水管路Pに接続された様子を模式的に示した図である。図2に示すこれらのマンホールは仮想的なものであり、各マンホールの断面積Aは各メッシュに存在するマンホールの総断面積(すなわちマンホール総面積)で表される。同様に、図2に示す下水管路Pも仮想的なものであり、各メッシュに対応する部分の管路の容量は、各メッシュに存在する管路の総容量(すなわち、最大貯留量)で表される。 FIG. 2 is a diagram illustrating the concept of the outflow analysis model used by the outflow analysis device 1 of the embodiment. FIG. 2 schematically shows a state in which two adjacent meshes (mesh i and mesh j) are connected to the sewer pipe P via manholes (manhole i and manhole j) existing in the respective meshes. FIG. These manhole shown in FIG. 2 are virtual in that the cross-sectional area A h of the manhole is represented by the total cross-sectional area of the manhole present in each mesh (i.e. manhole total area). Similarly, the sewer pipe P shown in FIG. 2 is also virtual, and the capacity of the pipe corresponding to each mesh is the total capacity of the pipes existing in each mesh (that is, the maximum storage amount). expressed.

このような仮想的なマンホール及び下水管路を想定したモデルでは、仮想下水管が満管になった場合に、溢れた水が即座に地表面に流出するのではなく、仮想的なマンホールに一時的に貯留され、仮想的なマンホールが満管になった時点で地表面に流出する状況を模擬することができる。すなわち、実施形態の流出解析モデルでは、動水勾配や径深、流積を計算するために必要なパラメータである水位を、雨水の流出経路(下水管路及びマンホール)が考慮された値として算出することができる。このような流出解析モデルによれば、従来の仮想タンクモデルよりも実際の水理学的現象をより正確に模擬することが可能になる。   In the model that assumes such a virtual manhole and sewer pipe, when the virtual sewer pipe is full, the overflowed water does not immediately flow to the ground surface, but temporarily flows into the virtual manhole. It is possible to simulate the situation of being stored and being discharged to the ground surface when the virtual manhole is full. That is, in the runoff analysis model of the embodiment, the water level, which is a parameter necessary for calculating the hydrodynamic gradient, diameter depth, and flow product, is calculated as a value that considers the runoff path (sewage pipe and manhole) of rainwater. can do. According to such a runoff analysis model, an actual hydraulic phenomenon can be simulated more accurately than a conventional virtual tank model.

なお、図2では簡単のため2つのメッシュ間の接続関係を示しているが、実際には、あるメッシュの周りには上下左右に隣接する4つのメッシュと、斜め方向に隣接する4つのメッシュとの合計8つのメッシュが存在する。そのため、図2のようなモデルで実際に流出解析を行う場合には、あるメッシュについて図2と同様の接続関係を隣接する8つのメッシュとの間に適用して流出解析を行うことになる。   In FIG. 2, the connection relationship between two meshes is shown for the sake of simplicity, but actually, there are four meshes that are vertically and horizontally adjacent to one mesh, and four meshes that are diagonally adjacent to each other. There are a total of 8 meshes. Therefore, when the outflow analysis is actually performed using the model as shown in FIG. 2, the outflow analysis is performed by applying the same connection relation as that of FIG.

例えば、メッシュiとメッシュjとの流量収支と、仮想水位との関係は次の式(1)〜式(3)によって表される。   For example, the relationship between the flow rate balance between the mesh i and the mesh j and the virtual water level is expressed by the following equations (1) to (3).

Figure 2019139455
Figure 2019139455

Figure 2019139455
Figure 2019139455

Figure 2019139455
Figure 2019139455

式(1)〜式(3)において、H(i,j)はメッシュjに隣接するメッシュiの仮想水位を表す。S(i)はメッシュiにおける貯留量を、Slim(i)はメッシュiの上限貯留量を、Smax(i)はメッシュiの最大貯留量をそれぞれ表す。D(i,j)はメッシュiとメッシュjとを接続する仮想下水管の管径を表す。A(i)はメッシュiのマンホール総面積を、Aは各メッシュの面積(以下「メッシュ面積」という。)をそれぞれ表す。 In Expressions (1) to (3), H (i, j) represents the virtual water level of the mesh i adjacent to the mesh j. S (i) represents the storage amount in mesh i, S lim (i) represents the upper limit storage amount of mesh i, and S max (i) represents the maximum storage amount of mesh i. D (i, j) represents the pipe diameter of the virtual sewer pipe connecting the mesh i and the mesh j. A h (i) is a manhole total area of the mesh i, A m respectively represent the area of each mesh (hereinafter referred to as "mesh area".).

式(1)は、あるメッシュの貯留量が仮想下水管の容量を超えない場合における仮想水位の計算式である。式(1)は、あるメッシュiにおける仮想水位が、仮想下水管の管径D(i,j)と、メッシュiの下水管内に貯留しうる最大貯留容量Smax(i)に対する実際の貯留量の割合との積によって得られることを表している。すなわち、あるメッシュiにおける仮想水位は、仮想下水管の管径D(i,j)をメッシュiにおける仮想下水管の満管率で按分することによって得られる。この場合、メッシュiにおける貯留量が最大貯留容量Smax(i)を超過するまでは、仮想水位は0〜D(i,j)範囲の値をとる。そのため、D(i,j)として適切な管径が設定されれば、仮想水位を実際の下水管の水位の近似値として算出することができる。 Formula (1) is a formula for calculating the virtual water level when the storage amount of a certain mesh does not exceed the capacity of the virtual sewer pipe. Equation (1) indicates that the virtual water level in a certain mesh i is the actual storage amount relative to the pipe diameter D (i, j) of the virtual sewer pipe and the maximum storage capacity S max (i) that can be stored in the sewer pipe of the mesh i. It is obtained by the product with the ratio of. That is, the virtual water level in a certain mesh i is obtained by dividing the pipe diameter D (i, j) of the virtual sewer pipe by the fullness ratio of the virtual sewer pipe in the mesh i. In this case, the virtual water level takes a value in the range of 0 to D (i, j) until the storage amount in the mesh i exceeds the maximum storage capacity S max (i). Therefore, if an appropriate pipe diameter is set as D (i, j), the virtual water level can be calculated as an approximate value of the actual water level of the sewer pipe.

さらに、仮想下水管を円管ではなく矩形管として想定した場合には、仮想水位に対してより明確な物理的解釈が可能になる。すなわち、D(i,j)を矩形管の高さ方向(z軸方向)の辺の長さとした場合、矩形管の水位は式(1)で算出した仮想水位に一致する(勾配が0の場合)。式(1)はこの考え方を拡張したものである。従来の仮想タンクモデルでは、最大貯留容量及び貯留量に対して底面積というパラメータを用いて仮想水位を算出していたのに対し、本実施形態では、最大貯留容量及び貯留量に対して仮想下水管の管径(又は仮想下水管の高さ)という実際の管路情報に基づいて定義される仮想下水管のパラメータを用いて仮想水位を算出する。これにより、仮想水位の数値を現実的な値の範囲(0〜D(i,j))で得られるようにすることができる。   Furthermore, when the virtual sewer pipe is assumed to be a rectangular pipe instead of a circular pipe, a clearer physical interpretation is possible with respect to the virtual water level. That is, when D (i, j) is the length of the side in the height direction (z-axis direction) of the rectangular tube, the water level of the rectangular tube matches the virtual water level calculated by Equation (1) (the gradient is 0). If). Equation (1) is an extension of this idea. In the conventional virtual tank model, the virtual water level is calculated using the parameter of the bottom area with respect to the maximum storage capacity and the storage amount. The virtual water level is calculated using the virtual sewer parameters defined based on the actual pipe line information called the pipe diameter of the water pipe (or the height of the virtual sewer pipe). Thereby, the numerical value of the virtual water level can be obtained within a realistic value range (0 to D (i, j)).

式(2)は、あるメッシュの貯留量が仮想下水管の最大貯留量を超え、かつ上限貯留量を越えない場合の仮想水位の計算式を表している。貯留量が下水管の容量を超えた場合、下水管は満管状態であるため圧力管として取り扱うことができる。図2のモデルでは、仮想下水管を圧力管として考えた場合、仮想下水管の水頭(水位)はマンホールの高さとして得られる。式(2)はこのような場合の水頭(水位)を表現した式である。このような式を用いることによって、あるメッシュの仮想下水管が満管(圧力管状態)になってから、溢れた水が地表面に流出するまでの現象をより正確に表現することができる。   Formula (2) represents a formula for calculating the virtual water level when the storage amount of a certain mesh exceeds the maximum storage amount of the virtual sewer pipe and does not exceed the upper limit storage amount. When the amount of storage exceeds the capacity of the sewage pipe, the sewage pipe can be handled as a pressure pipe because it is full. In the model of FIG. 2, when the virtual sewer pipe is considered as a pressure pipe, the head (water level) of the virtual sewer pipe is obtained as the height of the manhole. Expression (2) is an expression expressing the water head (water level) in such a case. By using such an expression, it is possible to more accurately represent a phenomenon from when a virtual sewage pipe of a certain mesh becomes full (pressure pipe state) until overflowing water flows out to the ground surface.

式(3)は、あるメッシュの貯留量が仮想下水管の上限貯留量を超えた場合の仮想水位の計算式を表している。式(3)は、溢れた水が地表面に流出し、実際に浸水が発生した後の現象を表現する式である。雨水が地表面に流出した後は、地表面の面積に応じて微小に仮想水位が上昇すると考えられる。メッシュ面積Aは、マンホール総面積Aよりも格段に大きい値をとる。そのため、式(3)は、地表面流出後に水位が急上昇することがないことを表現している。 Formula (3) represents the calculation formula of the virtual water level when the storage amount of a certain mesh exceeds the upper limit storage amount of the virtual sewer pipe. Expression (3) is an expression that expresses a phenomenon after overflowing water flows out to the ground surface and water is actually generated. After rainwater flows out to the ground surface, the virtual water level is considered to rise slightly depending on the area of the ground surface. Mesh area A m takes a much greater value than manhole total area A h. Therefore, Formula (3) expresses that the water level does not rise rapidly after the ground surface runoff.

以下、式(1)〜式(3)を用いて、各メッシュ間の水量の収支を計算する方法を説明する。本実施形態における水量の収支の計算方法は、基本的には従来の仮想タンクモデルと同様である。具体的には、各メッシュ間の水量の収支は、次の式(4)によって表される「連続の式」と、式(5)によって表される「マニング式」とを用いて算出することができる。   Hereinafter, a method for calculating the balance of the amount of water between the meshes will be described using Expressions (1) to (3). The calculation method of the water amount balance in this embodiment is basically the same as the conventional virtual tank model. Specifically, the balance of the amount of water between each mesh should be calculated using the “continuous formula” expressed by the following formula (4) and the “Manning formula” expressed by the formula (5). Can do.

Figure 2019139455
Figure 2019139455

Figure 2019139455
Figure 2019139455

式(4)及び式(5)において、Sはあるメッシュの貯留量を表す。貯留量Sについて、ここではiやj等のメッシュの識別子を省略して記載している。これは、以下で説明する他のパラメータについても同様である。Rはあるメッシュの有効降雨量を表す。Qはあるメッシュから隣接メッシュへの流出量を表す。Aは仮想下水管において、水が占める部分の断面積(流積)を表す。Rは仮想下水管の径深を表す。Iは仮想下水管の動水勾配を表す。nは仮想下水管の粗度を表す。ここで、式(5)を式(4)に代入することにより次の式(6)が得られる。式(6)において、流積A、径深R、動水勾配Iは、仮想下水管の水位に依存するパラメータである。 In Formula (4) and Formula (5), S represents the storage amount of a certain mesh. Here, the stored amount S is described by omitting mesh identifiers such as i and j. The same applies to other parameters described below. R e represents effective rainfall Mesh. Q represents an outflow amount from a certain mesh to an adjacent mesh. A d is in the virtual sewers, representing the cross-sectional area of a portion which water occupies (Nagareseki). R represents the diameter of the virtual sewer pipe. I d represents the hydraulic gradient of virtual sewer. n represents the roughness of the virtual sewer pipe. Here, the following equation (6) is obtained by substituting equation (5) into equation (4). In equation (6), the flow product A d , the radial depth R, and the dynamic water gradient I d are parameters that depend on the water level of the virtual sewer pipe.

Figure 2019139455
Figure 2019139455

図3は、仮想水位と、流積及び径深との関係を示す図である。ここで、仮想下水管を辺長Lの正方形の断面を持つ矩形管と想定した場合(図3(A))、流積A及び径深Rは仮想下水管の水位Hを用いて次の式(7)及び式(8)のように表される。また、仮想下水管を直径Dの円管と想定した場合(図3(B))、流積A及び径深Rは仮想下水管の水位Hを用いて次の式(9)〜式(11)のように表される。式(7)〜式(11)の各パラメータは、図3の各パラメータに対応している。 FIG. 3 is a diagram illustrating the relationship between the virtual water level, the flow volume, and the diameter depth. Here, assuming a rectangular tube having a square cross section of side length L virtual sewage pipe (FIG. 3 (A)), the Nagareseki A d and hydraulic radius R follows using the water level H of the virtual sewer It represents like Formula (7) and Formula (8). Further, when the virtual sewer pipe assumes a circular tube having a diameter D (FIG. 3 (B)), Nagareseki A d and hydraulic radius R by using the water level H of the virtual sewer following equation (9) to ( 11). Each parameter of Formula (7)-Formula (11) respond | corresponds to each parameter of FIG.

Figure 2019139455
Figure 2019139455

Figure 2019139455
Figure 2019139455

Figure 2019139455
Figure 2019139455

Figure 2019139455
Figure 2019139455

Figure 2019139455
Figure 2019139455

図4は、仮想水位と、動水勾配及び管路勾配との関係を示す図である。管路勾配は、仮想下水管自身の勾配を表し、動水勾配は仮想下水管を流下する水面の勾配を表す。図4に示す関係から、動水勾配Iは、仮想水位H及び管路勾配Iを用いて次の式(12)のように近似することができる。式(12)の各パラメータは、図4の各パラメータに対応している。 FIG. 4 is a diagram illustrating the relationship between the virtual water level, the dynamic water gradient, and the pipeline gradient. The pipe gradient represents the gradient of the virtual sewer pipe itself, and the dynamic water gradient represents the gradient of the water surface flowing down the virtual sewer pipe. From the relationship shown in FIG. 4, hydraulic gradient I d can be approximated as the following equation (12) using a virtual water level H and line slope I s. Each parameter of Expression (12) corresponds to each parameter in FIG.

Figure 2019139455
Figure 2019139455

式(12)において、Kはメッシュの辺長を表し、メッシュの辺長はすなわち隣接するメッシュ間の距離を表す。H及びH’は、あるメッシュでの仮想下水管の水位(H)と隣接するメッシュでの仮想下水管の水位(H’)を表す。H及びH’は、上記の式(1)〜式(3)によって算出される。なお、管路勾配I及び動水勾配Iについて、IがIに対して十分に小さいと仮定できる場合にはI≒0としても良い。 In Expression (12), K represents the side length of the mesh, and the side length of the mesh represents the distance between adjacent meshes. H and H ′ represent the water level (H ′) of the virtual sewer pipe in a mesh adjacent to the water level (H ′) of the virtual sewer pipe in a certain mesh. H and H ′ are calculated by the above formulas (1) to (3). Note that the line gradient I s and hydraulic gradient I d, if it can be assumed that I s is sufficiently small relative to the I d may be a I s ≒ 0.

以上説明した式(6)〜式(12)を組み合わせることにより流出解析モデルが得られる。ここでは、例として仮想下水管を円管と想定した場合の流出解析モデルについて説明する。円管の場合、式(9)〜式(12)が基本式となり、これらをまとめると次の式(13)で表される流出解析モデルが得られる。具体的には、式(13)は、式(9)〜式(12)を式(6)に代入して整理することによって得られる。   An outflow analysis model is obtained by combining the equations (6) to (12) described above. Here, as an example, an outflow analysis model when the virtual sewer pipe is assumed to be a circular pipe will be described. In the case of a circular pipe, equations (9) to (12) become basic equations, and when these are combined, an outflow analysis model represented by the following equation (13) is obtained. Specifically, Expression (13) is obtained by substituting Expression (9) to Expression (12) into Expression (6) and rearranging.

Figure 2019139455
Figure 2019139455

ここで、式(13)におけるθは上記の式(11)によって表される。式(13)の右辺第2項は、隣接するメッシュ間の流入出量を表している。水量収支の関係上、あるメッシュから流出(又は流入)する水の流量と、隣接するメッシュに流入(又は流出)する水の流量とは、絶対値が同じで、かつ符号が異なる量となる。そのため、流出解析の精度を高めるためには、各メッシュ間で水深角θに関する平均化処理などを行うことによって、各メッシュ間の水量収支が上記関係を満たすように調整されてもよい。   Here, θ in the equation (13) is expressed by the above equation (11). The second term on the right side of Equation (13) represents the inflow / outflow amount between adjacent meshes. In terms of the water balance, the flow rate of water flowing out (or flowing in) from a mesh and the flow rate of water flowing (or flowing out) into an adjacent mesh have the same absolute value but different signs. Therefore, in order to improve the accuracy of the outflow analysis, the water balance between the meshes may be adjusted so as to satisfy the above relationship by performing an averaging process on the water depth angle θ between the meshes.

なお、本実施形態では、流出解析モデルの構築にマニング式を用いているが、マニング式に代えて動水勾配を変数として含む他の式を用いても良い。例えば、クッター式を用いて、式(11)及び式(13)に相当する式を導出してもよい。   In this embodiment, the Manning equation is used to construct the outflow analysis model, but another equation including a hydraulic gradient as a variable may be used instead of the Manning equation. For example, equations corresponding to the equations (11) and (13) may be derived using the Kutta equation.

また、上記の実施形態では、仮想下水管として円管を想定した場合の流出解析モデルについて説明したが、これに代えて、仮想下水管として矩形管を想定した場合の流出解析モデルを導出してもよい。さらに、モデルの簡略化のため、径深や流積を仮想水位によって変動するパラメータとせずに固定値としてもよい。例えば径深=D/4、流積=π×D/8としても良い。 In the above embodiment, the outflow analysis model when a circular pipe is assumed as a virtual sewer pipe has been described. Instead, an outflow analysis model when a rectangular pipe is assumed as a virtual sewer pipe is derived. Also good. Furthermore, for simplification of the model, the diameter depth and the flow product may be fixed values without using parameters that vary depending on the virtual water level. For example, it is good also as diameter depth = D / 4 and flow product = (pi) * D < 2 > / 8.

すなわち、動水勾配が仮想水位に応じて変化するパラメータであること及び、仮想水位が式(1)〜式(3)のように貯留量に応じた算出式で求められることが表現されたモデルであれば、流出解析モデルには式(11)及び式(13)と異なる式が用いられても良い。   That is, a model that expresses that the dynamic water gradient is a parameter that changes according to the virtual water level, and that the virtual water level is obtained by a calculation formula according to the storage amount as in equations (1) to (3). If so, an expression different from Expression (11) and Expression (13) may be used for the outflow analysis model.

なお、上記説明では簡単のため、あるメッシュと、そのメッシュに隣接する1つのメッシュとの間の流出解析を行う場合を例に説明したが、実際の解析を行う場合には、あるメッシュに対して、隣接する全てのメッシュの影響を考慮する必要がある。具体的には、式(13)の第2項として、隣接する全てのメッシュ(例えば、上下左右斜めの全8メッシュ)について考慮した項を付加する必要がある。このように構築された流出解析モデルに対し、各メッシュの有効降雨量を入力として流量収支を計算することによって各メッシュの流出量(又は流入量)が算出される。   In the above description, for the sake of simplicity, an example has been described in which an outflow analysis between a mesh and one mesh adjacent to the mesh is performed. However, when an actual analysis is performed, Therefore, it is necessary to consider the influence of all adjacent meshes. Specifically, it is necessary to add a term that considers all adjacent meshes (for example, all eight meshes that are diagonally up, down, left, and right) as the second term of Expression (13). For the runoff analysis model constructed as described above, the runoff amount (or inflow amount) of each mesh is calculated by calculating the flow rate balance with the effective rainfall amount of each mesh as an input.

なお、仮想下水管の管径、仮想下水管の管路勾配、仮想下水管の粗度、マンホール総面積、最大貯留量、上限貯留量などのパラメータは、下水道管理台帳データを用いて算出可能である。これらのパラメータは、仮想下水管径算出部108、仮想下水管勾配算出部109、仮想下水管粗度算出部110、マンホール総面積算出部111、容量算出部112等によって、各メッシュ、又は互いに隣接するメッシュごとに予め算出される。   Parameters such as virtual sewage pipe diameter, virtual sewage pipe slope, virtual sewage pipe roughness, total manhole area, maximum storage volume, and maximum storage volume can be calculated using sewer management ledger data. is there. These parameters are determined by the virtual sewage pipe diameter calculation unit 108, the virtual sewage pipe gradient calculation unit 109, the virtual sewage pipe roughness calculation unit 110, the manhole total area calculation unit 111, the capacity calculation unit 112, and the like, or are adjacent to each other. It is calculated in advance for each mesh.

続いて、浸水リスクの評価について説明する。浸水リスク評価部104は、メッシュごとに、各メッシュについて算出された貯留量に基づいて以下のような評価を行う。例えば、あるメッシュの貯留量が最大貯留容量Smaxを超過した場合、浸水リスク評価部104は、そのメッシュについて「浸水の可能性がある」と判断する。また、あるメッシュの貯留量が上限貯留容量Slimを超過した場合、浸水リスク評価部104は、そのメッシュについて「浸水の可能性が極めて高い」と判断する。このような浸水リスクの判断基準は、診断対象地域の特性等に応じて適宜設定されてよい。 Next, the evaluation of flood risk will be described. The inundation risk evaluation unit 104 performs the following evaluation for each mesh based on the storage amount calculated for each mesh. For example, when the storage amount of a certain mesh exceeds the maximum storage capacity Smax , the inundation risk evaluation unit 104 determines that “there is a possibility of inundation” for the mesh. Further, when the storage amount of a certain mesh exceeds the upper limit storage capacity Slim , the inundation risk evaluation unit 104 determines that “the possibility of inundation is extremely high” for the mesh. Such criteria for judging inundation risk may be appropriately set according to the characteristics of the diagnosis target area.

診断結果表示部113は、このように診断された各メッシュの浸水リスク(流出解析の結果を含む)を視覚的に識別可能な態様で表示させる。例えば、診断結果表示部113は、浸水リスクを各メッシュに対応させたメッシュ状の分布として表示させる。具体的には、貯留量が最大貯留量を超過したメッシュを黄色で表示させ、上限貯留量を超過したメッシュを赤色で表示させても良い。また、最大貯留量又は上限貯留量と、算出された貯留量に基づいて仮想下水管の満管率を算出し、満管率に応じた色で各メッシュを表示させるようにしてもよい。診断結果表示部113は、流出解析の結果や診断リスクの診断結果等に基づいて得られる情報であれば、どのような情報を表示させてもよい。また、情報表示の態様も、上記のような色で識別可能にする態様のほか、文字や記号、図形、グラフ等の視覚的な識別を可能にする態様であれば、他のどのような態様であってもよい。   The diagnosis result display unit 113 displays the inundation risk (including the result of the outflow analysis) of each mesh diagnosed in this manner in a visually identifiable manner. For example, the diagnosis result display unit 113 displays the inundation risk as a mesh distribution corresponding to each mesh. Specifically, a mesh whose stored amount exceeds the maximum stored amount may be displayed in yellow, and a mesh whose stored amount exceeds the upper limit stored amount may be displayed in red. Further, the full fill rate of the virtual sewer pipe may be calculated based on the maximum storage amount or the upper limit storage amount and the calculated storage amount, and each mesh may be displayed in a color corresponding to the full pipe rate. The diagnosis result display unit 113 may display any information as long as the information is obtained based on the result of the outflow analysis, the diagnosis result of the diagnosis risk, and the like. Also, the information display mode is not limited to the above-described color identification mode, but any other mode as long as it allows visual identification of characters, symbols, figures, graphs, etc. It may be.

このように構成された実施形態の流出解析装置1は、従来の仮想タンクモデルを拡張し、仮想タンクを、仮想下水管と、仮想下水管と地表面とを接続する仮想マンホールとで表した流出解析モデルを用いて流出解析を行う構成を備える。このような構成を備えることにより、簡易な解析モデルを用いつつ、より現実に即した流出解析結果を得ることができる。   The runoff analysis device 1 of the embodiment configured as described above extends the conventional virtual tank model, and the virtual tank is represented by a virtual sewer pipe and a virtual manhole connecting the virtual sewer pipe and the ground surface. It has a configuration that performs runoff analysis using an analysis model. By providing such a configuration, it is possible to obtain a more realistic outflow analysis result while using a simple analysis model.

また、実施形態の流出解析装置1は、容易に取得可能な下水管路の物理情報を組み合わせて仮想タンクモデルを拡張した流出解析モデルで流出解析を行うため、分布型流出解析のような複雑な土木情報を不要としつつ、より精度の高い流出解析を行うことができる。具体的には、電子的に保存された下水道管理台帳データと、好ましくはGISデータとの情報を入手するだけで、分布型流出解析の解析結果や、解析結果から想定される実際の流出現象や浸水現象を、降雨レーダのメッシュサイズの解像度で近似することができる。   In addition, since the runoff analysis apparatus 1 of the embodiment performs the runoff analysis with the runoff analysis model obtained by combining the physical information of the sewer pipes that can be easily acquired and extending the virtual tank model, the runoff analysis apparatus 1 is complicated as in the distributed runoff analysis. It is possible to perform more accurate outflow analysis while making civil engineering information unnecessary. Specifically, just by obtaining information on electronically stored sewer management ledger data, and preferably GIS data, the analysis results of distributed runoff analysis, actual runoff phenomena assumed from the analysis results, The inundation phenomenon can be approximated with the resolution of the mesh size of the precipitation radar.

また、実施形態の流出解析装置1は、流出解析の結果や浸水リスクの診断結果を、メッシュごとにかつ視覚的に識別可能な態様で表示させることができる。そのため、利用者は、診断対象地域に関する総合的な浸水リスクに関する情報を容易に把握することができる。   Moreover, the outflow analysis apparatus 1 of the embodiment can display the result of the outflow analysis and the diagnosis result of the inundation risk for each mesh in a visually identifiable manner. Therefore, the user can easily grasp information on the overall inundation risk related to the diagnosis target area.

また、実施形態の流出解析装置1は、分布型流出解析のような複雑な水理モデルを用いた流出解析を行うのではなく、メッシュの個数と同数の1次の微分方程式を解くだけで流出解析を行える仮想タンクモデルを拡張したモデルに基づいて流出解析を行う。そのため、一周期の計算をより短い時間で行うことができ、リアルタイム性を向上させることができる。その結果、流出解析装置1は、流出解析や浸水リスクの診断の結果をリアルタイムハザードマップやリアルタイムで運用される雨水排水設備や雨水貯留設備運用の支援情報として提供することができる。   In addition, the outflow analysis device 1 of the embodiment does not perform outflow analysis using a complicated hydraulic model such as distributed outflow analysis, but only solves the same number of primary differential equations as the number of meshes. Runoff analysis is performed based on an extended model of a virtual tank model that can perform analysis. Therefore, one cycle can be calculated in a shorter time, and real-time performance can be improved. As a result, the runoff analysis apparatus 1 can provide the results of runoff analysis and inundation risk diagnosis as real-time hazard maps and support information for rainwater drainage equipment and rainwater storage equipment operation that are operated in real time.

以上説明した少なくともひとつの実施形態によれば、貯留量に応じて各メッシュの仮想水位を算出する流出解析部を持つことにより、簡易な解析モデルを用いつつ、より現実に即した流出解析結果を得ることができる。   According to at least one embodiment described above, by having an outflow analysis unit that calculates the virtual water level of each mesh according to the storage amount, a more realistic outflow analysis result can be obtained using a simple analysis model. Can be obtained.

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。   Although several embodiments of the present invention have been described, these embodiments are presented by way of example and are not intended to limit the scope of the invention. These embodiments can be implemented in various other forms, and various omissions, replacements, and changes can be made without departing from the spirit of the invention. These embodiments and their modifications are included in the scope and gist of the invention, and are also included in the invention described in the claims and the equivalents thereof.

1…流出解析装置、101…降雨分布入力部、102…流出係数算出部、103…流出解析部、104…浸水リスク評価部、105…地図情報記憶部、106…メッシュ面積設定部、107…下水道管理台帳データ記憶部、108…仮想下水管径算出部、109…仮想下水管勾配算出部、110…仮想下水管粗度算出部、111…マンホール総面積算出部、112…容量算出部、113…診断結果表示部 DESCRIPTION OF SYMBOLS 1 ... Runoff analysis apparatus, 101 ... Rainfall distribution input part, 102 ... Runoff coefficient calculation part, 103 ... Runoff analysis part, 104 ... Inundation risk evaluation part, 105 ... Map information storage part, 106 ... Mesh area setting part, 107 ... Sewer Management ledger data storage unit 108 ... virtual sewage pipe diameter calculation unit 109 ... virtual sewage pipe gradient calculation unit 110 ... virtual sewage pipe roughness calculation unit 111 ... manhole total area calculation unit 112 ... capacity calculation unit 113 ... Diagnostic result display

Claims (3)

複数のメッシュに分割された診断対象地域の降雨量を示す情報であって前記メッシュごとの降雨量を示す降雨分布情報を入力する降雨分布入力部と、
前記降雨分布情報に基づいて各メッシュ間における流量の収支バランスを計算することにより各メッシュの貯留量を算出する流出解析部と、
を備え、
前記流出解析部は、算出された各メッシュの貯留量と、前記メッシュごとに設定された仮想管路の最大貯留量及び前記仮想管路と地表面とを接続する仮想マンホールの最大貯留量とに基づいて、前記貯留量が前記仮想管路の最大貯留量未満であるメッシュの仮想水位を以下の式(1)で算出し、
仮想水位=貯留量÷最大貯留量×仮想管路直径 式(1)
前記貯留量が前記仮想管路の最大貯留量以上であり、かつ前記仮想管路の最大貯留量と前記仮想マンホールの最大貯留量との和である上限貯留量未満であるメッシュの仮想水位を以下の式(2)で算出し、
仮想水位=仮想管路直径+(貯留量−仮想管路の最大貯留量)÷仮想マンホールの断面積 式(2)
前記貯留量が上限貯留量以上であるメッシュの仮想水位を以下の式(3)で算出し、
仮想水位=仮想管路直径+(上限貯留量−最大貯留量)÷仮想マンホールの断面積+(貯留量−上限貯留量)÷メッシュの面積 式(3)
前記貯留量に応じて算出された前記仮想水位を用いて各メッシュの浸水リスクの評価に必要な諸量を算出する、
流出解析装置。
A rainfall distribution input unit that inputs rainfall distribution information indicating the rainfall amount of each mesh that is the rainfall amount of the diagnosis target area divided into a plurality of meshes;
An outflow analysis unit for calculating a storage amount of each mesh by calculating a balance of flow rate between each mesh based on the rainfall distribution information;
With
The outflow analysis unit calculates the storage amount of each mesh, the maximum storage amount of the virtual pipe set for each mesh, and the maximum storage amount of the virtual manhole connecting the virtual pipe and the ground surface. Based on the virtual water level of the mesh that the storage amount is less than the maximum storage amount of the virtual pipeline, according to the following formula (1),
Virtual water level = storage volume ÷ maximum storage volume × virtual pipe diameter formula (1)
The virtual water level of the mesh whose storage amount is equal to or greater than the maximum storage amount of the virtual pipeline and less than the upper limit storage amount, which is the sum of the maximum storage amount of the virtual pipeline and the maximum storage amount of the virtual manhole. Calculated by the equation (2)
Virtual water level = Virtual pipeline diameter + (Storage volume-Maximum storage volume of virtual pipeline) / Virtual manhole cross-sectional area formula (2)
The virtual water level of the mesh where the storage amount is equal to or greater than the upper limit storage amount is calculated by the following equation (3),
Virtual water level = virtual pipeline diameter + (upper limit storage amount-maximum storage amount) / virtual manhole cross-sectional area + (storage amount-upper limit storage amount) / mesh area Formula (3)
Using the virtual water level calculated according to the storage amount, various amounts necessary for evaluating the inundation risk of each mesh are calculated.
Runoff analyzer.
前記流出解析部は、各メッシュ間の動水勾配、各メッシュの仮想管路の径深又は流積の少なくとも1つを前記諸量として算出し、
前記諸量に基づいて、各メッシュの浸水リスクを評価する浸水リスク評価部をさらに備える、
請求項1に記載の流出解析装置。
The outflow analysis unit calculates, as the quantities, the hydrodynamic gradient between the meshes, the diameter of the virtual pipe line of each mesh, or the flow product,
Further comprising a flood risk evaluation unit for evaluating the flood risk of each mesh based on the various quantities,
The outflow analysis device according to claim 1.
複数のメッシュに分割された診断対象地域の降雨量を示す情報であって前記メッシュごとの降雨量を示す降雨分布情報を入力する降雨分布入力ステップと、
前記降雨分布情報に基づいて各メッシュ間における流量の収支バランスを計算することにより各メッシュの貯留量を算出する流出解析ステップと、
を備え、
前記流出解析ステップでは、算出された各メッシュの貯留量と、前記メッシュごとに設定された仮想管路の最大貯留量及び前記仮想管路と地表面とを接続する仮想マンホールの最大貯留量とに基づいて、前記貯留量が前記仮想管路の最大貯留量未満であるメッシュの仮想水位を以下の式(4)で算出し、
仮想水位=貯留量÷最大貯留量×仮想管路直径 式(4)
前記貯留量が前記仮想管路の最大貯留量以上であり、かつ前記仮想管路の最大貯留量と前記仮想マンホールの最大貯留量との和である上限貯留量未満であるメッシュの仮想水位を以下の式(5)で算出し、
仮想水位=仮想管路直径+(貯留量−仮想管路の最大貯留量)÷仮想マンホールの断面積 式(5)
前記貯留量が上限貯留量以上であるメッシュの仮想水位を以下の式(6)で算出し、
仮想水位=仮想管路直径+(上限貯留量−最大貯留量)÷仮想マンホールの断面積+(貯留量−上限貯留量)÷メッシュの面積 式(6)
前記貯留量に応じて算出された前記仮想水位を用いて各メッシュの浸水リスクの評価に必要な諸量を算出する、
流出解析方法。
A rainfall distribution input step for inputting rainfall distribution information indicating the rainfall amount of each of the meshes, which is information indicating the rainfall amount of the diagnosis target area divided into a plurality of meshes;
An outflow analysis step for calculating a storage amount of each mesh by calculating a balance of flow rate between each mesh based on the rainfall distribution information;
With
In the outflow analysis step, the calculated storage amount of each mesh, the maximum storage amount of the virtual conduit set for each mesh, and the maximum storage amount of the virtual manhole connecting the virtual conduit and the ground surface Based on the virtual water level of the mesh that the storage amount is less than the maximum storage amount of the virtual pipeline based on the following formula (4),
Virtual water level = storage volume ÷ maximum storage volume × virtual pipeline diameter formula (4)
The virtual water level of the mesh whose storage amount is equal to or greater than the maximum storage amount of the virtual pipeline and less than the upper limit storage amount, which is the sum of the maximum storage amount of the virtual pipeline and the maximum storage amount of the virtual manhole. (5)
Virtual water level = Virtual pipeline diameter + (Storage volume-Maximum storage volume of virtual pipeline) / Virtual manhole cross-sectional area formula (5)
The virtual water level of the mesh in which the storage amount is equal to or greater than the upper limit storage amount is calculated by the following equation (6),
Virtual water level = virtual pipeline diameter + (upper limit storage amount-maximum storage amount) / virtual manhole cross-sectional area + (storage amount-upper limit storage amount) / mesh area Formula (6)
Using the virtual water level calculated according to the storage amount, various amounts necessary for evaluating the inundation risk of each mesh are calculated.
Runoff analysis method.
JP2018021427A 2018-02-08 2018-02-08 Spill analyzer, spill analysis method, and computer program Active JP7013270B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018021427A JP7013270B2 (en) 2018-02-08 2018-02-08 Spill analyzer, spill analysis method, and computer program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018021427A JP7013270B2 (en) 2018-02-08 2018-02-08 Spill analyzer, spill analysis method, and computer program

Publications (2)

Publication Number Publication Date
JP2019139455A true JP2019139455A (en) 2019-08-22
JP7013270B2 JP7013270B2 (en) 2022-01-31

Family

ID=67694026

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018021427A Active JP7013270B2 (en) 2018-02-08 2018-02-08 Spill analyzer, spill analysis method, and computer program

Country Status (1)

Country Link
JP (1) JP7013270B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021166480A1 (en) 2020-02-18 2021-08-26 コベルコ建機株式会社 Disaster countermeasure support server, disaster countermeasure support system, and disaster countermeasure support method
JP7220276B1 (en) 2021-12-27 2023-02-09 アドソル日進株式会社 Inundation simulation device, inundation simulation method and program

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0836681A (en) * 1994-07-25 1996-02-06 Nohmi Bosai Ltd Fire detection equipment using optical fiber
JP2017194344A (en) * 2016-04-20 2017-10-26 株式会社東芝 Flooding risk diagnosis apparatus, flooding risk diagnostic method, controller, and computer program

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0836681A (en) * 1994-07-25 1996-02-06 Nohmi Bosai Ltd Fire detection equipment using optical fiber
JP2017194344A (en) * 2016-04-20 2017-10-26 株式会社東芝 Flooding risk diagnosis apparatus, flooding risk diagnostic method, controller, and computer program

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021166480A1 (en) 2020-02-18 2021-08-26 コベルコ建機株式会社 Disaster countermeasure support server, disaster countermeasure support system, and disaster countermeasure support method
JP7220276B1 (en) 2021-12-27 2023-02-09 アドソル日進株式会社 Inundation simulation device, inundation simulation method and program
JP2023096271A (en) * 2021-12-27 2023-07-07 アドソル日進株式会社 Inundation simulator, inundation simulation method, and program

Also Published As

Publication number Publication date
JP7013270B2 (en) 2022-01-31

Similar Documents

Publication Publication Date Title
JP7055654B2 (en) Outflow analysis device and outflow analysis parameter adjustment method
JP6207889B2 (en) Inundation prediction system, inundation prediction method and program
CN104851360B (en) The generation method and system of a kind of flood risk mapping
JP6716328B2 (en) Inundation risk diagnostic device, inundation risk diagnostic method, control device, and computer program
CN114997541A (en) Urban inland inundation prediction method and early warning platform based on digital twin technology
Kamel Application of a hydrodynamic MIKE 11 model for the Euphrates River in Iraq
JP7013270B2 (en) Spill analyzer, spill analysis method, and computer program
Parker The Control of Water: As Applied to Irrigation, Power and Town Water Supply Purposes
Azarang et al. Simulating the erosion and sedimentation of Karun Alluvial River in the Region of Ahvaz (southwest of Iran)
ShahiriParsa et al. Introduction to floodplain zoning simulation models through dimensional approach
Chen et al. Numerical simulation of potential inundation in a coastal zone
YANAGISAKI et al. Debris flow simulation by applying the Hyper KANAKO system for water and sediment runoff from overtopping erosion of a landslide dam
Ganjidoost et al. Shifting to a monetized quantitative approach for risk analysis using property damages
Pion et al. Dredging volumes prediction for the access channel of santos port considering different design depths
Ramsbottom et al. Developing a national programme of flood risk management measures: Moldova
Nardini et al. MODCEL vs. IBER: a comparison of flooding models in Riohacha, a coastal town of La Guajira, Colombia
Laz Morphological assessment of a selected reach of Jamuna river by using DELFT3D model
Milišić et al. Estimation of channel and flood plain roughness using HEC-RAS model: A case study of the Veseočica River, Bosnia and Herzegovina
Mustafa et al. Effects of tides on the hydrology and geometry of a freshwater channel
Hérault The impacts of varying 2D modelling strategies on flood hazard assessment in urban areas: case study of the Garonne River flood risk prevention plan
Tuukkanen et al. A GIS-based model for ditch erosion risk assessment in peatland forestry
Brabets et al. Geomorphology and river dynamics of the lower Copper River, Alaska
Kantarzhi Forecast of the Dynamics of a Sandy Beach in Complexed Hydrodynamic Conditions
Xia et al. High-Performance Integrated hydrodynamic Modelling of Storm Induced Floods at a Catchment Scale
Lulla Hydrodynamic optimization of wetland restoration measures at the Storelva river floodplain with its meanders and oxbow lakes in Southern Norway

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20201207

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210929

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20211005

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211206

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20211221

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220119

R150 Certificate of patent or registration of utility model

Ref document number: 7013270

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150