JP2016130435A - Method, program and system of determining sensing point in sewerage system - Google Patents

Method, program and system of determining sensing point in sewerage system Download PDF

Info

Publication number
JP2016130435A
JP2016130435A JP2015005385A JP2015005385A JP2016130435A JP 2016130435 A JP2016130435 A JP 2016130435A JP 2015005385 A JP2015005385 A JP 2015005385A JP 2015005385 A JP2015005385 A JP 2015005385A JP 2016130435 A JP2016130435 A JP 2016130435A
Authority
JP
Japan
Prior art keywords
location
time
series data
sewer network
sensing
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
JP2015005385A
Other languages
Japanese (ja)
Other versions
JP6503745B2 (en
Inventor
祐輔 樋田
Yusuke Toida
祐輔 樋田
鈴木 貴志
Takashi Suzuki
貴志 鈴木
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.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2015005385A priority Critical patent/JP6503745B2/en
Publication of JP2016130435A publication Critical patent/JP2016130435A/en
Application granted granted Critical
Publication of JP6503745B2 publication Critical patent/JP6503745B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

PROBLEM TO BE SOLVED: To provide a method and a system for determining proper sensing points in a sewerage system.SOLUTION: The method for determining proper sensing points in a sewerage system includes the steps of: calculating time series data of water level at plural points (manhole) in a sewerage system on the basis of various kinds of databases (S12); calculating a delay time τ at a specific point (lower most manhole) in the sewerage system with respect to other points of manholes in the time series data (S14, S16); calculating significance probability at each of the other points of manholes by means of testing of statistical hypothesis by using calculated delay time τ and the time series data at the other points of manholes (S18, S20); and determining optimum installation points (sensing point) of sensors on the basis of the calculated significance probability (S22 to S30).SELECTED DRAWING: Figure 6

Description

本発明は、下水道網のセンシング箇所決定方法、プログラム及び装置に関する。   The present invention relates to a sewer network sensing location determination method, program, and apparatus.

都市下水道は、降雨時に下水管路に流入してくる雨水を河川に排出したり、生活排水等の汚水を処理することを目的としており、下水道網の下流域にはポンプ場や下水処理場が設置されている。これらの下水道設備は降雨量や流入量等を加味して設計されているため、従来においては下水道管渠から溢水することはほとんどなかった。   Urban sewers are intended to discharge rainwater that flows into sewer pipes during rainfall into rivers and to treat sewage such as domestic wastewater. Pumps and sewage treatment plants are located downstream of the sewer network. is set up. Since these sewerage facilities are designed taking into account the rainfall and inflow, etc., there has been almost no overflow from sewerage pipes in the past.

しかし、近年においては計画降雨量を超える集中豪雨等による、都市での浸水被害が多発しているため、事前に浸水対策を講じる必要が出てきている。浸水対策の1つとして、例えば、下水道管の水位、流量、流速等の下水状況データをセンサにてリアルタイムで計測して下水道管路内の状況を把握し、ポンプ場や下水処理場での効率的な排水運用を行うことが知られている(例えば、特許文献1等参照)。   However, in recent years, inundation damage has frequently occurred in cities due to torrential rain that exceeds the planned rainfall, and it is necessary to take measures against inundation in advance. As one of the countermeasures against inundation, for example, sewage status data such as sewer pipe water level, flow rate, and flow velocity is measured in real time by a sensor to grasp the situation in the sewer pipe, and efficiency at the pump station and sewage treatment plant. It is known to perform a general drainage operation (see, for example, Patent Document 1).

特開2010−203964号公報JP 2010-203964 A 特開2006−18763号公報JP 2006-18763 A

一般的に、センサを多数設置すると計測データ量が膨大になるとともに、設置費用が高額になるため、センサ設置数はできるだけ少ないことが望ましい。しかしながら、下水道網の状況を正確に把握するためには、センサ設置数を無闇に減じることはできない。   Generally, when a large number of sensors are installed, the amount of measurement data becomes enormous and the installation cost becomes high. Therefore, it is desirable that the number of sensors installed be as small as possible. However, in order to accurately grasp the state of the sewer network, the number of sensors installed cannot be reduced without any darkness.

なお、特許文献2等には、予め設定した遅れ時間を利用して下水処理場に流入する下水量を予測する技術が開示されているが、本技術は、下水流入量を予測するためのものであって、センサの設置箇所や数を決定するためのものではない。   Patent Document 2 and the like disclose a technique for predicting the amount of sewage flowing into a sewage treatment plant using a preset delay time, but this technique is for predicting the amount of sewage inflow. However, it is not for determining the location and number of sensors.

1つの側面では、本発明は、下水道網において適切なセンシング箇所を決定することが可能な下水道網のセンシング箇所決定方法、プログラム及び装置を提供することを目的とする。   In one aspect, an object of the present invention is to provide a sewer network sensing location determination method, program, and apparatus capable of determining an appropriate sensing location in a sewer network.

一つの態様では、下水道網のセンシング箇所決定方法は、下水道網に関するデータと降雨量に関するデータに基づき、前記下水道網の複数の箇所における水位、流量、流速の少なくとも1つの時系列データを算出し、前記下水道網の特定箇所における前記時系列データの、該特定箇所以外の箇所における前記時系列データそれぞれに対する遅れ時間を算出し、算出した前記遅れ時間と、前記特定箇所以外の箇所それぞれにおける前記時系列データとを用いて、統計的仮説検定により前記特定箇所以外の箇所それぞれの有意確率を算出し、算出された前記有意確率に基づいてセンサを配置するセンシング箇所を決定する、処理をコンピュータが実行する下水道網のセンシング箇所決定方法である。   In one aspect, a method for determining a sensing location of a sewer network calculates at least one time series data of a water level, a flow rate, and a flow velocity at a plurality of locations of the sewer network based on data on the sewer network and data on rainfall. The time series data at a specific location of the sewer network is calculated with respect to each time series data at a location other than the specific location, and the calculated delay time and the time series at each location other than the specific location are calculated. The computer executes a process of calculating a significance probability of each location other than the specific location by statistical hypothesis testing using the data and determining a sensing location where the sensor is arranged based on the calculated significance probability This is a method for determining the sensing location of the sewer network.

下水道網において適切なセンシング箇所を決定することができる。   Appropriate sensing points can be determined in the sewer network.

一実施形態に係るセンシング箇所決定装置のハードウェア構成を概略的に示す図である。It is a figure which shows roughly the hardware constitutions of the sensing location determination apparatus which concerns on one Embodiment. センシング箇所決定装置の機能ブロック図である。It is a functional block diagram of a sensing location determination device. 図3(a)は、下水道人孔DBのデータ構造を示す図であり、図3(b)は、下水道管路DBのデータ構造を示す図である。Fig.3 (a) is a figure which shows the data structure of sewer manhole DB, FIG.3 (b) is a figure which shows the data structure of sewer line DB. 地形DBのデータ構造を示す図である。It is a figure which shows the data structure of landform DB. 降雨量DBのデータ構造を示す図である。It is a figure which shows the data structure of rainfall DB. センシング箇所決定装置の処理を示すフローチャートである。It is a flowchart which shows the process of a sensing location determination apparatus. 下水道網の概略図である。It is the schematic of a sewer network. 図8(a)は、人孔3の水位の時系列データの一例を示すグラフであり、図8(b)は、人孔9の水位の時系列データの一例を示すグラフであり、図8(c)は、人孔17の水位の時系列データの一例を示すグラフである。FIG. 8A is a graph showing an example of the time series data of the water level of the human hole 3, and FIG. 8B is a graph showing an example of the time series data of the water level of the human hole 9. (C) is a graph showing an example of time-series data of the water level of the human hole 17. 人孔3と人孔17の相互相関関数(相関係数)と、人孔9と人孔17の相互相関関数(相関係数)とを示すグラフである。It is a graph which shows the cross correlation function (correlation coefficient) of the human hole 3 and the human hole 17, and the cross correlation function (correlation coefficient) of the human hole 9 and the human hole 17. 図6のステップS18の処理を説明するための図である。It is a figure for demonstrating the process of step S18 of FIG. 検定結果表を示す図(その1)である。It is a figure (the 1) which shows a test result table. 図12(a)、図12(b)は、検定結果表を示す図(その2)である。FIG. 12A and FIG. 12B are diagrams (part 2) showing a test result table. センサの最適設置箇所の出力例を示す図である。It is a figure which shows the example of an output of the optimal installation location of a sensor. 本実施形態の効果を説明するための図である。It is a figure for demonstrating the effect of this embodiment. 図15(a)〜図15(c)は、変形例について説明するための図である。Fig.15 (a)-FIG.15 (c) are the figures for demonstrating a modification.

以下、センシング箇所決定装置の一実施形態について、図1〜図14に基づいて詳細に説明する。本実施形態のセンシング箇所決定装置は、下水道網に含まれる複数の人孔(マンホール)の中から、水位、流量、流速等の下水状況データを検出するセンサを設置するのに適切な人孔(適切なセンシング箇所)を決定するための装置である。   Hereinafter, an embodiment of a sensing location determination device will be described in detail with reference to FIGS. The sensing location determination apparatus of the present embodiment has a human hole (appropriate for installing a sensor for detecting sewage status data such as water level, flow rate, flow velocity, etc., from a plurality of human holes (manholes) included in the sewer network. This is a device for determining an appropriate sensing location.

図1には、本実施形態に係るセンシング箇所決定装置10のハードウェア構成が示されている。センシング箇所決定装置10は、PC(Personal Computer)等の端末であり、図1に示すように、CPU(Central Processing Unit)90、ROM(Read Only Memory)92、RAM(Random Access Memory)94、記憶部(ここではHDD(Hard Disk Drive))96、ネットワークインタフェース97、表示部93、入力部95及び可搬型記憶媒体用ドライブ99等を備えている。これらセンシング箇所決定装置10の各部は、バス98に接続されている。表示部93は、液晶ディスプレイ等を含み、入力部95は、キーボードやマウス等を含む。センシング箇所決定装置10では、CPU90が、ROM92あるいはHDD96に格納されているプログラム(センシング箇所決定プログラムを含む)、或いは可搬型記憶媒体用ドライブ99が可搬型記憶媒体91から読み取ったプログラム(センシング箇所決定プログラムを含む)を実行することにより、図2に示す各部として機能する。   FIG. 1 shows a hardware configuration of a sensing location determination device 10 according to the present embodiment. The sensing location determination device 10 is a terminal such as a PC (Personal Computer), and as shown in FIG. 1, a CPU (Central Processing Unit) 90, a ROM (Read Only Memory) 92, a RAM (Random Access Memory) 94, a storage A hard disk drive (HDD (Hard Disk Drive)) 96, a network interface 97, a display unit 93, an input unit 95, a portable storage medium drive 99, and the like. Each part of the sensing location determination device 10 is connected to the bus 98. The display unit 93 includes a liquid crystal display or the like, and the input unit 95 includes a keyboard, a mouse, or the like. In the sensing location determination device 10, the CPU 90 stores a program (including a sensing location determination program) stored in the ROM 92 or the HDD 96, or a program read from the portable storage medium 91 by the portable storage medium drive 99 (sensing location determination). By executing (including a program), it functions as each unit shown in FIG.

具体的には、CPU90がプログラムを実行することで、図2に示す、入力受付部12、第1〜第3算出部としての演算部14、判定部16、及び出力部18として機能する。なお、本実施形態では、判定部16と出力部18とを含んで決定部としての機能が実現される。なお、図2には、HDD96に格納されているデータベース群(下水道人孔DB22、下水道管路DB24、地形DB26、降雨量DB28を含む)も図示されている。   Specifically, when the CPU 90 executes the program, the CPU 90 functions as the input reception unit 12, the calculation unit 14 as the first to third calculation units, the determination unit 16, and the output unit 18 illustrated in FIG. 2. In the present embodiment, the function as the determination unit is realized including the determination unit 16 and the output unit 18. FIG. 2 also shows a database group (including a sewer manhole DB 22, a sewer line DB 24, a terrain DB 26, and a rainfall DB 28) stored in the HDD 96.

入力受付部12は、ユーザが入力部95を介して入力した情報を取得し、取得した情報を演算部14や判定部16に送信する。なお、入力受付部12は、ユーザが下水道網に設置するセンサの最大個数(最大センサ個数)Nを入力した場合に、当該情報(N)を判定部16に送信する。また、入力受付部12は、最大センサ個数Nの入力を受け付けた段階で、センシング箇所の決定処理を開始する旨の指示を演算部14に送信する。   The input reception unit 12 acquires information input by the user via the input unit 95 and transmits the acquired information to the calculation unit 14 and the determination unit 16. The input receiving unit 12 transmits the information (N) to the determination unit 16 when the user inputs the maximum number (maximum number of sensors) N of sensors installed in the sewer network. In addition, the input receiving unit 12 transmits an instruction to start the sensing location determination process to the calculation unit 14 when the input of the maximum number of sensors N is received.

演算部14は、各種データベースを利用して、下水道網の各箇所(センサ設置箇所候補:本実施形態では各人孔)の有意確率を計算し、有意確率を纏めた検定結果表(図11等参照)を作成する。なお、演算部14による処理の詳細については、後述する。   The calculation unit 14 calculates the significance probability of each location (sensor installation location candidate: each human hole in the present embodiment) of the sewer network using various databases, and the test result table summarizing the significance probabilities (FIG. 11 and the like) Browse). Details of the processing by the calculation unit 14 will be described later.

判定部16は、演算部14が作成した検定結果表に基づいて、適切な位置及び数のセンサ設置箇所(センシング箇所)が定まったか否かを判定する。判定の結果、適切な位置及び数のセンシング箇所が定まった場合には、判定部16は、出力部18にその旨を通知する。一方、適切な位置及び数のセンシング箇所が定まっていない場合には、判定部16は、演算部14にその旨を通知する。なお、演算部14は、判定部16から通知を受けると、センサ設置箇所候補を減らした上で、再度、センサ設置箇所候補それぞれの有意確率を計算し、検定結果表を作成する。   The determination unit 16 determines whether or not an appropriate position and number of sensor installation locations (sensing locations) have been determined based on the test result table created by the calculation unit 14. As a result of the determination, when an appropriate position and number of sensing locations are determined, the determination unit 16 notifies the output unit 18 to that effect. On the other hand, when an appropriate position and the number of sensing locations are not determined, the determination unit 16 notifies the calculation unit 14 to that effect. In addition, when the calculation part 14 receives notification from the determination part 16, after calculating a sensor installation location candidate, it calculates again the significance probability of each sensor installation location candidate, and produces a test result table.

出力部18は、表示部93等を介して、適切な位置及び数のセンサ設置箇所(センシング箇所)の情報を出力する。   The output unit 18 outputs information on an appropriate position and number of sensor installation locations (sensing locations) via the display unit 93 and the like.

ここで、図2に示す各種データベース(22,24,26,28)のデータ構造について、図3〜図5に基づいて、詳細に説明する。   Here, the data structure of the various databases (22, 24, 26, 28) shown in FIG. 2 will be described in detail with reference to FIGS.

図3(a)には、下水道人孔DB22のデータ構造が示されている。下水道人孔DB22には、下水道網に含まれる人孔(マンホール)に関する情報が格納されている。具体的には、下水道人孔DB22は、図3(a)に示すように、「人孔番号」、「緯度」、「経度」、「人孔地盤高(標高)」、「集水面積」、「不浸透面積率」、「等価粗度」、「斜面勾配」の各フィールドを有する。「人孔番号」のフィールドには、人孔に割り当てられた識別番号が格納され、「緯度」、「経度」及び「人孔地盤高(標高)」のフィールドには、人孔が設けられた位置及び標高の情報が格納される。また、「集水面積」のフィールドには、各人孔が受け持つ土地の面積が格納され、「不浸透面積率」のフィールドには、集水面積のうち不浸透な範囲(例えば建物の屋根、舗装道路など)が占める割合が格納される。また、「等価粗度」及び「斜面勾配」のフィールドには、各人孔が受け持つ土地の等価粗度と斜面勾配の値が格納される。   FIG. 3A shows the data structure of the sewer manhole DB 22. In the sewer personal hole DB 22, information related to human holes (manholes) included in the sewer network is stored. Specifically, as shown in FIG. 3 (a), the sewer manhole DB 22 is “person number”, “latitude”, “longitude”, “manhole ground height (elevation)”, “water collection area”. , “Impermeable area ratio”, “equivalent roughness”, and “slope gradient” fields. An identification number assigned to the human hole is stored in the “human hole number” field, and a human hole is provided in the “latitude”, “longitude” and “human hole ground height (elevation)” fields. Information on position and altitude is stored. The field of “Catchment area” stores the area of land that each manhole serves, and the field of “Impermeable area ratio” stores the impervious area of the catchment area (for example, the roof of a building, Stores the percentage occupied by paved roads). In the “equivalent roughness” and “slope gradient” fields, the values of the equivalent roughness and slope gradient of the land each manhole is responsible for are stored.

図3(b)には、下水道管路DB24のデータ構造が示されている。下水道管路DB24には、人孔と人孔との間に設けられている下水道管路に関する情報が格納されている。具体的には、下水道管路DB24は、図3(b)に示すように、「人孔番号(上流、下流)」、「管路形式」、「管路長」、「管路幅」、「管路高」、「粗度係数」、「管底高(標高)(上流、下流)」の各フィールドを有する。「人孔番号(上流、下流)」のフィールドには、下水道管路の上流側の人孔及び下流側の人孔の識別番号(人孔番号)が格納される。「管路形式」のフィールドには、管路の形式(矩形等)が格納され、「管路長」のフィールドには、管路の長さが格納され、「管路幅」のフィールドには、管路の幅が格納され、「管路高」のフィールドには、管路の高さが格納される。「粗度係数」のフィールドには、管路の粗度係数が格納され、「管底高(標高)(上流、下流)」のフィールドには、管路の上流側の底部が位置する標高と、管路の下流側の底部が位置する標高とが格納される。   FIG. 3B shows the data structure of the sewer line DB 24. In the sewer line DB 24, information on the sewer line provided between the manholes is stored. Specifically, as shown in FIG. 3 (b), the sewer pipe line DB 24 includes “manhole number (upstream, downstream)”, “pipe line type”, “pipe length”, “pipe width”, Each field includes “pipe height”, “roughness coefficient”, and “pipe bottom height (elevation) (upstream, downstream)”. In the “manhole number (upstream, downstream)” field, an identification number (manhole number) of the upstream and downstream manholes of the sewer pipe is stored. In the “Pipeline type” field, the pipeline type (rectangular etc.) is stored. In the “Pipeline length” field, the length of the pipeline is stored. In the “Pipeline width” field, The width of the pipe line is stored, and the height of the pipe line is stored in the “pipe line height” field. The “roughness coefficient” field stores the roughness coefficient of the pipeline, and the “pipe bottom height (elevation) (upstream, downstream)” field indicates the altitude at which the bottom of the upstream side of the pipeline is located. The altitude at which the bottom portion on the downstream side of the pipeline is located is stored.

図4には、地形DB26のデータ構造が示されている。地形DB26は、下水道網が設けられている地区をメッシュ状に分割した複数の領域それぞれの情報を管理するデータベースである。具体的には、地形DB26は、図4に示すように、「メッシュ番号」、「地盤高(標高)」、「粗度係数」、「建物占有率」の各フィールドを有する。「メッシュ番号」のフィールドには、各領域に割り当てられた識別番号が格納される。「地盤高(標高)」のフィールドには、各領域の標高が格納される。「粗度係数」及び「建物占有率」のフィールドには、各領域の粗度係数及び建物占有率が格納される。   FIG. 4 shows the data structure of the landform DB 26. The landform DB 26 is a database that manages information on each of a plurality of areas obtained by dividing a district where a sewer network is provided into a mesh shape. Specifically, as shown in FIG. 4, the landform DB 26 has fields of “mesh number”, “ground height (elevation)”, “roughness coefficient”, and “building occupancy rate”. The identification number assigned to each area is stored in the “mesh number” field. The altitude of each area is stored in the “ground height (elevation)” field. In the fields of “roughness coefficient” and “building occupancy”, the roughness coefficient and building occupancy of each area are stored.

図5には、降雨量DB28のデータ構造が示されている。降雨量DB28は、過去における降雨量の具体的な変化を格納するデータベースである。例えば、降雨量DB28は、図5に示すように、ある日時(2014/12/24の19:00〜22:00)におけるエリアA〜Cの降雨量の変化が格納されている。なお、エリアA〜Cは、図4のメッシュ状に分割した領域を複数含むエリアである。   FIG. 5 shows the data structure of the rainfall DB 28. The rainfall DB 28 is a database that stores specific changes in rainfall in the past. For example, as shown in FIG. 5, the rainfall DB 28 stores changes in rainfall in areas A to C at a certain date (2014/12/24 19:00:00 to 22:00). Areas A to C are areas including a plurality of areas divided into meshes in FIG.

なお、図3〜図5の各種データベースのデータ構造は、一例である。すなわち、各種データベースが有するフィールドの内容については、種々変更可能である。また、必要に応じて、他のデータベースを追加等してもよい。   The data structures of the various databases in FIGS. 3 to 5 are examples. That is, various changes can be made to the contents of the fields of various databases. Moreover, you may add another database as needed.

(センシング箇所決定装置10の処理)
次に、図6のフローチャートに沿って、その他図面を適宜参照しつつ、センシング箇所決定装置10の処理について、詳細に説明する。なお、本実施形態では、一例として、図7に示すような下水道網においてセンシング箇所(センサの最適設置箇所)を決定する場合について説明する。なお、図7において、下水道管が分岐している場合は、太い管(すなわち幹線)のみにセンサを置くと仮定する。このような条件で、仮に全ての人孔にセンサを設置すると、図7に番号を付して示すような最大17個のセンサ設置箇所候補がある。本実施形態では、センシング箇所決定装置10は、これらのセンサ設置箇所候補の中から、ユーザが定めた最大センサ個数以下の適切な箇所をセンシング箇所として決定するものとする。
(Processing of sensing location determination apparatus 10)
Next, the processing of the sensing location determination device 10 will be described in detail along the flowchart of FIG. 6 with reference to other drawings as appropriate. In addition, in this embodiment, the case where a sensing location (optimum installation location of a sensor) is determined in a sewer network as shown in FIG. 7 is demonstrated as an example. In FIG. 7, when the sewer pipe is branched, it is assumed that the sensor is placed only on the thick pipe (that is, the main line). Under such conditions, if sensors are installed in all human holes, there are a maximum of 17 sensor installation location candidates as shown in FIG. In this embodiment, the sensing location determination apparatus 10 shall determine an appropriate location below the maximum number of sensors determined by the user as a sensing location from these sensor installation location candidates.

図6の処理では、まず、ステップS10において、入力受付部12が、最大センサ個数Nの入力を受け付ける。この場合、ユーザが入力部95を介して最大センサ個数Nを入力すると、入力受付部12は、当該最大センサ個数Nを受け付け、判定部16に送信する。また、入力受付部12は、最大センサ個数Nの入力を受け付けた段階で、演算部14に対して、処理を開始する旨の指示を出す。なお、本実施形態では、ユーザは最大センサ個数Nとして「5」を入力したものとする。   In the process of FIG. 6, first, in step S <b> 10, the input receiving unit 12 receives an input of the maximum number of sensors N. In this case, when the user inputs the maximum sensor number N through the input unit 95, the input reception unit 12 receives the maximum sensor number N and transmits it to the determination unit 16. Further, the input receiving unit 12 gives an instruction to start the processing to the calculation unit 14 when the input of the maximum number of sensors N is received. In the present embodiment, it is assumed that the user inputs “5” as the maximum sensor number N.

次いで、ステップS12では、演算部14が、下水道管渠、地形、降雨量の情報に基づいて下水データを算出する。具体的には、演算部14は、下水道人孔DB22、下水道管路DB24、地形DB26、降雨量DB28の各データに基づいて、公知の演算方法により、下水道網の各箇所(各人孔)の水位、流量、流速等の時系列データの算出を行う。なお、本実施形態では、一例として、演算部14は、水位[m]の時系列データを算出したものとする。図8(a)〜図8(c)には、一例として、図7の人孔3、人孔9、人孔17における水位の時系列データが示されている。   Next, in step S12, the calculation unit 14 calculates sewage data based on information on sewer pipes, topography, and rainfall. Specifically, the calculation unit 14 calculates each location (each human hole) of the sewer network by a known calculation method based on each data of the sewer manhole DB22, the sewer pipe DB24, the topography DB26, and the rainfall DB28. Calculate time-series data such as water level, flow rate, and flow velocity. In the present embodiment, as an example, the calculation unit 14 calculates time-series data of the water level [m]. 8A to 8C show, as an example, time series data of water levels in the human hole 3, the human hole 9, and the human hole 17 in FIG.

次いで、ステップS14では、演算部14が、算出した時系列データを用いて各人孔と最下流の人孔(ここでは、人孔17とする)の相互相関関数R(τ)を計算する。なお、τは、各箇所(人孔)の時系列データと最下流の箇所(人孔17)の時系列データの変化の時間差を意味する。図9には、人孔3と人孔17の相互相関関数(相関係数)と、人孔9と人孔17の相互相関関数(相関係数)とが示されている。より詳細には、図9には、人孔3や人孔9の時系列データをタイムシフトさせた場合の、人孔17の時系列データとの相関の変化が示されている。   Next, in step S14, the computing unit 14 calculates a cross-correlation function R (τ) between each human hole and the most downstream human hole (here, human hole 17) using the calculated time series data. Note that τ means the time difference between changes in the time series data at each location (human hole) and the time series data at the most downstream location (human hole 17). FIG. 9 shows a cross-correlation function (correlation coefficient) between the human hole 3 and the human hole 17 and a cross-correlation function (correlation coefficient) between the human hole 9 and the human hole 17. More specifically, FIG. 9 shows a change in correlation with the time series data of the human hole 17 when the time series data of the human hole 3 and the human hole 9 are time-shifted.

次いで、ステップS16では、演算部14が、相互相関関数の極大値を計算する。すなわち、演算部14は、相関係数が正になり、かつ、時間差の相互相関関数の二階微分(d2R(τ)/dτ2)が0になる最も小さい時間差τを求める。なお、求められた時間差τは、最下流の人孔17の、他の人孔に対する遅れ時間を意味する。図9では、人孔3に対する人孔17の遅れ時間がT3、人孔9に対する人孔17の遅れ時間がT9として示されている。 Next, in step S16, the calculation unit 14 calculates the maximum value of the cross-correlation function. That is, the calculation unit 14 obtains the smallest time difference τ in which the correlation coefficient is positive and the second-order derivative (d 2 R (τ) / dτ 2 ) of the time difference cross-correlation function is zero. The obtained time difference τ means the delay time of the most downstream human hole 17 with respect to other human holes. In FIG. 9, the delay time of the human hole 17 with respect to the human hole 3 is shown as T3, and the delay time of the human hole 17 with respect to the human hole 9 is shown as T9.

次いで、ステップS18では、演算部14が、各時系列データを極大値(遅れ時間τ)だけシフトさせる。この場合、演算部14は、図10の上段の図に示すように、人孔3の時系列データを遅れ時間T3だけ右側にタイムシフトさせる(破線のグラフ参照)。また、演算部14は、図10の中段の図に示すように、人孔9の時系列データを遅れ時間T9だけ右側にタイムシフトさせる(破線のグラフ参照)。なお、その他の人孔の時系列データについても同様にタイムシフトさせる。   Next, in step S18, the calculation unit 14 shifts each time series data by the maximum value (delay time τ). In this case, as shown in the upper diagram of FIG. 10, the calculation unit 14 time-shifts the time series data of the human hole 3 to the right by the delay time T3 (see the broken line graph). Further, as shown in the middle diagram of FIG. 10, the calculation unit 14 time-shifts the time series data of the human hole 9 to the right by the delay time T9 (see the broken line graph). It should be noted that time shift is similarly performed for other time series data of human holes.

次いで、ステップS20では、演算部14が、有意確率を計算する。この場合、演算部14は、遅れ時間分だけタイムシフトした各人孔の時系列データに対し、統計的仮説検定を行い、有意確率を計算する。   Next, in step S20, the calculation unit 14 calculates a significance probability. In this case, the calculation unit 14 performs a statistical hypothesis test on the time series data of each human hole shifted by the delay time, and calculates a significance probability.

次いで、ステップS22では、演算部14が、センサ設置箇所候補(人孔)の検定結果表を作成する。図11には、検定結果表の一例が示されている。図11の検定結果表においては、センサ設置箇所候補(ここでは、人孔1〜16)のt値、P(Probability)−値が纏められている。   Next, in step S22, the computing unit 14 creates a test result table for sensor installation location candidates (human holes). FIG. 11 shows an example of the test result table. In the test result table of FIG. 11, t values and P (Probability) -values of sensor installation location candidates (here, human holes 1 to 16) are summarized.

次いで、ステップS24では、判定部16が、全ての有意確率が有意であるか否かを判断する。この場合、例えば、判定部16は、全てのP−値が、予め定められている閾値(例えば、0.05とする)以下であるか否かに基づいて、全ての有意確率が有意であるか否かを判断する。ここでの判断が否定された場合、すなわち、P−値が閾値を超えるセンサ設置箇所候補が存在する場合には、ステップS28に移行する。ステップS28では、演算部14が、センサ設置箇所候補を減じる処理を行う。この場合、演算部14は、ステップS28において閾値を超える最も大きなP−値を有する人孔(図11では、太破線枠で囲んだ人孔2)をセンサ設置箇所候補から除外した後、ステップS20に戻る。   Next, in step S24, the determination unit 16 determines whether or not all significant probabilities are significant. In this case, for example, the determination unit 16 determines that all significant probabilities are significant based on whether or not all P-values are equal to or less than a predetermined threshold value (for example, 0.05). Determine whether or not. If the determination is negative, that is, if there is a sensor installation location candidate whose P-value exceeds the threshold, the process proceeds to step S28. In step S28, the calculating part 14 performs the process which subtracts a sensor installation location candidate. In this case, the calculation unit 14 excludes the human hole having the largest P-value exceeding the threshold value in step S28 (in FIG. 11, the human hole 2 surrounded by the thick broken line frame) from the sensor installation location candidates, and then performs step S20. Return to.

ステップS20に戻ると、演算部14は、残りのセンサ設置箇所候補(ここでは、人孔1、人孔3〜16)の有意確率を再度計算する。そして、次のステップS22では、演算部14が、残りのセンサ設置箇所候補(人孔)の検定結果表を再度作成する。この場合、図12(a)に示すように、人孔1、人孔3〜16に関する検定結果表が作成される。次いで、ステップS24では、判定部16が、全ての有意確率が有意であるか否かを再度判断する。このステップS24の判断が否定されると、ステップS28に移行する。なお、ステップS28において図12(a)の場合、演算部14は、太破線枠で囲んだ人孔10を除外し、ステップS20に戻る。   Returning to step S20, the calculation unit 14 calculates again the significance probability of the remaining sensor installation location candidates (here, human hole 1, human holes 3 to 16). In the next step S22, the calculation unit 14 again creates a test result table for the remaining sensor installation location candidates (human holes). In this case, as shown in FIG. 12A, a test result table for the human hole 1 and the human holes 3 to 16 is created. Next, in step S24, the determination unit 16 determines again whether or not all significant probabilities are significant. If the determination in step S24 is negative, the process proceeds to step S28. In the case of FIG. 12A in step S28, the calculation unit 14 excludes the human hole 10 surrounded by the thick broken line frame and returns to step S20.

なお、ステップS24の判断が肯定された場合、すなわち、全ての有意確率が有意であった場合には、判定部16は、ステップS26に移行する。ステップS26に移行すると、判定部16は、センサ設置箇所候補がN以下であるか否かを判断する。ここでの判断が否定された場合には、ステップS28に移行し、検定結果表に含まれるセンサ設置箇所候補のうち、最も有意確率(P−値)が大きい人孔を除外し、ステップS20に戻る。   If the determination in step S24 is affirmative, that is, if all significant probabilities are significant, the determination unit 16 proceeds to step S26. If transfering it to step S26, the determination part 16 will judge whether a sensor installation location candidate is N or less. If the determination is negative, the process proceeds to step S28, and the human hole having the highest significance probability (P-value) is excluded from the sensor installation location candidates included in the test result table, and the process proceeds to step S20. Return.

以降、ステップS26の判断が肯定されるまで、ステップS20〜S28の処理・判断が繰り返し実行される。なお、ステップS26の判断が肯定される段階における検定結果表の例が、図12(b)に示されている。図12(b)においては、全ての有意確率(P−値)が閾値(0.05)未満となり、センサ設置箇所候補(人孔)の数がN(=5)以下になっている。   Thereafter, the processes and determinations in steps S20 to S28 are repeatedly executed until the determination in step S26 is affirmed. An example of the test result table at the stage where the determination in step S26 is affirmed is shown in FIG. In FIG. 12B, all significant probabilities (P-values) are less than the threshold (0.05), and the number of sensor installation location candidates (human holes) is N (= 5) or less.

ステップS26の判断が肯定され、ステップS30に移行すると、出力部18は、センサの最適設置箇所を決定し、出力する。この場合、出力部18は、図12(b)の検定計画表に含まれている人孔(人孔6,9,11,13,16)をセンサの最適設置箇所として、表示部93等を介して出力する。この場合、出力部18は、図13に示すような下水道網を示す図を用いて、センサの最適設置箇所(図13の黒丸参照)を出力するようにしてもよい。   When the determination in step S26 is affirmed and the process proceeds to step S30, the output unit 18 determines and outputs the optimum installation location of the sensor. In this case, the output unit 18 uses the human holes (human holes 6, 9, 11, 13, 16) included in the test plan table of FIG. Output via. In this case, the output unit 18 may output an optimum installation location of the sensor (see a black circle in FIG. 13) using a diagram showing a sewer network as shown in FIG.

以上、詳細に説明したように、本実施形態によると、演算部14は、各種データベース(22,24,26,28)に基づき、下水道網の複数の箇所(人孔)における水位の時系列データを算出し(S12)、下水道網の特定箇所(最下流の人孔17)における時系列データの、その他の人孔における時系列データそれぞれに対する遅れ時間τを算出し(S14,S16)、算出した遅れ時間τと、その他の人孔それぞれにおける時系列データとを用いて、統計的仮説検定によりその他の人孔それぞれの有意確率を算出する(S18、S20)。そして、判定部16及び出力部18は、算出された有意確率に基づいてセンサの最適設置箇所(センシング箇所)を決定する(S22〜S30)。このように、本実施形態では、遅れ時間を考慮して、センシング箇所を決定することで、適切なセンシング箇所を決定することができる。図14には、本実施形態と比較例(遅れ時間を考慮せずに統計的仮説検定を行った例)における、設置センサ数による標準誤差の推移が示されている。なお、標準誤差は、設置したセンサを用いて最下流の人孔17の水位を予測した場合の誤差を意味する。本実施形態と比較例とを比較すると、本実施形態のほうが比較例よりも、全体的にモデル精度が向上している(標準誤差が小さくなっている)ことが確認できる。また、比較例では、センサ数を削減した場合に標準誤差が急激に大きくなるが、本実施形態では緩やかに標準誤差が大きくなることがわかる。これらより、本実施形態では、モデルが安定しており、たとえセンサの数を減らしても、モデル精度への影響が小さいことがわかる。このため、本実施形態では、図12(a)や図13に示すようにセンシング箇所をユーザが入力した最大センサ個数N(=5)以下に減らしても、精度良く最下流の人孔の水位を推定することが可能となる。したがって、本実施形態では、センサの数を極力減らすことにより、下水道網の状況を予測する際に処理する必要のあるデータ量を削減することができるとともに、センサ設置費用の削減を図ることが可能となる。   As described above in detail, according to this embodiment, the calculation unit 14 is based on various databases (22, 24, 26, 28), and time-series data of water levels at a plurality of locations (human holes) in the sewer network. (S12), and calculating the delay time τ of the time series data in the specific part of the sewer network (the most downstream human hole 17) with respect to each of the time series data in the other human holes (S14, S16). Using the delay time τ and the time series data in each of the other foramen, the significance probability for each of the other foramen is calculated by a statistical hypothesis test (S18, S20). And the determination part 16 and the output part 18 determine the optimal installation location (sensing location) of a sensor based on the calculated significant probability (S22-S30). Thus, in this embodiment, an appropriate sensing location can be determined by determining the sensing location in consideration of the delay time. FIG. 14 shows the transition of the standard error depending on the number of installed sensors in the present embodiment and a comparative example (an example in which a statistical hypothesis test is performed without considering delay time). The standard error means an error when the water level of the most downstream human hole 17 is predicted using the installed sensor. Comparing the present embodiment with the comparative example, it can be confirmed that the model accuracy is generally improved (standard error is smaller) in the present embodiment than in the comparative example. In the comparative example, the standard error increases rapidly when the number of sensors is reduced, but it can be seen that the standard error gradually increases in this embodiment. From these, it can be seen that in this embodiment, the model is stable, and even if the number of sensors is reduced, the influence on the model accuracy is small. For this reason, in this embodiment, as shown in FIG. 12A and FIG. 13, even if the sensing location is reduced to the number N (= 5) or less of the maximum number of sensors input by the user, the water level of the most downstream human hole is accurately obtained. Can be estimated. Therefore, in this embodiment, by reducing the number of sensors as much as possible, it is possible to reduce the amount of data that needs to be processed when predicting the status of the sewer network, and to reduce the cost of installing the sensors It becomes.

また、本実施形態では、遅れ時間τを算出する際、最下流以外の人孔17における時系列データをタイムシフトさせた場合に、最下流の人孔17における時系列データとの相関が最も大きくなるシフト時間を遅れ時間として算出する。これにより、遅れ時間τとして適切な値を算出することができる。   Further, in the present embodiment, when calculating the delay time τ, when the time series data in the human hole 17 other than the most downstream is time-shifted, the correlation with the time series data in the most downstream human hole 17 is the largest. Is calculated as a delay time. Thereby, an appropriate value can be calculated as the delay time τ.

なお、上記実施形態では、演算部14は、ステップS12において、下水道網のセンサ設置箇所候補における水位の時系列データを算出する場合について説明したが、これに限られるものではない。例えば、演算部14は、流量の時系列データを算出することとしてもよい。この場合、図15(a)〜図15(c)に示すように、水位の時系列データ(図8(a)〜図8(c)参照)と同様のデータが得られるので、演算部14及び判定部16は、上記実施形態と同様の処理により、センサの最適設置箇所を決定するようにすればよい。また、演算部14は、ステップS12において、流速の時系列データを算出することとしてもよい。この場合、演算部14は、人孔間それぞれの距離を人孔間における下水の流速で除した値(下水が人孔間を流れる時間)を求め、当該値を最下流の人孔まで積算することで、遅れ時間を算出するようにすればよい。なお、上記実施形態では、水位、流量、流速のうち複数の時系列データを算出し、各時系列データを用いて、図6の処理をそれぞれ実行し、各時系列データから求まるセンサの最適設置箇所から、最終的なセンサの最適設置箇所を決定するようにしてもよい。   In addition, although the calculating part 14 demonstrated the case where the calculation part 14 calculated the time series data of the water level in the sensor installation location candidate of a sewer network in step S12, it is not restricted to this. For example, the calculation unit 14 may calculate time-series data of the flow rate. In this case, as shown in FIGS. 15A to 15C, data similar to the water level time-series data (see FIGS. 8A to 8C) is obtained. And the determination part 16 should just determine the optimal installation location of a sensor by the process similar to the said embodiment. Moreover, the calculating part 14 is good also as calculating the time series data of the flow velocity in step S12. In this case, the calculation unit 14 obtains a value obtained by dividing the distance between the human holes by the flow rate of the sewage between the human holes (the time during which the sewage flows between the human holes), and integrates the value to the most downstream human hole. Thus, the delay time may be calculated. In the above embodiment, a plurality of time series data is calculated from the water level, the flow rate, and the flow velocity, and the processing shown in FIG. 6 is executed using each time series data, and the optimum sensor installation determined from each time series data is performed. The final optimal installation location of the sensor may be determined from the location.

なお、上記実施形態では、ステップS24において、判定部16は、P−値を用いて各人孔の有意性を判定する場合について説明したが、これに限らず、t値を用いることとしてもよい。また、P−値やt値に限らず、仮説検定で用いられるF値や標準誤差などを用いてもよい。また、図6では、ユーザが最大センサ個数Nを入力し、センシング箇所決定装置10は、センサの最適設置箇所がN箇所以下となるように決定する場合について説明したが、これに限られるものではない。例えば、検定結果表において、全変数が有意となった場合(ステップS24が肯定された場合)に、検定結果表に残っている人孔をセンサの最適設置箇所として決定するようにしてもよい。   In the above embodiment, the determination unit 16 has described the case where the significance of each human hole is determined using the P-value in step S24. However, the present invention is not limited to this, and a t value may be used. . Further, not only the P-value and the t value, but also an F value and standard error used in the hypothesis test may be used. In FIG. 6, the user inputs the maximum number of sensors N, and the sensing location determination device 10 has been described to determine that the optimal installation location of the sensor is N or less. However, the present invention is not limited to this. Absent. For example, when all the variables are significant in the test result table (when step S24 is affirmed), the human hole remaining in the test result table may be determined as the optimum installation location of the sensor.

なお、上記実施形態では、センシング箇所決定装置10がインターネットなどのネットワークに接続されたサーバであってもよい。この場合、センシング箇所決定装置10は、ネットワークに接続された端末(クライアント)からの入力に応じて、図6の処理を実行し、処理結果を端末(クライアント)に送信するようにしてもよい。   In the above embodiment, the sensing location determination device 10 may be a server connected to a network such as the Internet. In this case, the sensing location determination device 10 may execute the processing of FIG. 6 in accordance with an input from a terminal (client) connected to the network, and transmit the processing result to the terminal (client).

なお、上記の処理機能は、コンピュータによって実現することができる。その場合、処理装置が有すべき機能の処理内容を記述したプログラムが提供される。そのプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体(ただし、搬送波は除く)に記録しておくことができる。   The above processing functions can be realized by a computer. In that case, a program describing the processing contents of the functions that the processing apparatus should have is provided. By executing the program on a computer, the above processing functions are realized on the computer. The program describing the processing contents can be recorded on a computer-readable recording medium (except for a carrier wave).

プログラムを流通させる場合には、例えば、そのプログラムが記録されたDVD(Digital Versatile Disc)、CD−ROM(Compact Disc Read Only Memory)などの可搬型記録媒体の形態で販売される。また、プログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することもできる。   When the program is distributed, for example, it is sold in the form of a portable recording medium such as a DVD (Digital Versatile Disc) or a CD-ROM (Compact Disc Read Only Memory) on which the program is recorded. It is also possible to store the program in a storage device of a server computer and transfer the program from the server computer to another computer via a network.

プログラムを実行するコンピュータは、例えば、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、自己の記憶装置に格納する。そして、コンピュータは、自己の記憶装置からプログラムを読み取り、プログラムに従った処理を実行する。なお、コンピュータは、可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することもできる。また、コンピュータは、サーバコンピュータからプログラムが転送されるごとに、逐次、受け取ったプログラムに従った処理を実行することもできる。   The computer that executes the program stores, for example, the program recorded on the portable recording medium or the program transferred from the server computer in its own storage device. Then, the computer reads the program from its own storage device and executes processing according to the program. The computer can also read the program directly from the portable recording medium and execute processing according to the program. Further, each time the program is transferred from the server computer, the computer can sequentially execute processing according to the received program.

上述した実施形態は本発明の好適な実施の例である。但し、これに限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変形実施可能である。   The above-described embodiment is an example of a preferred embodiment of the present invention. However, the present invention is not limited to this, and various modifications can be made without departing from the scope of the present invention.

なお、以上の実施形態の説明に関して、更に以下の付記を開示する。
(付記1) 下水道網に関するデータと降雨量に関するデータに基づき、前記下水道網の複数の箇所における水位、流量、流速の少なくとも1つの時系列データを算出し、
前記下水道網の特定箇所における前記時系列データの、該特定箇所以外の箇所における前記時系列データそれぞれに対する遅れ時間を算出し、
算出した前記遅れ時間と、前記特定箇所以外の箇所それぞれにおける前記時系列データとを用いて、統計的仮説検定により前記特定箇所以外の箇所それぞれの有意確率を算出し、
算出された前記有意確率に基づいてセンサを配置するセンシング箇所を決定する、
処理をコンピュータが実行することを特徴とする下水道網のセンシング箇所決定方法。
(付記2) 前記決定する処理では、算出された前記有意確率が満たすべき条件と、予め設定されたセンシング箇所の最大数と、に基づいて、前記センシング箇所を決定する、ことを特徴とする付記1に記載の下水道網のセンシング箇所決定方法。
(付記3) 前記遅れ時間を算出する処理では、前記特定箇所以外の箇所における前記時系列データをタイムシフトさせた場合に、前記特定箇所における前記時系列データとの相関が最も大きくなるシフト時間を遅れ時間として算出する、ことを特徴とする付記1又は2に記載の下水道網のセンシング箇所決定方法。
(付記4) 前記時系列データを算出する処理において、前記下水道網の複数の箇所における流速の時系列データを算出する場合に、
前記遅れ時間を算出する処理では、前記流速の時系列データと、前記複数の箇所相互間の距離とに基づいて、遅れ時間を算出することを特徴とする付記1又は2に記載の下水道網のセンシング箇所決定方法。
(付記5) 下水道網に関するデータと降雨量に関するデータに基づき、前記下水道網の複数の箇所における水位、流量、流速の少なくとも1つの時系列データを算出し、
前記下水道網の特定箇所における前記時系列データの、該特定箇所以外の箇所における前記時系列データそれぞれに対する遅れ時間を算出し、
算出した前記遅れ時間と、前記特定箇所以外の箇所それぞれにおける前記時系列データとを用いて、統計的仮説検定により前記特定箇所以外の箇所それぞれの有意確率を算出し、
算出された前記有意確率に基づいてセンサを配置するセンシング箇所を決定する、
処理をコンピュータに実行させることを特徴とする下水道網のセンシング箇所決定プログラム。
(付記6) 前記決定する処理では、算出された前記有意確率が満たすべき条件と、予め設定されたセンシング箇所の最大数と、に基づいて、前記センシング箇所を決定する、ことを特徴とする付記5に記載の下水道網のセンシング箇所決定プログラム。
(付記7) 前記遅れ時間を算出する処理では、前記特定箇所以外の箇所における前記時系列データをタイムシフトさせた場合に、前記特定箇所における前記時系列データとの相関が最も大きくなるシフト時間を遅れ時間として算出する、ことを特徴とする付記5又は6に記載の下水道網のセンシング箇所決定プログラム。
(付記8) 前記時系列データを算出する処理において、前記下水道網の複数の箇所における流速の時系列データを算出する場合に、
前記遅れ時間を算出する処理では、前記流速の時系列データと、前記複数の箇所相互間の距離とに基づいて、遅れ時間を算出することを特徴とする付記5又は6に記載の下水道網のセンシング箇所決定プログラム。
(付記9) 下水道網に関するデータと降雨量に関するデータに基づき、前記下水道網の複数の箇所における水位、流量、流速の少なくとも1つの時系列データを算出する第1算出部と、
前記下水道網の特定箇所における前記時系列データの、該特定箇所以外の箇所における前記時系列データそれぞれに対する遅れ時間を算出する第2算出部と、
算出した前記遅れ時間と、前記特定箇所以外の箇所それぞれにおける前記時系列データとを用いて、統計的仮説検定により前記特定箇所以外の箇所それぞれの有意確率を算出する第3算出部と、
算出された前記有意確率に基づいてセンサを配置するセンシング箇所を決定する決定部と、を備える下水道網のセンシング箇所決定装置。
(付記10) 前記決定部は、算出された前記有意確率が満たすべき条件と、予め設定されたセンシング箇所の最大数と、に基づいて、前記センシング箇所を決定する、ことを特徴とする付記9に記載の下水道網のセンシング箇所決定装置。
(付記11) 前記第2算出部は、前記特定箇所以外の箇所における前記時系列データをタイムシフトさせた場合に、前記特定箇所における前記時系列データとの相関が最も大きくなるシフト時間を遅れ時間として算出する、ことを特徴とする付記9又は10に記載の下水道網のセンシング箇所決定装置。
(付記12) 前記第1算出部が、前記下水道網の複数の箇所における流速の時系列データを算出する場合に、
前記第2算出部は、前記流速の時系列データと、前記複数の箇所相互間の距離とに基づいて、遅れ時間を算出することを特徴とする付記9又は10に記載の下水道網のセンシング箇所決定装置。
In addition, the following additional remarks are disclosed regarding description of the above embodiment.
(Supplementary note 1) Based on the data on the sewer network and the data on the rainfall, at least one time series data of the water level, the flow rate and the flow velocity in the plurality of locations of the sewer network is calculated.
Calculate the delay time for each of the time-series data in a location other than the specific location of the time-series data in a specific location of the sewer network,
Using the calculated delay time and the time-series data in each location other than the specific location, calculate the significance probability of each location other than the specific location by statistical hypothesis test,
Determining a sensing location to place the sensor based on the calculated significance;
A method for determining a sensing location of a sewer network, wherein the processing is performed by a computer.
(Additional remark 2) In the said process to determine, the said sensing location is determined based on the conditions which the calculated said significant probability should satisfy | fill, and the preset maximum number of sensing locations. The sensing location determination method of the sewer network of 1.
(Additional remark 3) In the process which calculates the said delay time, when the time series data in parts other than the said specific part are time-shifted, the shift time when the correlation with the said time series data in the said specific part becomes the largest is used. The sensing location determination method of the sewer network according to appendix 1 or 2, characterized in that it is calculated as a delay time.
(Additional remark 4) In the process which calculates the said time series data, when calculating the time series data of the flow velocity in the several location of the said sewer network,
In the process of calculating the delay time, the delay time is calculated based on the time-series data of the flow velocity and the distance between the plurality of locations. Sensing location determination method.
(Appendix 5) Based on the data on the sewer network and the data on the rainfall, calculate at least one time-series data of the water level, flow rate, and flow velocity at a plurality of locations of the sewer network,
Calculate the delay time for each of the time-series data in a location other than the specific location of the time-series data in a specific location of the sewer network,
Using the calculated delay time and the time-series data in each location other than the specific location, calculate the significance probability of each location other than the specific location by statistical hypothesis test,
Determining a sensing location to place the sensor based on the calculated significance;
A sewer network sensing location determination program that causes a computer to execute processing.
(Additional remark 6) In the said process to determine, the said sensing location is determined based on the conditions which the calculated said significant probability should satisfy | fill, and the preset maximum number of sensing locations, It is characterized by the above-mentioned. 5. Sewer network sensing location determination program according to 5.
(Additional remark 7) In the process which calculates the said delay time, when the time series data in parts other than the said specific part are time-shifted, the shift time when the correlation with the said time series data in the said specific part becomes the largest is used. The sewer network sensing location determination program according to appendix 5 or 6, characterized in that it is calculated as a delay time.
(Additional remark 8) In the process which calculates the said time series data, when calculating the time series data of the flow velocity in the several location of the said sewer network,
In the process of calculating the delay time, the delay time is calculated based on the time-series data of the flow velocity and the distance between the plurality of locations. Sensing location determination program.
(Additional remark 9) The 1st calculation part which calculates at least 1 time series data of the water level in the several location of the said sewer network, flow volume, and the flow velocity based on the data regarding a sewer network, and the data regarding rainfall,
A second calculation unit that calculates a delay time for each of the time series data in a place other than the specific part of the time series data in the specific part of the sewer network;
Using the calculated delay time and the time series data in each of the locations other than the specific location, a third calculation unit for calculating the significance probability of each location other than the specific location by a statistical hypothesis test;
A sewer network sensing location determination apparatus, comprising: a determination unit that determines a sensing location where a sensor is arranged based on the calculated significance probability.
(Additional remark 10) The said determination part determines the said sensing location based on the conditions which the calculated said significant probability should satisfy | fill, and the preset maximum number of sensing locations, The additional remark 9 characterized by the above-mentioned. Sensing location determination device for sewer network described in 1.
(Additional remark 11) When the said 2nd calculation part carries out time shift of the said time series data in places other than the said specific location, the shift time when the correlation with the said time series data in the said specific location becomes the largest is a delay time. The sewer network sensing location determination device according to appendix 9 or 10, characterized in that it is calculated as:
(Supplementary Note 12) When the first calculation unit calculates time-series data of flow velocities at a plurality of locations of the sewer network,
The second calculation unit calculates a delay time based on the time-series data of the flow velocity and the distance between the plurality of locations, and the sensing location of the sewer network according to appendix 9 or 10, Decision device.

10 センシング箇所決定装置
14 演算部(第1〜第3算出部)
16 判定部(決定部の一部)
18 出力部(決定部の一部)
90 コンピュータ
10 Sensing Location Determination Device 14 Calculation Unit (First to Third Calculation Units)
16 Judgment part (part of decision part)
18 Output part (part of the decision part)
90 computers

Claims (6)

下水道網に関するデータと降雨量に関するデータに基づき、前記下水道網の複数の箇所における水位、流量、流速の少なくとも1つの時系列データを算出し、
前記下水道網の特定箇所における前記時系列データの、該特定箇所以外の箇所における前記時系列データそれぞれに対する遅れ時間を算出し、
算出した前記遅れ時間と、前記特定箇所以外の箇所それぞれにおける前記時系列データとを用いて、統計的仮説検定により前記特定箇所以外の箇所それぞれの有意確率を算出し、
算出された前記有意確率に基づいてセンサを配置するセンシング箇所を決定する、
処理をコンピュータが実行することを特徴とする下水道網のセンシング箇所決定方法。
Based on the data on the sewer network and the data on the rainfall, calculate at least one time-series data of the water level, the flow rate and the flow velocity at a plurality of locations of the sewer network,
Calculate the delay time for each of the time-series data in a location other than the specific location of the time-series data in a specific location of the sewer network,
Using the calculated delay time and the time-series data in each location other than the specific location, calculate the significance probability of each location other than the specific location by statistical hypothesis test,
Determining a sensing location to place the sensor based on the calculated significance;
A method for determining a sensing location of a sewer network, wherein the processing is performed by a computer.
前記決定する処理では、算出された前記有意確率が満たすべき条件と、予め設定されたセンシング箇所の最大数と、に基づいて、前記センシング箇所を決定する、ことを特徴とする請求項1に記載の下水道網のセンシング箇所決定方法。   The said determination process WHEREIN: The said sensing location is determined based on the conditions which the calculated said significant probability should satisfy | fill, and the preset maximum number of sensing locations, The said sensing location is characterized by the above-mentioned. To determine the sensing location of the sewer network. 前記遅れ時間を算出する処理では、前記特定箇所以外の箇所における前記時系列データをタイムシフトさせた場合に、前記特定箇所における前記時系列データとの相関が最も大きくなるシフト時間を遅れ時間として算出する、ことを特徴とする請求項1又は2に記載の下水道網のセンシング箇所決定方法。   In the process of calculating the delay time, when the time-series data at a location other than the specific location is time-shifted, a shift time that maximizes the correlation with the time-series data at the specific location is calculated as the delay time. The sensing location determination method of the sewer network of Claim 1 or 2 characterized by the above-mentioned. 前記時系列データを算出する処理において、前記下水道網の複数の箇所における流速の時系列データを算出する場合に、
前記遅れ時間を算出する処理では、前記流速の時系列データと、前記複数の箇所相互間の距離とに基づいて、遅れ時間を算出することを特徴とする請求項1又は2に記載の下水道網のセンシング箇所決定方法。
In the process of calculating the time series data, when calculating the time series data of the flow velocity at a plurality of locations of the sewer network,
The sewer network according to claim 1 or 2, wherein, in the process of calculating the delay time, the delay time is calculated based on the time-series data of the flow velocity and the distance between the plurality of locations. Sensing location determination method.
下水道網に関するデータと降雨量に関するデータに基づき、前記下水道網の複数の箇所における水位、流量、流速の少なくとも1つの時系列データを算出し、
前記下水道網の特定箇所における前記時系列データの、該特定箇所以外の箇所における前記時系列データそれぞれに対する遅れ時間を算出し、
算出した前記遅れ時間と、前記特定箇所以外の箇所それぞれにおける前記時系列データとを用いて、統計的仮説検定により前記特定箇所以外の箇所それぞれの有意確率を算出し、
算出された前記有意確率に基づいてセンサを配置するセンシング箇所を決定する、
処理をコンピュータに実行させることを特徴とする下水道網のセンシング箇所決定プログラム。
Based on the data on the sewer network and the data on the rainfall, calculate at least one time-series data of the water level, the flow rate and the flow velocity at a plurality of locations of the sewer network,
Calculate the delay time for each of the time-series data in a location other than the specific location of the time-series data in a specific location of the sewer network,
Using the calculated delay time and the time-series data in each location other than the specific location, calculate the significance probability of each location other than the specific location by statistical hypothesis test,
Determining a sensing location to place the sensor based on the calculated significance;
A sewer network sensing location determination program that causes a computer to execute processing.
下水道網に関するデータと降雨量に関するデータに基づき、前記下水道網の複数の箇所における水位、流量、流速の少なくとも1つの時系列データを算出する第1算出部と、
前記下水道網の特定箇所における前記時系列データの、該特定箇所以外の箇所における前記時系列データそれぞれに対する遅れ時間を算出する第2算出部と、
算出した前記遅れ時間と、前記特定箇所以外の箇所それぞれにおける前記時系列データとを用いて、統計的仮説検定により前記特定箇所以外の箇所それぞれの有意確率を算出する第3算出部と、
算出された前記有意確率に基づいてセンサを配置するセンシング箇所を決定する決定部と、を備える下水道網のセンシング箇所決定装置。
A first calculation unit that calculates at least one time-series data of a water level, a flow rate, and a flow velocity at a plurality of locations of the sewer network based on data on the sewer network and rainfall data;
A second calculation unit that calculates a delay time for each of the time series data in a place other than the specific part of the time series data in the specific part of the sewer network;
Using the calculated delay time and the time series data in each of the locations other than the specific location, a third calculation unit for calculating the significance probability of each location other than the specific location by a statistical hypothesis test;
A sewer network sensing location determination apparatus, comprising: a determination unit that determines a sensing location where a sensor is arranged based on the calculated significance probability.
JP2015005385A 2015-01-14 2015-01-14 Method, program and apparatus for determining sensing location of sewer network Expired - Fee Related JP6503745B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2015005385A JP6503745B2 (en) 2015-01-14 2015-01-14 Method, program and apparatus for determining sensing location of sewer network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015005385A JP6503745B2 (en) 2015-01-14 2015-01-14 Method, program and apparatus for determining sensing location of sewer network

Publications (2)

Publication Number Publication Date
JP2016130435A true JP2016130435A (en) 2016-07-21
JP6503745B2 JP6503745B2 (en) 2019-04-24

Family

ID=56415828

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015005385A Expired - Fee Related JP6503745B2 (en) 2015-01-14 2015-01-14 Method, program and apparatus for determining sensing location of sewer network

Country Status (1)

Country Link
JP (1) JP6503745B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2016132533A1 (en) * 2015-02-20 2017-12-28 富士通株式会社 Thermoelectric conversion module, sensor module, and information processing system
JP7389526B1 (en) * 2023-07-12 2023-11-30 日本インフラ計測株式会社 Optimization method for water level gauge placement, method for predicting manhole water level, and manhole water level prediction system

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05265513A (en) * 1992-03-19 1993-10-15 Hitachi Ltd Method for controlling, number of pumps in operation at place of drainage pump and device therefor
JPH0762719A (en) * 1993-08-27 1995-03-07 Toshiba Corp Measuring system
US20020156590A1 (en) * 2001-04-20 2002-10-24 Schutzbach James S. Flow transport analysis method and system
JP2003150656A (en) * 2001-11-12 2003-05-23 Toshiba Corp Sewage pipe design supporting device and program
JP2007280161A (en) * 2006-04-10 2007-10-25 Yamatake Corp Inundation occurrence estimating device and program
JP2010133191A (en) * 2008-12-08 2010-06-17 Metawater Co Ltd Anti-inundation measure system of sewage facility
JP2010203964A (en) * 2009-03-04 2010-09-16 Toshiba Corp Monitoring control system of sewerage facility
JP2014182018A (en) * 2013-03-19 2014-09-29 Fujitsu Ltd Correction method of measurement data, measurement data monitoring system, and measurement data monitoring program

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05265513A (en) * 1992-03-19 1993-10-15 Hitachi Ltd Method for controlling, number of pumps in operation at place of drainage pump and device therefor
JPH0762719A (en) * 1993-08-27 1995-03-07 Toshiba Corp Measuring system
US20020156590A1 (en) * 2001-04-20 2002-10-24 Schutzbach James S. Flow transport analysis method and system
JP2003150656A (en) * 2001-11-12 2003-05-23 Toshiba Corp Sewage pipe design supporting device and program
JP2007280161A (en) * 2006-04-10 2007-10-25 Yamatake Corp Inundation occurrence estimating device and program
JP2010133191A (en) * 2008-12-08 2010-06-17 Metawater Co Ltd Anti-inundation measure system of sewage facility
JP2010203964A (en) * 2009-03-04 2010-09-16 Toshiba Corp Monitoring control system of sewerage facility
JP2014182018A (en) * 2013-03-19 2014-09-29 Fujitsu Ltd Correction method of measurement data, measurement data monitoring system, and measurement data monitoring program

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2016132533A1 (en) * 2015-02-20 2017-12-28 富士通株式会社 Thermoelectric conversion module, sensor module, and information processing system
JP7389526B1 (en) * 2023-07-12 2023-11-30 日本インフラ計測株式会社 Optimization method for water level gauge placement, method for predicting manhole water level, and manhole water level prediction system
JP7453719B1 (en) 2023-07-12 2024-03-21 日本インフラ計測株式会社 Optimization method for water level gauge placement, manhole water level prediction method, and manhole water level prediction system

Also Published As

Publication number Publication date
JP6503745B2 (en) 2019-04-24

Similar Documents

Publication Publication Date Title
JP6207889B2 (en) Inundation prediction system, inundation prediction method and program
Liguori et al. Using probabilistic radar rainfall nowcasts and NWP forecasts for flow prediction in urban catchments
Aronica et al. Probabilistic evaluation of flood hazard in urban areas using Monte Carlo simulation
Nie et al. Impacts of climate change on urban drainage systems–a case study in Fredrikstad, Norway
Russo et al. Analysis of extreme flooding events through a calibrated 1D/2D coupled model: the case of Barcelona (Spain)
Yang et al. Linking a storm water management model to a novel two-dimensional model for urban pluvial flood modeling
Gallegos et al. Structural damage prediction in a high-velocity urban dam-break flood: Field-scale assessment of predictive skill
Hussain et al. Modeling the effects of land-use and climate change on the performance of stormwater sewer system using SWMM simulation: case study
Sriwastava et al. Quantifying uncertainty in simulation of sewer overflow volume
JP6794899B2 (en) Liquid level measurement system, liquid level measurement method, and liquid level measurement program
Wu et al. A novel approach for determining integrated water discharge from the ground surface to trunk sewer networks for fast prediction of urban floods
Li et al. Effects of nonstationarity in urban land cover and rainfall on historical flooding intensity in a semiarid catchment
Li Exploring the potential of utilizing unsupervised machine learning for urban drainage sensor placement under future rainfall uncertainty
JP2007280161A (en) Inundation occurrence estimating device and program
JP6503745B2 (en) Method, program and apparatus for determining sensing location of sewer network
Sowby et al. A practical statistical method to differentiate inflow and infiltration in sanitary sewer systems
Wang et al. Analysis of drainage efficiency under extreme precipitation events based on numerical simulation
Nasrin et al. Modelling impact of extreme rainfall on sanitary sewer system by predicting rainfall derived infiltration/inflow
Bean et al. Predicting Utah ground snow loads with PRISM
Sumer et al. Real-time detection of sanitary sewer overflows using neural networks and time series analysis
JP5270322B2 (en) Sewerage facility inundation countermeasure system
Shrestha Impact of climate change on urban flooding in Sukhumvit area of Bangkok
Mugume Modelling and resilience-based evaluation of urban drainage and flood management systems for future cities
Freni et al. Role of modeling uncertainty in the estimation of climate and socioeconomic impact on river water quality
JP2013231313A5 (en)

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20171215

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180925

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20181016

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20181030

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: 20190226

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190311

R150 Certificate of patent or registration of utility model

Ref document number: 6503745

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees