JP6503745B2 - Method, program and apparatus for determining sensing location of sewer network - Google Patents

Method, program and apparatus for determining sensing location of sewer network Download PDF

Info

Publication number
JP6503745B2
JP6503745B2 JP2015005385A JP2015005385A JP6503745B2 JP 6503745 B2 JP6503745 B2 JP 6503745B2 JP 2015005385 A JP2015005385 A JP 2015005385A JP 2015005385 A JP2015005385 A JP 2015005385A JP 6503745 B2 JP6503745 B2 JP 6503745B2
Authority
JP
Japan
Prior art keywords
time
location
series data
sewer network
data
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.)
Expired - Fee Related
Application number
JP2015005385A
Other languages
Japanese (ja)
Other versions
JP2016130435A (en
Inventor
祐輔 樋田
祐輔 樋田
鈴木 貴志
貴志 鈴木
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

Description

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

都市下水道は、降雨時に下水管路に流入してくる雨水を河川に排出したり、生活排水等の汚水を処理することを目的としており、下水道網の下流域にはポンプ場や下水処理場が設置されている。これらの下水道設備は降雨量や流入量等を加味して設計されているため、従来においては下水道管渠から溢水することはほとんどなかった。   The purpose of urban sewerage is to discharge rainwater that flows into the sewage pipeline into the river when it rains, and to treat sewage such as domestic drainage, and there are pumping stations and sewage treatment plants in the lower reaches of the sewer network. is set up. Since these sewerage systems are designed in consideration of the amount of rainfall, inflow, etc., in the past, there was almost no overflow from the sewer pipe.

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

特開2010−203964号公報Unexamined-Japanese-Patent No. 2010-203964 特開2006−18763号公報JP, 2006-18763, A

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

なお、特許文献2等には、予め設定した遅れ時間を利用して下水処理場に流入する下水量を予測する技術が開示されているが、本技術は、下水流入量を予測するためのものであって、センサの設置箇所や数を決定するためのものではない。   In addition, although the technique which estimates the amount of sewage which flows in into a sewage treatment plant using the delay time set beforehand is disclosed by patent document 2 grade | etc., This technique is for estimating the amount of sewage inflow. It is not for determining the installation location or number of sensors.

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

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

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

一実施形態に係るセンシング箇所決定装置のハードウェア構成を概略的に示す図である。It is a figure showing roughly the hardware constitutions of the sensing part determination device concerning 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 sewer hole DB, FIG.3 (b) is a figure which shows the data structure of sewer pipeline DB. 地形DBのデータ構造を示す図である。It is a figure which shows the data structure of topography DB. 降雨量DBのデータ構造を示す図である。It is a figure which shows the data structure of rainfall amount DB. センシング箇所決定装置の処理を示すフローチャートである。It is a flowchart which shows a process of a sensing location determination apparatus. 下水道網の概略図である。It is the schematic of a sewer network. 図8(a)は、人孔3の水位の時系列データの一例を示すグラフであり、図8(b)は、人孔9の水位の時系列データの一例を示すグラフであり、図8(c)は、人孔17の水位の時系列データの一例を示すグラフである。8 (a) is a graph showing an example of time series data of the water level of the human hole 3, and FIG. 8 (b) is a graph showing an example of time series data of the water level of the human hole 9, (C) is a graph which shows an example of the time series data of the water level of the person hole 17. FIG. 人孔3と人孔17の相互相関関数(相関係数)と、人孔9と人孔17の相互相関関数(相関係数)とを示すグラフである。It is a graph which shows the cross correlation function (correlation coefficient) of the person hole 3 and the person hole 17, and the cross correlation function (correlation coefficient) of the person hole 9 and the person hole 17. 図6のステップS18の処理を説明するための図である。It is a figure for demonstrating the process of FIG.6 S18. 検定結果表を示す図(その1)である。It is a figure (the 1) which shows a test result table. 図12(a)、図12(b)は、検定結果表を示す図(その2)である。12 (a) and 12 (b) 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 figures for demonstrating a modification.

以下、センシング箇所決定装置の一実施形態について、図1〜図14に基づいて詳細に説明する。本実施形態のセンシング箇所決定装置は、下水道網に含まれる複数の人孔(マンホール)の中から、水位、流量、流速等の下水状況データを検出するセンサを設置するのに適切な人孔(適切なセンシング箇所)を決定するための装置である。   Hereinafter, one embodiment of a sensing location determination device will be described in detail based on FIGS. 1 to 14. The sensing location determination device according to the present embodiment is a human hole suitable for installing a sensor that detects sewage condition data such as water level, flow rate, flow velocity, etc. from a plurality of human holes (manholes) included in the sewer network It is an apparatus 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 the hardware configuration of a sensing location determination apparatus 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 central processing unit (CPU) 90, a read only memory (ROM) 92, a random access memory (RAM) 94, and a storage. And 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 and the like, and the input unit 95 includes a keyboard, a mouse and 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 by the portable storage medium drive 99 from the portable storage medium 91 (sensing location determination Functions as each unit shown in FIG. 2 by executing the program.

具体的には、CPU90がプログラムを実行することで、図2に示す、入力受付部12、第1〜第3算出部としての演算部14、判定部16、及び出力部18として機能する。なお、本実施形態では、判定部16と出力部18とを含んで決定部としての機能が実現される。なお、図2には、HDD96に格納されているデータベース群(下水道人孔DB22、下水道管路DB24、地形DB26、降雨量DB28を含む)も図示されている。   Specifically, when the CPU 90 executes a program, the CPU 90 functions as the input reception unit 12, the arithmetic unit 14 as the first to third calculation units, the determination unit 16, and the output unit 18 shown 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. Note that FIG. 2 also shows a database group (including a sewer personal hole DB 22, a sewer pipeline DB 24, a topography DB 26, and a rainfall amount 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. In addition, the input reception part 12 transmits the said information (N) to the determination part 16, when the user inputs the largest number (maximum number of sensors) N of the sensor installed in a sewer network. Further, at the stage of receiving the input of the maximum number of sensors N, the input reception unit 12 transmits an instruction to start the process of determining the sensing location to the calculation unit 14.

演算部14は、各種データベースを利用して、下水道網の各箇所(センサ設置箇所候補:本実施形態では各人孔)の有意確率を計算し、有意確率を纏めた検定結果表(図11等参照)を作成する。なお、演算部14による処理の詳細については、後述する。   The calculation unit 14 calculates significance probability of each location (sensor installation location candidate: each human hole in this embodiment) of the sewer network using various databases, and a test result table (FIG. 11 etc.) Create a reference). The 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 a sensor installation location (sensing location) of an appropriate position and number has been determined based on the test result table created by the computation unit 14. As a result of the determination, when the sensing position of the appropriate position and number is determined, the determination unit 16 notifies the output unit 18 of that. On the other hand, when the sensing location of the appropriate position and number is not determined, the determination unit 16 notifies the calculation unit 14 to that effect. In addition, when receiving the notification from the determination unit 16, the calculation unit 14 reduces the sensor installation location candidates, and then calculates the significance probability of each of the sensor installation location candidates again to create a test result table.

出力部18は、表示部93等を介して、適切な位置及び数のセンサ設置箇所(センシング箇所)の情報を出力する。   The output unit 18 outputs the information of the sensor installation location (sensing location) of the appropriate position and number 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 based on FIGS.

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

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

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

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

なお、図3〜図5の各種データベースのデータ構造は、一例である。すなわち、各種データベースが有するフィールドの内容については、種々変更可能である。また、必要に応じて、他のデータベースを追加等してもよい。   In addition, the data structure of the various databases of FIGS. 3-5 is an example. That is, the contents of the fields of the various databases can be variously changed. In addition, other databases may be added as needed.

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

図6の処理では、まず、ステップS10において、入力受付部12が、最大センサ個数Nの入力を受け付ける。この場合、ユーザが入力部95を介して最大センサ個数Nを入力すると、入力受付部12は、当該最大センサ個数Nを受け付け、判定部16に送信する。また、入力受付部12は、最大センサ個数Nの入力を受け付けた段階で、演算部14に対して、処理を開始する旨の指示を出す。なお、本実施形態では、ユーザは最大センサ個数Nとして「5」を入力したものとする。   In the process of FIG. 6, first, in step S10, the input receiving unit 12 receives an input of the maximum number of sensors N. In this case, when the user inputs the maximum number of sensors N via the input unit 95, the input reception unit 12 receives the maximum number of sensors N and transmits it to the determination unit 16. Further, at the stage of receiving the input of the maximum number of sensors N, the input reception unit 12 instructs the calculation unit 14 to start processing. In the present embodiment, it is assumed that the user inputs “5” as the maximum number of sensors 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 the information on the sewer pipe, the topography, and the rainfall amount. Specifically, based on each data of the sewer human-hole DB 22, the sewer pipeline DB 24, the topography DB 26, and the rainfall amount DB 28, the calculation unit 14 calculates each part (each human hole) of the sewer network by a known calculation method. Calculate time series data such as water level, flow rate, and flow velocity. In the present embodiment, as an example, it is assumed that the calculation unit 14 calculates time series data of the water level [m]. The time-series data of the water level in the human hole 3, the human hole 9, and the human hole 17 of FIG. 7 are shown in FIGS. 8A to 8C as an example.

次いで、ステップS14では、演算部14が、算出した時系列データを用いて各人孔と最下流の人孔(ここでは、人孔17とする)の相互相関関数R(τ)を計算する。なお、τは、各箇所(人孔)の時系列データと最下流の箇所(人孔17)の時系列データの変化の時間差を意味する。図9には、人孔3と人孔17の相互相関関数(相関係数)と、人孔9と人孔17の相互相関関数(相関係数)とが示されている。より詳細には、図9には、人孔3や人孔9の時系列データをタイムシフトさせた場合の、人孔17の時系列データとの相関の変化が示されている。   Next, in step S14, the calculation unit 14 calculates the cross-correlation function R (τ) of each human hole and the most downstream human hole (here, the human hole 17) using the calculated time-series data. In addition, (tau) means the time difference of a change of the time-series data of each location (human hole) and the time-series data of the most downstream location (human hole 17). In FIG. 9, 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 are shown. 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 is 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 local maximum value of the cross correlation function. That is, the calculation unit 14 obtains the smallest time difference τ in which the correlation coefficient becomes positive and the second derivative (d 2 R (τ) / dτ 2 ) of the cross-correlation function of the time difference becomes 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 part of FIG. 10, the calculation unit 14 shifts the time series data of the human hole 3 to the right by the delay time T3 (see the broken line graph). In addition, as shown in the middle part of FIG. 10, the calculation unit 14 shifts the time series data of the human hole 9 to the right by the delay time T9 (see the broken line graph). The time shift is similarly performed for time series data of other human holes.

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

次いで、ステップS22では、演算部14が、センサ設置箇所候補(人孔)の検定結果表を作成する。図11には、検定結果表の一例が示されている。図11の検定結果表においては、センサ設置箇所候補(ここでは、人孔1〜16)のt値、P(Probability)−値が纏められている。   Next, in step S22, the calculation unit 14 creates a test result table of sensor installation location candidates (human holes). An example of the test result table is shown in FIG. 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 collected.

次いで、ステップ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 all the significant probabilities are significant. In this case, for example, the determination unit 16 determines that all significant probabilities are significant based on whether all P-values are equal to or less than a predetermined threshold (for example, 0.05). Determine if it is or not. If the determination here 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 calculation unit 14 performs a process of subtracting the sensor installation location candidate. In this case, the operation unit 14 excludes the human hole having the largest P-value exceeding the threshold value (the human hole 2 surrounded by a thick broken line frame in FIG. 11) in step S28 from the sensor installation location candidate. Return to

ステップS20に戻ると、演算部14は、残りのセンサ設置箇所候補(ここでは、人孔1、人孔3〜16)の有意確率を再度計算する。そして、次のステップS22では、演算部14が、残りのセンサ設置箇所候補(人孔)の検定結果表を再度作成する。この場合、図12(a)に示すように、人孔1、人孔3〜16に関する検定結果表が作成される。次いで、ステップS24では、判定部16が、全ての有意確率が有意であるか否かを再度判断する。このステップS24の判断が否定されると、ステップS28に移行する。なお、ステップS28において図12(a)の場合、演算部14は、太破線枠で囲んだ人孔10を除外し、ステップS20に戻る。   If it returns to step S20, the calculating part 14 will calculate again the significant probability of the remaining sensor installation location candidate (here human hole 1, human holes 3-16). Then, in the next step S22, the computing unit 14 creates again the verification result table of the remaining sensor installation location candidates (human holes). In this case, as shown in FIG. 12 (a), a verification result table regarding the human hole 1 and the human holes 3 to 16 is created. Next, in step S24, the determination unit 16 again determines whether all the 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 affirmed, that is, if all the significant probabilities are significant, the determination unit 16 proceeds to step S26. If it transfers to step S26, the determination part 16 will determine whether a sensor installation location candidate is N or less. If the determination here is denied, the process proceeds to step S28, and the person 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 process / determination in steps S20 to S28 is repeatedly performed until the determination in step S26 is affirmed. In addition, the example of the test result table | surface in the step in which the judgment of step S26 is affirmed is shown by FIG.12 (b). 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 less than or equal to N (= 5).

ステップ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 sets the display unit 93 or the like as the optimum installation location of the sensor, with the human holes (human holes 6, 9, 11, 13, 16) included in the verification plan table of FIG. Output through. In this case, the output unit 18 may output the optimum installation location (see black circle in FIG. 13) of the sensor 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 the present embodiment, the operation unit 14 is based on various databases (22, 24, 26, 28), and time series data of water level in a plurality of locations (human holes) in the sewer network Were calculated (S12), and the delay time τ with respect to each of the time-series data in the other human holes of the time-series data in the specific location of the sewer network (the most downstream man hole 17) was calculated (S14, S16) Using the delay time τ and time series data in each of the other human holes, the statistical probability test is performed to calculate the significance probability of each of the other human holes (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 the present 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 according to the number of installed sensors in the present embodiment and the comparative example (example in which the statistical hypothesis test is performed without considering the delay time). In addition, a standard error means the difference | error at the time of estimating the water level of the most downstream manhole 17 using the installed sensor. Comparing the present embodiment with the comparative example, it can be confirmed that the model accuracy is generally improved (the standard error is smaller) in the present embodiment than in the comparative example. Further, in the comparative example, the standard error rapidly increases when the number of sensors is reduced, but it is understood that the standard error gradually increases in the present embodiment. From these, it is understood that in the present 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 the present embodiment, as shown in FIG. 12A and FIG. 13, even if the number of sensing points is reduced to the number N (= 5) or less of the maximum number of sensors input by the user, It is possible to estimate Therefore, in the present 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 situation of the sewer network, and to reduce the sensor installation cost. It becomes.

また、本実施形態では、遅れ時間τを算出する際、最下流以外の人孔17における時系列データをタイムシフトさせた場合に、最下流の人孔17における時系列データとの相関が最も大きくなるシフト時間を遅れ時間として算出する。これにより、遅れ時間τとして適切な値を算出することができる。   Further, in the present embodiment, when the delay time τ is calculated, 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 largest. Shift time is calculated as the delay time. Thus, 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 said embodiment demonstrated the case where the calculating 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. 15 (a) to 15 (c), data similar to time series data of water level (see FIGS. 8 (a) to 8 (c)) can be obtained. The determination unit 16 may determine the optimum installation location of the sensor by the same process as the above embodiment. In addition, the calculation unit 14 may calculate time series data of the flow velocity in step S12. In this case, the calculation unit 14 obtains a value (time for sewage to flow between the human holes) obtained by dividing each distance between the human holes by the flow velocity of the sewage between the human holes, and integrates the value to the most downstream human hole. Therefore, the delay time may be calculated. In the above embodiment, a plurality of time series data of water level, flow rate, and flow velocity are calculated, and the processing of FIG. 6 is executed using each time series data, and the optimum installation of the sensor determined from each time series data From the location, the final installation location of the sensor may be determined.

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

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

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

プログラムを実行するコンピュータは、例えば、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、自己の記憶装置に格納する。そして、コンピュータは、自己の記憶装置からプログラムを読み取り、プログラムに従った処理を実行する。なお、コンピュータは、可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することもできる。また、コンピュータは、サーバコンピュータからプログラムが転送されるごとに、逐次、受け取ったプログラムに従った処理を実行することもできる。   The computer executing 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 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 in accordance with the program. The computer can also execute processing in accordance with the received program each time the program is transferred from the server computer.

上述した実施形態は本発明の好適な実施の例である。但し、これに限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変形実施可能である。   The embodiments described above are examples of preferred implementations of the 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) At least one time-series data of water level, flow rate, and flow velocity at a plurality of locations of the sewer network is calculated based on the data on the sewer network and the data on the rainfall,
Calculating a delay time of each of the time-series data at a specific location of the sewer network with respect to each of the time-series data at a location other than the specific location,
Using the calculated delay time and the time-series data at each location other than the specific location, the statistical probability test is performed to calculate the significance probability of each location other than the specific location,
Determining a sensing location where a sensor is to be placed based on the calculated significance probability
A method for determining a sensing location of a sewer network, wherein a computer executes the processing.
(Supplementary Note 2) In the process of determining, the sensing location is determined based on the condition to be satisfied by the calculated significant probability and the maximum number of sensing locations set in advance. The method of determining the sensing location of the sewer network described in 1.
(Supplementary Note 3) 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 at which the correlation with the time-series data at the specific location is maximized It calculates as a delay time, The sensing location determination method of the sewer network as described in Additional remark 1 or 2 characterized by the above-mentioned.
(Supplementary Note 4) In the process of calculating the time-series data, in the case of calculating time-series data of the flow velocity at a plurality of locations of the sewer network,
In the processing for 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 places. How to determine sensing location.
(Supplementary Note 5) At least one time-series data of water level, flow rate, and flow velocity at a plurality of locations of the sewer network is calculated based on the data on the sewer network and the data on the rainfall,
Calculating a delay time of each of the time-series data at a specific location of the sewer network with respect to each of the time-series data at a location other than the specific location,
Using the calculated delay time and the time-series data at each location other than the specific location, the statistical probability test is performed to calculate the significance probability of each location other than the specific location,
Determining a sensing location where a sensor is to be placed based on the calculated significance probability
The sensing location determination program of the sewer network characterized by making a computer perform a process.
(Supplementary Note 6) In the process of determining, the sensing location is determined based on the condition to be satisfied by the calculated significant probability and the maximum number of sensing locations set in advance. The sensing location determination program for the sewerage network described in 5.
(Supplementary Note 7) In the process of calculating the delay time, the shift time at which the correlation with the time-series data at the specific location is maximized when the time-series data at the location other than the specific location is time-shifted It calculates as a delay time, The sensing location determination program of the sewer network as described in Additional remark 5 or 6 characterized by the above-mentioned.
(Supplementary Note 8) In the process of calculating the time-series data, in the case of calculating time-series data of the flow velocity at a plurality of locations of the sewer network,
In the processing for calculating the delay time, the delay time is calculated based on time series data of the flow velocity and a distance between the plurality of places. Sensing location determination program.
(Supplementary Note 9) A first calculation unit that calculates at least one time-series data of water level, flow rate, and flow velocity at a plurality of locations of the sewer network based on data on the sewer network and data on rainfall.
A second calculator configured to calculate a delay time of each of the time-series data at a specific location of the sewer network with respect to each of the time-series data at a location other than the specific location;
A third calculator configured to calculate a significance probability of each place other than the specific place by a statistical hypothesis test using the calculated delay time and the time-series data in each place other than the specific place;
A determination unit of a sewer network, comprising: a determination unit configured to determine a sensing location at which a sensor is arranged based on the calculated significance probability.
(Supplementary Note 10) The determination part determines the sensing location based on a condition that the calculated significance probability should satisfy and the maximum number of sensing locations set in advance. The sensing location determination device of the sewer network as described in.
(Supplementary Note 11) When the second calculation unit time-shifts the time-series data at a place other than the specific place, a shift time at which the correlation with the time-series data at the specific place becomes largest is a delay time It calculates as, The sensing location determination apparatus of the sewer network as described in Additional remark 9 or 10 characterized by the above-mentioned.
(Supplementary Note 12) When the first calculation unit calculates time-series data of flow velocity at a plurality of locations of the sewerage network,
The sensing location of the sewer network according to Supplementary Note 9 or 10, wherein the second calculation unit calculates a delay time based on time series data of the flow velocity and a distance between the plurality of locations. Decision device.

10 センシング箇所決定装置
14 演算部(第1〜第3算出部)
16 判定部(決定部の一部)
18 出力部(決定部の一部)
90 コンピュータ
10 Sensing location determination device 14 Arithmetic unit (first to third calculation units)
16 Judgment part (part of decision part)
18 Output unit (part of the decision unit)
90 computer

Claims (6)

下水道網に関するデータと降雨量に関するデータに基づき、前記下水道網の複数の箇所における水位、流量、流速の少なくとも1つの時系列データを算出し、
前記下水道網の特定箇所における前記時系列データの、該特定箇所以外の箇所における前記時系列データそれぞれに対する遅れ時間を算出し、
算出した前記遅れ時間と、前記特定箇所以外の箇所それぞれにおける前記時系列データとを用いて、統計的仮説検定により前記特定箇所以外の箇所それぞれの有意確率を算出し、
算出された前記有意確率に基づいてセンサを配置するセンシング箇所を決定する、
処理をコンピュータが実行することを特徴とする下水道網のセンシング箇所決定方法。
Calculating at least one time series data of water level, flow rate and flow velocity at a plurality of locations of the sewer network based on data on the sewer network and data on rainfall amount,
Calculating a delay time of each of the time-series data at a specific location of the sewer network with respect to each of the time-series data at a location other than the specific location,
Using the calculated delay time and the time-series data at each location other than the specific location, the statistical probability test is performed to calculate the significance probability of each location other than the specific location,
Determining a sensing location where a sensor is to be placed based on the calculated significance probability
A method for determining a sensing location of a sewer network, wherein a computer executes the processing.
前記決定する処理では前記有意確率が満たすべき条件と、予め設定されたセンシング箇所の最大数と、に基づいて、前記センシング箇所を決定する、ことを特徴とする請求項1に記載の下水道網のセンシング箇所決定方法。 In the process of the determination, the condition should the significance probability satisfies the maximum number of preset sensing portion, based on the sewer network of claim 1, wherein determining the sensing portion, characterized in that How to determine the location of 前記遅れ時間を算出する処理では、前記特定箇所以外の箇所における前記時系列データをタイムシフトさせた場合に、前記特定箇所における前記時系列データとの相関が最も大きくなるシフト時間を遅れ時間として算出する、ことを特徴とする請求項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, the shift time at which the correlation with the time-series data at the specific location is maximized is calculated as the delay time The sensing location determination method of the sewer network according to claim 1 or 2, characterized in that. 前記時系列データを算出する処理において、前記下水道網の複数の箇所における流速の時系列データを算出する場合に、
前記遅れ時間を算出する処理では、前記流速の時系列データと、前記複数の箇所相互間の距離とに基づいて、遅れ時間を算出することを特徴とする請求項1又は2に記載の下水道網のセンシング箇所決定方法。
In the process of calculating the time-series data, in the case of calculating time-series data of 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 time series data of the flow velocity and a distance between the plurality of places. How to determine the location of
下水道網に関するデータと降雨量に関するデータに基づき、前記下水道網の複数の箇所における水位、流量、流速の少なくとも1つの時系列データを算出し、
前記下水道網の特定箇所における前記時系列データの、該特定箇所以外の箇所における前記時系列データそれぞれに対する遅れ時間を算出し、
算出した前記遅れ時間と、前記特定箇所以外の箇所それぞれにおける前記時系列データとを用いて、統計的仮説検定により前記特定箇所以外の箇所それぞれの有意確率を算出し、
算出された前記有意確率に基づいてセンサを配置するセンシング箇所を決定する、
処理をコンピュータに実行させることを特徴とする下水道網のセンシング箇所決定プログラム。
Calculating at least one time series data of water level, flow rate and flow velocity at a plurality of locations of the sewer network based on data on the sewer network and data on rainfall amount,
Calculating a delay time of each of the time-series data at a specific location of the sewer network with respect to each of the time-series data at a location other than the specific location,
Using the calculated delay time and the time-series data at each location other than the specific location, the statistical probability test is performed to calculate the significance probability of each location other than the specific location,
Determining a sensing location where a sensor is to be placed based on the calculated significance probability
The sensing location determination program of the sewer network characterized by making a computer perform a process.
下水道網に関するデータと降雨量に関するデータに基づき、前記下水道網の複数の箇所における水位、流量、流速の少なくとも1つの時系列データを算出する第1算出部と、
前記下水道網の特定箇所における前記時系列データの、該特定箇所以外の箇所における前記時系列データそれぞれに対する遅れ時間を算出する第2算出部と、
算出した前記遅れ時間と、前記特定箇所以外の箇所それぞれにおける前記時系列データとを用いて、統計的仮説検定により前記特定箇所以外の箇所それぞれの有意確率を算出する第3算出部と、
算出された前記有意確率に基づいてセンサを配置するセンシング箇所を決定する決定部と、を備える下水道網のセンシング箇所決定装置。
A first calculation unit that calculates at least one time-series data of water level, flow rate, and velocity at a plurality of locations of the sewer network based on data on the sewer network and data on rainfall amount;
A second calculator configured to calculate a delay time of each of the time-series data at a specific location of the sewer network with respect to each of the time-series data at a location other than the specific location;
A third calculator configured to calculate a significance probability of each place other than the specific place by a statistical hypothesis test using the calculated delay time and the time-series data in each place other than the specific place;
A determination unit of a sewer network, comprising: a determination unit configured to determine a sensing location at which 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 JP2016130435A (en) 2016-07-21
JP6503745B2 true 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)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016132533A1 (en) * 2015-02-20 2016-08-25 富士通株式会社 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

Family Cites Families (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
US6757623B2 (en) * 2001-04-20 2004-06-29 Ads Corporation Flow transport analysis method and system
JP2003150656A (en) * 2001-11-12 2003-05-23 Toshiba Corp Sewage pipe design supporting device and program
JP4828991B2 (en) * 2006-04-10 2011-11-30 株式会社山武 Overflow generation estimation device and program
JP5270322B2 (en) * 2008-12-08 2013-08-21 メタウォーター株式会社 Sewerage facility inundation countermeasure system
JP2010203964A (en) * 2009-03-04 2010-09-16 Toshiba Corp Monitoring control system of sewerage facility
JP6244637B2 (en) * 2013-03-19 2017-12-13 富士通株式会社 Measurement data correction method, measurement data monitoring system, and measurement data monitoring program

Also Published As

Publication number Publication date
JP2016130435A (en) 2016-07-21

Similar Documents

Publication Publication Date Title
Nie et al. Impacts of climate change on urban drainage systems–a case study in Fredrikstad, Norway
JP6207889B2 (en) Inundation prediction system, inundation prediction method and program
Aronica et al. Probabilistic evaluation of flood hazard in urban areas using Monte Carlo simulation
Mancipe-Munoz et al. Calibration of rainfall-runoff model in urban watersheds for stormwater management assessment
Yang et al. Linking a storm water management model to a novel two-dimensional model for urban pluvial flood modeling
Hussain et al. Modeling the effects of land-use and climate change on the performance of stormwater sewer system using SWMM simulation: case study
CN115471078B (en) Flood risk point assessment method and device based on urban water affair system
Zhou et al. Adaptation to urbanization impacts on drainage in the city of Hohhot, China
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
Pereira et al. Improving operational management of wastewater systems. A case study
JP6503745B2 (en) Method, program and apparatus for determining sensing location of sewer network
JP2007280161A (en) Inundation occurrence estimating device and program
Yap et al. Assessment on inflow and infiltration in sewerage systems of Kuantan, Pahang
Zeng et al. A new probability-embodied model for simulating variable contributing areas and hydrologic processes dominated by surface depressions
Sowby et al. A practical statistical method to differentiate inflow and infiltration in sanitary sewer systems
JP3839361B2 (en) Rainwater runoff coefficient prediction method, rainwater inflow prediction method, rainwater runoff coefficient prediction program, and rainwater inflow forecast program
CN113282577A (en) Sewage pipe network monitoring method and device, electronic equipment and storage medium
Chen et al. Modelling the effect of rainfall patterns on the runoff control performance of permeable pavements
Nasrin et al. Modelling impact of extreme rainfall on sanitary sewer system by predicting rainfall derived infiltration/inflow
Wang et al. Analysis of drainage efficiency under extreme precipitation events based on numerical simulation
Bean et al. Predicting Utah ground snow loads with PRISM
Wang et al. Spotting strategic storm drain inlets in flat urban catchments
Wang et al. Uncertainty analysis of a pollutant-hydrograph model in assessing inflow and infiltration of sanitary sewer systems
JP5925585B2 (en) Main line manhole emergency simulation system, method and program
JP5270322B2 (en) Sewerage facility inundation countermeasure system

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