JP2013182448A - Maximum supply amount prediction system, and maximum supply amount prediction method - Google Patents

Maximum supply amount prediction system, and maximum supply amount prediction method Download PDF

Info

Publication number
JP2013182448A
JP2013182448A JP2012046159A JP2012046159A JP2013182448A JP 2013182448 A JP2013182448 A JP 2013182448A JP 2012046159 A JP2012046159 A JP 2012046159A JP 2012046159 A JP2012046159 A JP 2012046159A JP 2013182448 A JP2013182448 A JP 2013182448A
Authority
JP
Japan
Prior art keywords
supply amount
prediction
unit
period
maximum
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.)
Pending
Application number
JP2012046159A
Other languages
Japanese (ja)
Inventor
Shinnosuke Kimura
新之介 木村
Takafumi Aoki
崇文 青木
Kazuyoshi Takahashi
一喜 高橋
Toshinori Sasaya
俊徳 笹谷
Hidetaka Shinozaki
英孝 篠崎
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.)
Tokyo Gas Co Ltd
Original Assignee
Tokyo Gas Co 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 Tokyo Gas Co Ltd filed Critical Tokyo Gas Co Ltd
Priority to JP2012046159A priority Critical patent/JP2013182448A/en
Publication of JP2013182448A publication Critical patent/JP2013182448A/en
Pending legal-status Critical Current

Links

Images

Landscapes

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

Abstract

PROBLEM TO BE SOLVED: To provide a maximum supply amount prediction system capable of predicting a maximum supply amount of energy with high accuracy.SOLUTION: A temperature calculating part 22 acquires past supply amount record data, past average temperature record data, scheduled supply amount data for each intended purpose, and the like on the basis of a prediction unit and a prediction period received by an input reception part 21. On the basis of the acquired data, the temperature calculating part 22 calculates a temperature for the prediction unit and the prediction period; a supply amount calculating part 23 calculates a supply amount in the prediction unit and the prediction period for each month, for each day of the week, and for each intended purpose, an in-period maximum supply amount calculating part 24 calculates a supply amount in the prediction unit and the prediction period. A predicted maximum supply amount calculating part 25 calculates a maximum supply amount which can be generated at a constant possibility on the basis of a simulation result of supply amount distribution in the prediction unit and the prediction period.

Description

本発明は、エネルギーの最大供給量予測システム等に関するものである。   The present invention relates to a maximum energy supply amount prediction system and the like.

従来、ガスや電気等のエネルギー供給企業では、エネルギーセキュリティを高める上で必要となる設備稼働計画や設備投資計画のために、精度良く、エネルギーの最大供給量を予測しておく必要がある。   Conventionally, an energy supply company such as gas or electricity needs to accurately predict the maximum energy supply amount for an equipment operation plan and an equipment investment plan which are necessary for enhancing energy security.

例えば、特許文献1には、熟練オペレータの需要予測技術を複数定式化し、実績との相違で成績を付けて、予測精度を高める技術が提案されている。   For example, Patent Literature 1 proposes a technique for formulating a plurality of demand prediction techniques for skilled operators, giving a result based on a difference from the actual results, and improving the prediction accuracy.

特開2007−241360号公報JP 2007-241360 A

しかしながら、特許文献1の技術では、過去の実績データ及び気象情報が考慮されていることに留まる。エネルギーの最大供給量を精度良く算出するためには、例えば、エネルギーの用途ごとに異なる「気温に対する感度」等を考慮する必要がある。   However, in the technique of Patent Document 1, past performance data and weather information are considered. In order to calculate the maximum supply amount of energy with high accuracy, for example, it is necessary to consider “sensitivity to air temperature” and the like that differ for each use of energy.

また、特許文献1の技術では、特定の1日における需要の予測を行うことに留まる。設備稼働計画や設備投資計画のために有用な情報は、ある期間(例えば、ガス需要量のピーク期である「12月〜2月」)におけるエネルギーの最大供給量である。   Moreover, with the technique of patent document 1, it will remain only in the prediction of the demand in a specific day. Information useful for the facility operation plan and the facility investment plan is the maximum energy supply amount in a certain period (for example, “December to February”, which is the peak period of gas demand).

本発明は、前述した問題点に鑑みてなされたもので、その目的とすることは、設備稼働計画や設備投資計画のために有用なエネルギーの最大供給量を精度良く算出することが可能な最大供給量予測システム等を提供することである。   The present invention has been made in view of the above-described problems, and the purpose of the present invention is to maximize the maximum supply amount of energy useful for an equipment operation plan and an equipment investment plan. It is to provide a supply amount prediction system and the like.

前述した目的を達成するために、第1の発明は、コンピュータによって構成され、エネルギーの最大供給量を予測する最大供給量予測システムであって、前記最大供給量の予測を行う予測単位及び予測期間の入力を受け付ける入力受付手段と、前記入力受付手段により受け付けられた前記予測単位及び予測期間に基づく、過去の供給量実績データを取得する取得手段と、前記予測単位及び予測期間における気温を算出する第1の算出手段と、前記取得手段により取得された前記過去の供給量実績データおよび前記第1の算出手段により算出された前記気温に基づいて、前記予測単位及び予測期間における供給量の分布を算出する第2の算出手段と、前記第2の算出手段により算出された前記予測単位及び予測期間における供給量の分布から、前記最大供給量を算出する第3の算出手段と、を備えることを特徴とする最大供給量予測システムである。第1の発明によって、特定の期間、特定の時間、特定の時刻といった予測単位及び予測期間における最大供給量を精度良く予測することができる。   In order to achieve the above-described object, the first invention is a maximum supply amount prediction system configured by a computer and predicting the maximum supply amount of energy, and a prediction unit and a prediction period for predicting the maximum supply amount Input receiving means for receiving the input, acquisition means for acquiring past supply amount actual data based on the prediction unit and prediction period received by the input receiving means, and calculating the temperature in the prediction unit and prediction period Based on the first calculation means, the past supply amount actual data acquired by the acquisition means, and the temperature calculated by the first calculation means, the distribution of the supply amount in the prediction unit and the prediction period is calculated. From the second calculation means for calculating, and the distribution of the supply amount in the prediction unit and the prediction period calculated by the second calculation means, A third calculating means for calculating the serial maximum supply amount, a maximum supply amount prediction system comprising: a. According to the first invention, the maximum supply amount in the prediction unit and the prediction period such as the specific period, the specific time, and the specific time can be accurately predicted.

前記第2の算出手段は、更に、用途別に、前記予測単位及び予測期間における供給量の分布を算出する。用途別に算出する理由は、例えば、家庭用の用途であれば、気温に対する感度が高いが、工業用の用途であれば、気温に対する感度はあまり高くないと考えられるため、用途別に算出することによって、気温に対する感度の違いを考慮することができるからである。   The second calculation means further calculates the distribution of the supply amount in the prediction unit and the prediction period for each use. The reason for calculating by use is, for example, that for home use, the sensitivity to temperature is high, but for industrial use, it is considered that the sensitivity to temperature is not so high. This is because the difference in sensitivity to temperature can be taken into account.

前記第2の算出手段は、前記予測単位及び予測期間における供給量の分布の最大値のみを記憶部に記憶し、前記第3の算出手段は、前記予測単位及び予測期間における供給量の分布の最大値から前記最大供給量を算出する。これによって、使用メモリを最小限にすることができる。   The second calculation means stores only a maximum value of the supply amount distribution in the prediction unit and the prediction period in the storage unit, and the third calculation means calculates the supply amount distribution in the prediction unit and the prediction period. The maximum supply amount is calculated from the maximum value. As a result, the memory used can be minimized.

前記第1の算出手段は、前記予測単位及び予測期間における気温を、過去の実績値から取得する。これによって、現実に即した計算ができる。   The first calculation means acquires the temperature in the prediction unit and the prediction period from past performance values. This makes it possible to calculate according to reality.

第2の発明は、コンピュータが実行し、エネルギーの最大供給量を予測する最大供給量予測方法であって、前記コンピュータが、前記最大供給量の予測を行う予測単位及び予測期間の入力を受け付ける入力受付ステップと、前記入力受付ステップにより受け付けられた前記予測単位及び予測期間に基づく、過去の供給量実績データを取得する取得ステップと、前記予測単位及び予測期間における気温を算出する第1の算出ステップと、前記取得ステップにより取得された前記過去の供給量実績データおよび前記第1の算出ステップにより算出された前記気温に基づいて、前記予測単位及び予測期間における供給量の分布を算出する第2の算出ステップと、前記第2の算出ステップにより算出された前記予測単位及び予測期間における供給量の分布から、前記最大供給量を算出する第3の算出ステップと、を含むことを特徴とする最大供給量予測方法である。第2の発明によって、特定の期間、特定の時間、特定の時刻といった予測単位及び予測期間における最大供給量を精度良く予測することができる。   A second invention is a maximum supply amount prediction method executed by a computer and predicting a maximum supply amount of energy, wherein the computer receives an input of a prediction unit for performing prediction of the maximum supply amount and a prediction period. A receiving step, an acquisition step of acquiring past supply amount actual data based on the prediction unit and the prediction period received by the input receiving step, and a first calculation step of calculating the temperature in the prediction unit and the prediction period And calculating a distribution of supply amount in the prediction unit and prediction period based on the past supply amount actual data acquired in the acquisition step and the temperature calculated in the first calculation step. The amount of supply in the prediction unit and the prediction period calculated in the calculation step and the second calculation step. From the fabric, the maximum supply amount prediction method characterized by including: a third calculation step of calculating the maximum supply amount. According to the second invention, the maximum supply amount in the prediction unit and the prediction period such as the specific period, the specific time, and the specific time can be accurately predicted.

本発明により、設備稼働計画や設備投資計画のために有用なエネルギーの最大供給量を精度良く予測することが可能な最大供給量予測システムおよび最大供給量予測方法を提供することができる。   According to the present invention, it is possible to provide a maximum supply amount prediction system and a maximum supply amount prediction method capable of accurately predicting the maximum supply amount of energy useful for facility operation plans and facility investment plans.

最大供給量予測システムの概要を示す図である。It is a figure which shows the outline | summary of the maximum supply amount prediction system. 端末(サーバ)を実現するコンピュータのハードウエア構成図である。It is a hardware block diagram of the computer which implement | achieves a terminal (server). 端末内の制御部の機能構成例を示すブロック図である。It is a block diagram which shows the function structural example of the control part in a terminal. 最大供給量予測システム1の最大供給量算出処理を説明するフローチャートである。It is a flowchart explaining the maximum supply amount calculation process of the maximum supply amount prediction system 1. 予測単位及び予測期間における気温を算出する例を説明するための正規分布である。It is a normal distribution for demonstrating the example which calculates the temperature in a prediction unit and a prediction period. 予測単位及び予測期間における気温(平均、標準偏差)の入力値の例を示す図である。It is a figure which shows the example of the input value of the temperature (average, standard deviation) in a prediction unit and a prediction period. 月別、曜日別の予測単位及び予測期間における供給量の分布を算出する処理の流れを説明するためのデータフローである。It is a data flow for demonstrating the flow of the process which calculates distribution of the supply amount in the prediction unit according to a month and a day of the week, and a prediction period. 用途毎にt℃における供給量の分布のグラフデータの出力例である。It is an example of an output of graph data of distribution of supply amount in t ° C for every use. 全用途のt℃における供給量の分布のグラフデータの出力例である。It is an output example of the graph data of the distribution of the supply amount at t ° C. for all uses. 予測単位及び予測期間における予測最大供給量を算出する例を説明するためのグラフである。It is a graph for demonstrating the example which calculates the prediction maximum supply amount in a prediction unit and a prediction period.

以下、図面に基づいて、本発明の実施形態を詳細に説明する。以下では、最大供給量予測の対象エネルギーとして、ガスを例に挙げて説明する。但し、本発明の適用範囲はガスに限らない。本発明は、前述した課題があるエネルギーに対して適用すれば、同様の効果を奏する。   Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings. Hereinafter, gas will be described as an example of target energy for maximum supply amount prediction. However, the application range of the present invention is not limited to gas. The present invention has the same effect when applied to energy having the above-described problems.

図1は、最大供給量予測システム1の概要を示す図である。図1に示すように、最大供給量予測システム1は、例えば、端末2とサーバ3がネットワーク6を介して接続されている。ネットワーク6は、例えば、インターネット又はLAN(Local Area Network)等である。端末2は、例えば、PC(Personal
Computer、以下「コンピュータ」)や携帯端末(携帯電話、スマートフォン、タブレット端末等)等であり、ネットワーク6に接続し、データ通信(HTTP通信、TCP/IP通信など)が可能であれば、どのような機器でも良い。サーバ3も、端末2と同様に、ネットワーク6に接続し、データ通信が可能であれば良いが、望ましくは高性能なサーバ用コンピュータが良い。
FIG. 1 is a diagram showing an outline of the maximum supply amount prediction system 1. As shown in FIG. 1, in the maximum supply amount prediction system 1, for example, a terminal 2 and a server 3 are connected via a network 6. The network 6 is, for example, the Internet or a LAN (Local Area Network). The terminal 2 is, for example, a PC (Personal
Computer, hereinafter referred to as “computer”) or a mobile terminal (cell phone, smartphone, tablet terminal, etc.), etc. Simple equipment may be used. Similarly to the terminal 2, the server 3 may be connected to the network 6 and can perform data communication, but is preferably a high-performance server computer.

端末2には、本発明の一形態である供給量予測プログラム4がインストールされている。本発明の実施形態では、端末2が、供給量予測プログラム4に従って各種の手段として機能し、予測されるガスの最大需要量に対して必要なガス供給量を算出し、ユーザに提示する。端末2は、必要に応じて、データの要求命令をサーバ3に送信する。   A supply amount prediction program 4 that is an embodiment of the present invention is installed in the terminal 2. In the embodiment of the present invention, the terminal 2 functions as various means according to the supply amount prediction program 4, calculates a necessary gas supply amount for the predicted maximum demand amount of gas, and presents it to the user. The terminal 2 transmits a data request command to the server 3 as necessary.

サーバ3には、本発明の実施形態において利用される各種のデータを記憶するデータベース(以下「DB」)5が構築されている。データベース5には、例えば、過去の用途別供給量データ、過去の供給量実績データ、過去の平均気温実績データ、用途別計画供給量データ等が記憶されている。計画供給量データとは、過去のデータに基づいて統計的に算出されたデータである。サーバ3は、データの要求に対してDB5を検索し、要求されたデータを端末2に送信する。   In the server 3, a database (hereinafter “DB”) 5 for storing various data used in the embodiment of the present invention is constructed. The database 5 stores, for example, past supply data by usage, past supply actual data, past average temperature actual data, usage-specific planned supply data, and the like. The planned supply amount data is data statistically calculated based on past data. The server 3 searches the DB 5 for the data request and transmits the requested data to the terminal 2.

尚、最大供給量予測システム1の構成は、図1に示す例に限らない。例えば、最大供給量予測システム1は、端末2のみで構成されても良い。つまり、端末2が、DB5を備えるようにしても良い。   The configuration of the maximum supply amount prediction system 1 is not limited to the example shown in FIG. For example, the maximum supply amount prediction system 1 may be configured by only the terminal 2. That is, the terminal 2 may include the DB 5.

また、供給量予測プログラム4は、サーバ3にインストールされていても良い。つまり、サーバ3が、供給量予測プログラム4に従って各種の手段として機能し、予測されるガスの最大需要量に対して必要なガス供給量を算出し、ユーザに提示するようにしても良い。この場合、端末2は、ユーザとのインターフェースの役割を果たす。つまり、端末2は、ユーザから入力されるデータをサーバ3に送信し、サーバ3から受信するデータを出力(表示や印刷など)する。   The supply amount prediction program 4 may be installed in the server 3. That is, the server 3 may function as various means according to the supply amount prediction program 4 to calculate a necessary gas supply amount for the predicted maximum gas demand and present it to the user. In this case, the terminal 2 serves as an interface with the user. That is, the terminal 2 transmits data input from the user to the server 3 and outputs (displays, prints, etc.) data received from the server 3.

また、DB5に代えて、単なるファイルとしてデータを記憶しても良い。また、DB5に記憶されているデータは、外部のサーバから取得するようにしても良い。   Further, instead of DB5, data may be stored as a simple file. The data stored in the DB 5 may be acquired from an external server.

図2は、端末2(サーバ3)を実現するコンピュータのハードウエア構成図である。尚、図2のハードウエア構成は一例であり、用途、目的に応じて様々な構成を採ることが可能である。   FIG. 2 is a hardware configuration diagram of a computer that realizes the terminal 2 (server 3). Note that the hardware configuration in FIG. 2 is an example, and various configurations can be adopted depending on the application and purpose.

端末2(サーバ3)を実現するコンピュータは、制御部11、記憶部12、メディア入出力部13、通信制御部14、入力部15、表示部16、周辺機器I/F部17等が、バス18を介して接続される。   A computer that realizes the terminal 2 (server 3) includes a control unit 11, a storage unit 12, a media input / output unit 13, a communication control unit 14, an input unit 15, a display unit 16, a peripheral device I / F unit 17, and the like. 18 is connected.

制御部11は、CPU(Central Processing Unit)、ROM(Read Only Memory)、RAM(Random Access Memory)等で構成される。CPUは、記憶部12、ROM、記録媒体等に格納されるプログラムをRAM上のワークメモリ領域に呼び出して実行し、バス18を介して接続された各装置を駆動制御し、コンピュータが行う後述する処理を実現する。ROMは、不揮発性メモリであり、コンピュータのブートプログラムやBIOS等のプログラム、データ等を恒久的に保持している。RAMは、揮発性メモリであり、記憶部12、ROM、記録媒体等からロードしたプログラム、データ等を一時的に保持するとともに、制御部11が各種処理を行う為に使用するワークエリアを備える。   The control unit 11 includes a CPU (Central Processing Unit), a ROM (Read Only Memory), a RAM (Random Access Memory), and the like. The CPU calls a program stored in the storage unit 12, ROM, recording medium, etc. to a work memory area on the RAM and executes it, drives and controls each device connected via the bus 18, and will be described later. Realize processing. The ROM is a non-volatile memory and permanently holds a computer boot program, a program such as BIOS, data, and the like. The RAM is a volatile memory, and temporarily stores programs, data, and the like loaded from the storage unit 12, ROM, recording medium, and the like, and includes a work area used by the control unit 11 for performing various processes.

記憶部12は、例えば、HDD(Hard Disk Drive)であり、制御部11が実行するプログラム、プログラム実行に必要なデータ、OS(Operating System)等が格納される。プログラムに関しては、OSに相当する制御プログラムや、後述する処理をコンピュータに実行させるためのアプリケーションプログラム(供給量予測プログラム4)が格納されている。これらの各プログラムコードは、制御部11により必要に応じて読み出されてRAMに移され、CPUに読み出されて各種の手段として実行される。   The storage unit 12 is, for example, an HDD (Hard Disk Drive), and stores a program executed by the control unit 11, data necessary for program execution, an OS (Operating System), and the like. As for the program, a control program corresponding to the OS and an application program (supply amount prediction program 4) for causing a computer to execute processing to be described later are stored. Each of these program codes is read by the control unit 11 as necessary, transferred to the RAM, read by the CPU, and executed as various means.

メディア入出力部13(ドライブ装置)は、データの入出力を行い、例えば、CDドライブ(−ROM、−R、−RW等)、DVDドライブ(−ROM、−R、−RW等)等のメディア入出力装置を有する。通信制御部14は、通信制御装置、通信ポート等を有し、コンピュータとネットワーク6間の通信を媒介する通信インターフェースであり、ネットワーク6を介して、他のコンピュータ間との通信制御を行う。ネットワーク6は、有線、無線を問わない。   The media input / output unit 13 (drive device) inputs / outputs data, for example, media such as a CD drive (-ROM, -R, -RW, etc.), DVD drive (-ROM, -R, -RW, etc.) Has input / output devices. The communication control unit 14 includes a communication control device, a communication port, and the like, and is a communication interface that mediates communication between the computer and the network 6, and controls communication with other computers via the network 6. The network 6 may be wired or wireless.

入力部15は、データの入力を行い、例えば、キーボード、マウス等のポインティングデバイス、テンキー等の入力装置を有する。入力部15を介して、コンピュータに対して、操作指示、動作指示、データ入力等を行うことができる。表示部16は、液晶パネル、有機EL等のディスプレイ装置、ディスプレイ装置と連携してコンピュータのビデオ機能を実現するための論理回路等(ビデオアダプタ等)を有する。尚、入力部15及び表示部16は、タッチパネルディスプレイのように、一体となっていても良い。   The input unit 15 inputs data and includes, for example, a keyboard, a pointing device such as a mouse, and an input device such as a numeric keypad. An operation instruction, an operation instruction, data input, and the like can be performed on the computer via the input unit 15. The display unit 16 includes a liquid crystal panel, a display device such as an organic EL, and a logic circuit or the like (video adapter or the like) for realizing a video function of a computer in cooperation with the display device. The input unit 15 and the display unit 16 may be integrated like a touch panel display.

周辺機器I/F(インターフェース)部17は、コンピュータに周辺機器を接続させるためのポートであり、周辺機器I/F部17を介してコンピュータは周辺機器とのデータの送受信を行う。周辺機器I/F部17は、USB等で構成されており、通常複数の周辺機器I/Fを有する。周辺機器との接続形態は有線、無線を問わない。バス18は、各装置間の制御信号、データ信号等の授受を媒介する経路である。   The peripheral device I / F (interface) unit 17 is a port for connecting a peripheral device to the computer, and the computer transmits and receives data to and from the peripheral device via the peripheral device I / F unit 17. The peripheral device I / F unit 17 is configured by a USB or the like, and usually includes a plurality of peripheral devices I / F. The connection form with the peripheral device may be wired or wireless. The bus 18 is a path that mediates transmission / reception of control signals, data signals, and the like between the devices.

図3は、端末2内の制御部11の機能構成例を示すブロック図である。図3に示す機能部のうちの少なくとも一部は、図2の制御部11により供給量予測プログラム4が実行されることによって実現される。   FIG. 3 is a block diagram illustrating a functional configuration example of the control unit 11 in the terminal 2. At least a part of the functional units shown in FIG. 3 is realized by the supply amount prediction program 4 being executed by the control unit 11 of FIG.

入力受付部21は、表示部16に表示されている入力画面(図示せず)において、利用者が入力部15を操作することで、最大供給量を求めたい単位(日、時刻等)の入力を受け付けるとともに、最大供給量を予測する期間の入力を受け付ける。以下、最大供給量を求めたい単位を「予測単位」と定義し、最大供給量を予測する期間を「予測期間」と定義する。   The input receiving unit 21 inputs a unit (day, time, etc.) for which a maximum supply amount is to be obtained by a user operating the input unit 15 on an input screen (not shown) displayed on the display unit 16. And an input of a period for predicting the maximum supply amount. Hereinafter, a unit for which the maximum supply amount is desired is defined as “prediction unit”, and a period during which the maximum supply amount is predicted is defined as “prediction period”.

入力受付部21は、利用者から受け付けた入力情報(予測単位及び予測期間)を、気温算出部22に供給する。入力受付部21は、予測単位と予測期間を独立に受け付ける。例えば、予測期間が「12月〜2月」、予測単位が「日」である場合は、12月〜2月の中での最大日供給量を求めることになる。また、予測期間が「12月」、予測単位が「8〜9時」である場合は、12月の中での8〜9時の最大供給量を求めことになる。後述する処理に用いられる過去の用途別供給量実績データ及び過去の供給量実績データについては、予測単位に合わせて、入力データの単位を変える。   The input reception unit 21 supplies the input information (prediction unit and prediction period) received from the user to the temperature calculation unit 22. The input receiving unit 21 receives the prediction unit and the prediction period independently. For example, when the prediction period is “December to February” and the prediction unit is “day”, the maximum daily supply amount in December to February is obtained. When the prediction period is “December” and the prediction unit is “8 to 9 o'clock”, the maximum supply amount at 8 to 9 o'clock in December is obtained. For past actual supply amount actual data and past supply amount actual data used for processing to be described later, the unit of input data is changed in accordance with the prediction unit.

気温算出部22は、ネットワーク6を介してサーバ3のDB5にアクセスし、入力受付部21から供給された入力情報(予測単位及び予測期間)に基づく、過去の用途別供給量実績データ、過去の供給量実績データ、過去の平均気温実績データ、用途別計画供給量データ等を取得する。気温算出部22は、予測単位及び予測期間における過去の平均気温実績データから、予測単位及び予測期間における気温を算出し、その算出結果を、予測単位及び予測期間における過去の用途別供給量実績データ、過去の供給量実績データ、用途別計画供給量データとともに供給量算出部23に供給する。   The air temperature calculation unit 22 accesses the DB 5 of the server 3 via the network 6, and based on the input information (prediction unit and prediction period) supplied from the input reception unit 21, the past supply amount actual data by use, the past Acquire actual supply amount data, past average temperature actual result data, usage-specific planned supply amount data, and the like. The temperature calculation unit 22 calculates the temperature in the prediction unit and the prediction period from the past average temperature actual data in the prediction unit and the prediction period, and the calculation result is used as the past supply amount actual data for each usage in the prediction unit and the prediction period. The supply amount is supplied to the supply amount calculation unit 23 together with past supply amount actual data and usage-specific planned supply amount data.

供給量算出部23は、気温算出部22から供給された予測単位及び予測期間における気温、過去の用途別供給量実績データ、過去の供給量実績データ、用途別の計画供給量データから、月別、曜日別、用途別の予測単位における供給量の分布を算出する。そして、供給量算出部23は、用途間の相関を考慮して全用途のデータを合算し、用途合算された単一の供給量を算出する。   The supply amount calculation unit 23 calculates the monthly temperature from the prediction unit and the temperature in the prediction period supplied from the temperature calculation unit 22, the past supply amount actual data by use, the past supply amount actual data, and the planned supply amount data by use, Calculate the distribution of supply in forecast units by day of the week and usage. Then, the supply amount calculation unit 23 adds the data for all uses in consideration of the correlation between uses, and calculates a single supply amount that is added up.

気温算出部22および供給量算出部23は、予測単位及び予測期間における供給量の算出処理を繰り返すことで、予測単位及び予測期間における全用途の供給量を求めることができ、その算出結果を、期間中の最大供給量算出部24に出力する。   The temperature calculation unit 22 and the supply amount calculation unit 23 can obtain supply amounts for all uses in the prediction unit and the prediction period by repeating the calculation process of the supply amount in the prediction unit and the prediction period, Output to the maximum supply amount calculation unit 24 during the period.

期間中の最大供給量算出部24は、供給量算出部23から供給された予測単位及び予測期間における全用途の供給量の算出処理を、シミュレーションにより、さらに所定回数繰り返し、その最大値を算出する。期間中の最大供給量算出部24は、予測単位及び予測期間における供給量の分布のシミュレーション結果を、予測最大供給量算出部25に出力する。   The maximum supply amount calculation unit 24 during the period repeats the prediction unit supplied from the supply amount calculation unit 23 and the supply amount of all uses in the prediction period by a predetermined number of times, and calculates the maximum value. . The maximum supply amount calculation unit 24 during the period outputs a simulation result of the prediction unit and the distribution of the supply amount in the prediction period to the predicted maximum supply amount calculation unit 25.

予測最大供給量算出部25は、期間中の最大供給量算出部24から供給された、予測単位及び予測期間における最大供給量のシミュレーション結果から、多数の最大供給量を算出し、さらにその分布を算出する。予測最大供給量算出部25は、多数の最大供給量の分布から、一定の確率で発生しうる最大供給量を算出する。   The predicted maximum supply amount calculation unit 25 calculates a number of maximum supply amounts from the simulation result of the maximum supply amount in the prediction unit and the prediction period supplied from the maximum supply amount calculation unit 24 during the period, and further distributes the distribution. calculate. The predicted maximum supply amount calculation unit 25 calculates the maximum supply amount that can occur with a certain probability from the distribution of a large number of maximum supply amounts.

次に、図4のフローチャートを参照して、最大供給量予測システム1の最大供給量算出処理について説明する。   Next, the maximum supply amount calculation process of the maximum supply amount prediction system 1 will be described with reference to the flowchart of FIG.

ステップS1において、入力受付部21は、表示部16に表示されている入力画面(図示せず)において、利用者から、予測単位及び予測期間の入力を受け付ける。   In step S <b> 1, the input receiving unit 21 receives an input of a prediction unit and a prediction period from a user on an input screen (not shown) displayed on the display unit 16.

例えば、予測期間が「12月〜2月」、予測単位が「日(例えば、12月1日、12月2日、・・・)」である。   For example, the prediction period is “December to February”, and the prediction unit is “day (for example, December 1, December 2,...)”.

なお、予測単位や予測期間の他に、「予測地域」というパラメータを用いることも可能である。予測地域は、予測対象の地域を限定して予測する場合に、利用者により指定される。これによって、予測対象の地域における予測最大供給量を算出することができ、ガスの導管設備の健全性を検証することもできる。   In addition to the prediction unit and the prediction period, a parameter called “prediction region” can also be used. The prediction area is specified by the user when the prediction target area is limited. Thereby, the predicted maximum supply amount in the prediction target area can be calculated, and the soundness of the gas conduit facility can be verified.

ステップS2において、気温算出部22は、ネットワーク6を介してサーバ3のDB5にアクセスし、ステップS1の処理で受け付けられた予測単位及び予測期間における、過去の用途別供給量実績データ、過去の供給量実績データ、過去の平均気温実績データ、用途別計画供給量データ等を取得する。そして気温算出部22は、取得した過去の平均気温実績データから、予測単位及び予測期間における気温を算出する。   In step S2, the air temperature calculation unit 22 accesses the DB 5 of the server 3 via the network 6, and the past usage amount actual data and past supply in the prediction unit and the prediction period received in the process of step S1. Acquire actual volume data, past average temperature actual data, usage-specific planned supply volume data, etc. And the temperature calculation part 22 calculates the temperature in a prediction unit and a prediction period from the acquired past average temperature actual result data.

図5は、予測単位及び予測期間における気温を算出する例を説明するための正規分布である。図5において、横軸は気温であり、縦軸は発生確率である。この例では、予測単位が「日」、予測期間が「12月〜2月」として気温を算出する場合について説明する。   FIG. 5 is a normal distribution for explaining an example of calculating the temperature in the prediction unit and the prediction period. In FIG. 5, the horizontal axis is the temperature, and the vertical axis is the occurrence probability. In this example, a case where the temperature is calculated assuming that the prediction unit is “day” and the prediction period is “December to February” will be described.

気温算出部22は、12月〜2月の各月それぞれにおいて、予測単位及び予測期間における気温を、正規分布に従うとしてシミュレーションにより算出する。例えば、正規分布の平均は、当該月の予想平均気温であり、標準偏差は、当該月の実績気温の標準偏差の過去5年分のデータを用いて算出する。   The temperature calculation unit 22 calculates the temperature in the prediction unit and the prediction period by simulation in accordance with the normal distribution in each month from December to February. For example, the average of the normal distribution is the expected average temperature of the month, and the standard deviation is calculated using data for the past five years of the standard deviation of the actual temperature of the month.

図6は、予測単位及び予測期間における気温の入力の例を示す図である。   FIG. 6 is a diagram illustrating an example of the input of the temperature in the prediction unit and the prediction period.

図6の例では、12月に、計画気温「8.9℃」、標準偏差「2.6℃」を入力し、1月に、計画気温「6.5℃」、標準偏差「1.9℃」を入力し、2月に、計画気温「7.2℃」、標準偏差「2.4℃」を入力している。   In the example of FIG. 6, the planned temperature “8.9 ° C.” and the standard deviation “2.6 ° C.” are input in December, and the planned temperature “6.5 ° C.” and the standard deviation “1.9” are input in January. "Centigrade" was entered, and in February, the planned temperature "7.2 C" and the standard deviation "2.4 C" were entered.

図4の説明に戻る。ステップS3において、供給量算出部23は、ステップS2の処理で算出された予測単位及び予測期間における気温、過去の用途別供給量実績データ、過去の供給量実績データ、用途別の計画供給量データから、月別、曜日別、用途別の予測単位及び予測期間における供給量の分布を算出する。   Returning to the description of FIG. In step S3, the supply amount calculation unit 23 calculates the temperature in the prediction unit and the prediction period calculated in the process of step S2, past supply amount actual data by use, past supply amount actual data, and planned supply amount data by use. Then, the distribution of the supply amount in the prediction unit and prediction period for each month, each day of the week, and each use is calculated.

例えば、「曜日別」とは、「月、金」や「火、水、木」等の曜日単位で区分することであり、「用途別」とは、「A用途、B用途、C用途、・・・」等の用途単位で区分することである。なお、曜日別において、「月、金」と「火、水、木」とに区分する理由は、一般に、月曜日や金曜日等において供給量が下がる傾向にあるためである。   For example, “by day of the week” means to classify by day of the week such as “Monday, Friday” or “Tue, Wed, Thu”, etc., and “by use” means “A use, B use, C use, "..." and so on. In addition, the reason for dividing into “Monday, Friday” and “Tue, Wed, Thu” by day of the week is that the supply amount generally tends to decrease on Mondays, Fridays, and the like.

図7は、月別、曜日別の予測単位における供給量の分布を算出する処理の流れを説明するためのデータフローである。   FIG. 7 is a data flow for explaining the flow of processing for calculating the distribution of the supply amount in the prediction unit by month and day of the week.

供給量算出部23は、曜日31、気温t℃32、用途別計画供給量33、過去の用途別供給量実績データ34、過去の供給量実績データ(全用途合算データ)35、過去の平均気温実績データ36を用いて、月別、曜日別、用途別に、気温と供給量の関係を示すグラフデータ37を作成し、出力する。   The supply amount calculation unit 23 includes a day of the week 31, an air temperature t ° C. 32, a planned supply amount 33 according to usage, a past supply amount actual data 34 according to past usage, a past supply amount actual data (all usage combined data) 35, and a past average temperature. Using the record data 36, graph data 37 indicating the relationship between the temperature and the supply amount is created and output for each month, each day of the week, and each use.

ここで、過去の用途別供給量実績データ34の求め方について、具体例を挙げて説明する。例えば、過去の供給量実績データ(全用途合算データ)35は有しているが、過去の用途別供給量実績データ34は有していない場合を考える。この場合、用途別の需要者に対してサンプリング調査を行い、サンプリング調査の結果及び過去の供給量実績データ(全用途合算データ)35から、過去の用途別供給量実績データ34を推測する。例えば、以下のように、12月1日のサンプリング調査データの合計値が得られたとする。
用途A(サンプリング値):10
用途B(サンプリング値):20
用途C(サンプリング値):30
用途全体合計値X(サンプリング値):60
Here, a method for obtaining the past supply amount data 34 by usage will be described with a specific example. For example, consider a case where the past supply amount actual data (total usage total data) 35 is included, but the past supply amount actual data 34 is not included. In this case, a sampling survey is performed with respect to the consumer for each usage, and the past supply volume actual data 34 for each usage is estimated from the result of the sampling survey and the past supply volume actual data (total usage total data) 35. For example, it is assumed that the total value of sampling survey data on December 1 is obtained as follows.
Application A (sampling value): 10
Application B (sampling value): 20
Application C (sampling value): 30
Total use total value X (sampling value): 60

これに対して、サンプリング調査データの合計値が、母集団全体のデータとなるように、各合計値を膨らませる。例えば、サンプル数:母集団数=1:10の場合、以下のようになる。
用途A(サンプリング値の単純増幅値):100
用途B(サンプリング値の単純増幅値):200
用途C(サンプリング値の単純増幅値):300
用途全体合計値X(サンプリング値の単純増幅値):600
On the other hand, each total value is expanded so that the total value of the sampling survey data becomes the data of the entire population. For example, when the number of samples: the number of populations = 1: 10, it is as follows.
Application A (simple amplification value of sampling value): 100
Application B (simple amplified value of sampling value): 200
Use C (simple amplified value of sampling value): 300
Total use total value X (simple amplification value of sampling value): 600

さらに、用途全体合計値X(サンプリング値の単純増幅値)が、過去の供給量実績データ(全用途合算データ)35と等しくなるように、用途A(サンプリング値の単純増幅値)、用途B(サンプリング値の単純増幅値)、用途C(サンプリング値の単純増幅値)を補正し、推測値を求める。例えば、過去の供給量実績データ(全用途合算データ)35が660の場合、以下のようになる。
用途A(推測値):110
用途B(推測値):220
用途C(推測値):330
用途全体合計値X(実績値):660
Furthermore, application A (simple amplification value of sampling value), application B (simple amplification value of sampling value) so that total use value X (simple amplification value of sampling value) becomes equal to past supply amount actual data (total application total data) 35. The estimated value is obtained by correcting the simple amplification value of the sampling value) and the application C (simple amplification value of the sampling value). For example, when the past supply amount actual data (total usage total data) 35 is 660, it is as follows.
Application A (estimated value): 110
Application B (estimated value): 220
Application C (estimated value): 330
Total use value X (actual value): 660

このように、一部をサンプリング調査のデータに頼る場合であっても、実績値に基づいてサンプリング調査のデータを補正することによって、実態に即した予測が可能となる。   As described above, even when a part of the data is based on the sampling survey data, the actual survey can be predicted by correcting the sampling survey data based on the actual value.

供給量算出部23は、以上のような算出結果を、月別、曜日別、気温別にまとめることによって、月別、曜日別、用途別の気温と供給量の分布の関係を示すグラフデータ37を作成し、出力することができる。   The supply amount calculation unit 23 creates graph data 37 indicating the relationship between the distribution of the temperature and the supply amount by month, by day of the week, and by use by summarizing the above calculation results by month, day of the week, and temperature. Can be output.

用途別に算出する理由は、例えば、家庭用の用途であれば、気温に対する感度が高いが、工業用の用途であれば、気温に対する感度はあまり高くないと考えられるため、用途別に算出することによって、気温に対する感度の違いを考慮することができるからである。   The reason for calculating by use is, for example, that for home use, the sensitivity to temperature is high, but for industrial use, it is considered that the sensitivity to temperature is not so high. This is because the difference in sensitivity to temperature can be taken into account.

尚、例えば供給量が日供給量の場合、日供給量を月間実績供給量で除算するようにしても良い。理由は、サンプル数を増やすために複数年度のデータを使用しており、年度毎のばらつきを正規化するためである。なおその場合はグラフデータ38を算出する際、月間値の用途別計画供給量33をかける必要がある。またこのことにより、当該予測期間における景気変動等による過去との供給量全体のベース量の差を考慮することもできる。   For example, when the supply amount is the daily supply amount, the daily supply amount may be divided by the monthly actual supply amount. The reason is to use data for a plurality of years in order to increase the number of samples, and to normalize the variation for each year. In this case, when calculating the graph data 38, it is necessary to multiply the monthly planned supply amount 33 by usage. In addition, this makes it possible to take into account the difference in the base amount of the entire supply amount from the past due to economic fluctuations in the forecast period.

例えば、用途毎に気温と供給量の分布の関係を示すグラフ(グラフデータ37)を作成すると、気温が低いほど供給量が多くなる(気温感度がある)用途と、気温に依らず供給量が一様に分布している(気温感度がない)用途が存在することが分かる。   For example, if a graph (graph data 37) showing the relationship between the temperature and the supply amount distribution is created for each application, the supply amount increases as the temperature decreases (the temperature is sensitive), and the supply amount does not depend on the temperature. It can be seen that there are applications that are uniformly distributed (no temperature sensitivity).

図8は、用途毎にt℃における供給量の分布38のグラフデータ39の出力例である。図8において、横軸は供給量であり、縦軸は頻度である。   FIG. 8 is an output example of the graph data 39 of the supply amount distribution 38 at t ° C. for each application. In FIG. 8, the horizontal axis is the supply amount, and the vertical axis is the frequency.

図8(a)は、用途A、気温5℃における供給量の分布を示すグラフである。図8(b)は、用途B、気温5℃における供給量の分布を示すグラフである。   FIG. 8A is a graph showing the distribution of the supply amount at application A and at an air temperature of 5 ° C. FIG. 8B is a graph showing the distribution of the supply amount at application B and at an air temperature of 5 ° C.

供給量算出部23は、算出した用途毎のt℃における供給量の分布38から、正規分布を仮定して用途別の相関を考慮しながら、用途合算の供給量の分布40を算出し、グラフデータ41を出力する。   The supply amount calculation unit 23 calculates the distribution 40 of the combined supply amount from the calculated supply amount distribution 38 at t ° C. for each application, assuming a normal distribution and considering the correlation for each application. Data 41 is output.

用途別の相関を考慮する場合、月間計画供給量Si,mは、次式(1)に従って算出される。S´i,mは、補正前の月間計画供給量であり、n´mは、m月の日数である。 When considering the correlation for each application, the monthly planned supply amount S i, m is calculated according to the following equation (1). S ′ i, m is the monthly planned supply amount before correction, and n ′ m is the number of days in m month.

Figure 2013182448
Figure 2013182448

上記式(1)に示すように、月間計画供給量Si,mは、月の日数で補正したm月の用途iの月間計画供給量(m)となる。 As shown in the above equation (1), the monthly planned supply amount S i, m is the monthly planned supply amount (m 3 ) of the use i for m months corrected by the number of days of the month.

このとき、合算日供給量の分布Ym,n(m月、曜日区分nの全用途の日供給量合計の分布)は、正規分布N(μY,t,m,σY,t,m)に従う分布となる。平均μY,t,mは、次式(2)に示す通り定義され、標準偏差σY,t,mは、次式(3)に示す通り定義される。なお、Iは、気温感度がある用途の集合、Jは、気温感度がない用途の集合、tは、m月の日平均気温の分布である。ただし、係数ai,n、係数bi,n、共分散σesn,etnは、過去の実績データから予め算出しておく。 At this time, the distribution Y m, n of the total daily supply amount (the distribution of the total daily supply amount for all uses in m month and day of week classification n) is the normal distribution N (μ Y, t, m , σ Y, t, m ). The average μ Y, t, m is defined as shown in the following equation (2), and the standard deviation σ Y, t, m is defined as shown in the following equation (3). Here, I is a set of uses with temperature sensitivity, J is a set of uses without temperature sensitivity, and t m is a distribution of daily average temperatures in m months. However, the coefficient a i, n , the coefficient b i, n , and the covariance σ esn, etn are calculated in advance from past performance data.

Figure 2013182448
Figure 2013182448

供給量算出部23は、以上のような数式を用いることによって、正規分布を仮定して、用途別の相関を考慮しながら、用途合算の供給量の分布40を算出し、グラフデータ41を出力することができる。そしてその分布からランダムに1つの供給量を算出する。   The supply amount calculation unit 23 assumes the normal distribution by using the above formula, calculates the distribution 40 of the combined use amount while considering the correlation for each use, and outputs the graph data 41. can do. Then, one supply amount is randomly calculated from the distribution.

図9は、全用途のt℃における供給量の分布のグラフデータ41の出力例である。図9において、横軸は供給量であり、縦軸は頻度である。特に、t=5℃における供給量の分布を示すグラフである。   FIG. 9 is an output example of the graph data 41 of the distribution of the supply amount at t ° C. for all uses. In FIG. 9, the horizontal axis is the supply amount, and the vertical axis is the frequency. In particular, it is a graph showing the distribution of the supply amount at t = 5 ° C.

図4の説明に戻る。ステップS4において、供給量算出部23は、予測単位及び予測期間における供給量の分布を算出したか否かを判定し、まだ算出していないと判定した場合、12月〜2月の各月それぞれにおいて、ステップS2、S3の処理を繰り返し実行する。そして、ステップS4において、供給量算出部23は、予測単位及び予測期間における全用途の供給量の分布を算出したと判定した場合、その分布からランダムに1つの供給量を算出しステップS5に進む。   Returning to the description of FIG. In step S4, the supply amount calculation unit 23 determines whether or not the distribution of the supply amount in the prediction unit and the prediction period has been calculated. If it is determined that the distribution has not yet been calculated, each month from December to February is determined. In Steps S2 and S3, the processes are repeatedly executed. In step S4, when the supply amount calculation unit 23 determines that the distribution of supply amounts for all uses in the prediction unit and the prediction period is calculated, one supply amount is randomly calculated from the distribution, and the process proceeds to step S5. .

ステップS5において、期間中の最大供給量算出部24は、ステップS3の処理で算出された全用途のt℃における供給量の分布40からランダムに1つの供給量を算出し、期間中のシミュレーション結果から単純に最大値をとることで、期間中の最大供給量を算出する。例えば、前のフェーズで、12月1日の供給量が「10」、12月2日の供給量が「15」、12月3日の供給量が「12」、・・・の算出結果が得られた場合、最大供給量は、MAX(10、15、12、・・・)によって算出される。   In step S5, the maximum supply amount calculation unit 24 during the period randomly calculates one supply amount from the supply amount distribution 40 at t ° C. of all uses calculated in the process of step S3, and the simulation result during the period. The maximum supply amount during the period is calculated by simply taking the maximum value. For example, in the previous phase, the calculation result of the supply amount on December 1 is “10”, the supply amount on December 2 is “15”, the supply amount on December 3 is “12”,. If obtained, the maximum supply amount is calculated by MAX (10, 15, 12,...).

ステップS6において、期間中の最大供給量算出部24は、期間中の最大供給量を算出する処理を所定回数(例えば、100万回)行ったか否かを判定し、まだ所定回数行っていないと判定した場合、ステップS2に戻り、上述した処理を繰り返し実行する。つまり、期間中の最大供給量算出部24は、ステップS2〜S6の処理を所定回数繰り返し実行(シミュレーション)することにより、予測単位及び予測期間における最大供給量の分布を算出する。   In step S6, the maximum supply amount calculation unit 24 during the period determines whether or not the process of calculating the maximum supply amount during the period has been performed a predetermined number of times (for example, 1 million times). When it determines, it returns to step S2 and performs the process mentioned above repeatedly. That is, the maximum supply amount calculation unit 24 during the period calculates the distribution of the maximum supply amount in the prediction unit and the prediction period by repeatedly executing (simulating) the processes of steps S2 to S6 a predetermined number of times.

そして、ステップS6において、期間中の最大供給量算出部24は、期間中の最大供給量を算出する処理を所定回数行ったと判定した場合、ステップS7に進む。   In step S6, if the maximum supply amount calculation unit 24 during the period determines that the process of calculating the maximum supply amount during the period has been performed a predetermined number of times, the process proceeds to step S7.

ステップS7において、予測最大供給量算出部25は、ステップS2〜S6の処理を所定回数繰り返し実行(シミュレーション)することにより、予測単位及び予測期間における多数の最大供給量の分布を算出する。そして、予測最大供給量算出部25は、多数の最大供給量の分布から、一定の確率で発生しうる最大供給量を算出する。   In step S <b> 7, the predicted maximum supply amount calculation unit 25 repeatedly calculates (simulates) the processes in steps S <b> 2 to S <b> 6 a predetermined number of times, thereby calculating a distribution of a large number of maximum supply amounts in the prediction unit and the prediction period. Then, the predicted maximum supply amount calculation unit 25 calculates the maximum supply amount that can occur with a certain probability from the distribution of a large number of maximum supply amounts.

図10は、予測単位及び予測期間における予測最大供給量を算出する例を説明するためのグラフである。図10において、横軸は予測最大供給量であり、縦軸は発生確率である。   FIG. 10 is a graph for explaining an example of calculating the predicted maximum supply amount in the prediction unit and the prediction period. In FIG. 10, the horizontal axis represents the predicted maximum supply amount, and the vertical axis represents the occurrence probability.

図10(a)は、予測単位及び予測期間における予測最大供給量の分布のシミュレーション結果(ステップS7の処理による結果)を示すグラフである。   FIG. 10A is a graph showing a simulation result of the distribution of the predicted maximum supply amount in the prediction unit and the prediction period (result of the process in step S7).

このような、予測単位及び予測期間における予測最大供給量を算出するシミュレーションを所定回数(無限回)繰り返した結果は、矢印A1の先に示す図10(b)のグラフのようになる。   The result of repeating such a simulation for calculating the predicted maximum supply amount in the prediction unit and the prediction period a predetermined number of times (infinite times) is as shown in the graph of FIG.

予測最大供給量算出部25は、図10(b)に示すグラフから、一定の確率(例えば、90%の信頼区間)で発生しうる最大供給量を算出する。予測最大供給量算出部25は、算出結果(予測結果)を表示部16に表示させるとともに、記憶部12に記憶させる。   The predicted maximum supply amount calculation unit 25 calculates the maximum supply amount that can occur with a certain probability (for example, a 90% confidence interval) from the graph shown in FIG. The predicted maximum supply amount calculation unit 25 displays the calculation result (prediction result) on the display unit 16 and also stores it in the storage unit 12.

このように、ある期間における供給量の最大値を予測することが可能となる。   Thus, it becomes possible to predict the maximum value of the supply amount in a certain period.

また、シミュレーションによる100万回の繰り返し計算において、各日(例えば、12月1日、12月2日、・・・等)の予測結果を全て記憶部12に記憶するのではなく、12月〜2月の予測結果の最大値のみを記憶部12に記憶することによって、メモリ使用量を少なくできる。   Further, in the repeated calculation of 1 million times by simulation, not all the prediction results of each day (for example, December 1, December 2,...) Are stored in the storage unit 12 but from December to By storing only the maximum value of the prediction result for February in the storage unit 12, the memory usage can be reduced.

また、前述の説明では、該当日の気温が正規分布に従うと仮定してシミュレーションするものとしたが、本発明はこの例に限定されない。例えば、正規分布以外の分布を用いても良い。また、例えば、予測単位が「8〜9時」の場合、「8〜9時」の気温を用いる等、予測単位、予測期間に相関が強い時間帯の気温を用いても良い。また、例えば、該当日の気温は過去の実績値からランダムに取得してシミュレーションすることで、現実的には起こり得ない気温(例えば、−50°など)を考慮することがなく、現実に即した計算ができる。   In the above description, the simulation is performed assuming that the temperature of the day follows a normal distribution, but the present invention is not limited to this example. For example, a distribution other than the normal distribution may be used. Further, for example, when the prediction unit is “8 to 9 o'clock”, the temperature in the time zone having a strong correlation between the prediction unit and the prediction period may be used, such as using the temperature of “8 to 9 o'clock”. In addition, for example, the temperature of the corresponding day is randomly obtained from the past actual values and simulated, so that it is not necessary to take into consideration a temperature that cannot actually occur (for example, −50 ° etc.). Can be calculated.

[発明の実施の形態における効果]
1.利用者によって指定された特定の期間、特定の時間、特定の時刻といった予測単位における最大供給量を精度良く算出(予測)することができる。
2.用途別に、最大供給量の分布を算出することによって、気温に対する感度の違いを考慮することができる。
3.予測単位及び予測期間における最大供給量の分布の算出結果として、最大値のみを記憶部に記憶するため、使用メモリを最小限にすることができる。
4.過去の実績値を用いてシミュレーションすることで、現実に即した計算ができる。
[Effects of the embodiment of the invention]
1. It is possible to accurately calculate (predict) the maximum supply amount in a prediction unit such as a specific period, a specific time, and a specific time designated by the user.
2. By calculating the distribution of maximum supply for each application, the difference in sensitivity to temperature can be taken into account.
3. Since only the maximum value is stored in the storage unit as the calculation result of the distribution of the maximum supply amount in the prediction unit and the prediction period, the memory used can be minimized.
4). By using the past actual values, simulations can be performed according to reality.

[変形例]
以上においては、シミュレーションによって最大供給量を算出するようにしたが、これに限らず、解析的な計算によって最大供給量を算出することも可能である。
[Modification]
In the above description, the maximum supply amount is calculated by simulation. However, the present invention is not limited to this, and the maximum supply amount can also be calculated by analytical calculation.

ここで、解析的な計算による最大供給量の算出方法について説明する。まず、最大供給量の算出処理に用いる計算式の変数について定義する。
i,m,n:m月、曜日区分nの用途iの日供給量の分布(m
m,n:m月、曜日区分nの、全用途の日供給量合計の分布(合算日供給量の分布)(m
y,m,n(x):m月、曜日区分nの、合算日供給量の分布Ym,nの累積確率密度関数
i,m:月の日数で補正したm月の用途iの月間計画供給量(m
S´i,m:月の日数を補正する前のm月の用途iの月間計画供給量(m
I:気温感度がある用途の集合
J:気温感度がない用途の集合
:m月の日平均気温の分布(正規分布N(μt、m、σt、m)に従う分布とする。)
:日最大供給量の分布(m
Y*(x):日最大供給量の分布Yの累積確率密度関数
M:最大供給量が発生する月の集合
N:曜日区分の集合
m,n:m月、曜日区分nの平日の日数
n´m:m月の日数
Here, a method of calculating the maximum supply amount by analytical calculation will be described. First, the variables of the calculation formula used for the calculation process of the maximum supply amount are defined.
y i, m, n : Distribution of daily supply for use i in m month and day of week division n (m 3 )
Y m, n : Distribution of the total daily supply for all uses in m month, day of week classification n (distribution of combined daily supply) (m 3 )
F y, m, n (x): cumulative probability density function Y m, n distribution of total daily supply for m month, day of week division n Si, m : m month usage i corrected by the number of days in the month Monthly planned supply (m 3 )
S ′ i, m : Monthly planned supply amount for use i in m month before correcting the number of days in the month (m 3 )
I: Set of uses with temperature sensitivity J: Set of uses without temperature sensitivity t m : Distribution of daily average temperature of m month (with a distribution according to normal distribution N (μ t, m , σ t, m ))
Y * : Distribution of maximum daily supply (m 3 )
F Y * (x): Cumulative probability density function of daily maximum supply distribution Y * M: Set of months in which maximum supply occurs N: Set of days of week division n m, n : Weekdays of m month and day of week division n Days n ' m : days in m month

まず、m月、曜日区分nの、合算日供給量の分布Ym,nを、次式(4)の通り定義し、分布(平均、標準偏差)を算出する。
ただし、気温感度がある用途では、
i,m,n=Si,m×(ai,n×t+bi,n+ei,n);ei,nは正規分布N(0,σei,m)に従う。
一方、気温感度がない用途では、
j,m,n=Sj,m×(bj,n+ej,n);ej,nは正規分布N(0,σej,n)に従う。
尚、es,nとet,nの共分散をσesn,etn(∀s∈I,J ,∀t∈I,J)とする。
First, the distribution Y m, n of the total day supply amount for m month and day of week classification n is defined as the following equation (4), and the distribution (average, standard deviation) is calculated.
However, in applications with temperature sensitivity,
y i, m, n = S i, m × (a i, n × t + b i, n + e i, n ); e i, n follows a normal distribution N (0, σ ei, m ).
On the other hand, in applications where there is no temperature sensitivity,
y j, m, n = S j, m × (b j, n + e j, n ); e j, n follows a normal distribution N (0, σ ej, n ).
Note that the covariance between e s, n and e t, n is σ esn , etn (∀sεI, J, ∀tεI, J).

Figure 2013182448
Figure 2013182448

このとき、月の日数で補正したm月の用途iの月間計画供給量Si,mは、上記式(1)で示したように、日数補正前の月間計画供給量S´i,mを補正することで求められる。また係数ai,n、係数bi,n、共分散σesn,etnは、過去の実績データを用いて算出される。 At this time, monthly plan supply amount S i applications i of m months corrected in days of the month, m, as shown in the above formula (1), before the days correction monthly plan supply amount S'i, the m It is calculated by correcting. The coefficient a i, n , coefficient b i, n , and covariance σ esn, etn are calculated using past performance data.

そして、m月の日平均気温の分布tの平均μt,m、標準偏差σt,m、日数補正済みのm月の用途iの月間計画供給量Si,m(∀i∈I,J)を与えることにより、合算供給量の分布Ym,nを算出する。 Then, the average μ t, m , standard deviation σ t, m of the daily average temperature distribution t in m month, and the monthly planned supply amount S i, m (∀i∈I, J ) To calculate the distribution Y m, n of the combined supply amount.

以上算出した分布は、ピーク期平日の一日当たりの日供給量の分布である。この分布から、ピーク期間(最大供給量が発生する月の集合M)での最大供給量の分布を算出する。m月、曜日区分nの合算日供給量の分布の累積確率密度関数Fy,m,n(x)、およびm月、曜日区分nの平日の日数nm,nを与えた時の、ピーク期間での最大供給量の分布は、次式(5)の通り定義される。 The distribution calculated above is the distribution of daily supply per day during peak weekdays. From this distribution, the distribution of the maximum supply amount in the peak period (a set M of months in which the maximum supply amount occurs) is calculated. Peak when the cumulative probability density function F y, m, n (x) of the distribution of the total daily supply of m month and day of week division n and the number of weekdays n m, n of m month and day of week division n are given The distribution of the maximum supply amount in the period is defined as the following equation (5).

Figure 2013182448
Figure 2013182448

以上のように、解析的な計算では、例えば、「12月、1月、2月の日平均気温の分布が正規分布に従うこと」を仮定し、最大供給量の分布を算出することで、計算処理時間の短縮を図ることが可能となる。   As described above, in the analytical calculation, for example, assuming that the distribution of daily average temperatures in December, January, and February follows a normal distribution, the distribution of the maximum supply is calculated. It is possible to shorten the processing time.

以上、添付図面を参照しながら、本発明に係る画像出力受付端末等の好適な実施形態について説明したが、本発明はかかる例に限定されない。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。   The preferred embodiments of the image output receiving terminal and the like according to the present invention have been described above with reference to the accompanying drawings, but the present invention is not limited to such examples. It will be apparent to those skilled in the art that various changes or modifications can be conceived within the scope of the technical idea disclosed in the present application, and these naturally belong to the technical scope of the present invention. Understood.

1………最大供給量予測システム
2………端末
4………供給量予測プログラム
11………制御部
12………記憶部
15………入力部
16………表示部
21………入力受付部
22………該当日の気温算出部
23………該当日の日供給量算出部
24………期間中の最大供給量算出部
25………予測最大供給量算出部
1 ......... Maximum supply amount prediction system 2 ......... Terminal 4 ......... Supply amount prediction program 11 ......... Control unit 12 ......... Storage unit 15 ......... Input unit 16 ......... Display unit 21 ......... Input accepting unit 22 ......... Temperature calculation unit for the day 23 ......... Daily supply amount calculation unit for the day 24 ......... Maximum supply amount calculation unit during the period 25 ......... Predicted maximum supply amount calculation unit

Claims (5)

コンピュータによって構成され、エネルギーの最大供給量を予測する最大供給量予測システムであって、
前記最大供給量の予測を行う予測単位及び予測期間の入力を受け付ける入力受付手段と、
前記入力受付手段により受け付けられた前記予測単位及び予測期間に基づく、過去の供給量実績データを取得する取得手段と、
前記予測単位及び予測期間における気温を算出する第1の算出手段と、
前記取得手段により取得された前記過去の供給量実績データおよび前記第1の算出手段により算出された前記気温に基づいて、前記予測単位及び予測期間における供給量の分布を算出する第2の算出手段と、
前記第2の算出手段により算出された前記予測単位及び予測期間における供給量の分布から、前記最大供給量を算出する第3の算出手段と、
を備えることを特徴とする最大供給量予測システム。
A maximum supply amount prediction system configured by a computer and predicting the maximum supply amount of energy,
Input receiving means for receiving an input of a prediction unit and a prediction period for predicting the maximum supply amount;
An acquisition means for acquiring past supply amount actual data based on the prediction unit and the prediction period received by the input reception means;
First calculating means for calculating the temperature in the prediction unit and the prediction period;
Second calculation means for calculating a distribution of the supply amount in the prediction unit and the prediction period based on the past supply amount actual data acquired by the acquisition means and the temperature calculated by the first calculation means. When,
Third calculation means for calculating the maximum supply amount from the prediction unit calculated by the second calculation means and the distribution of supply amount in the prediction period;
A maximum supply amount prediction system comprising:
前記第2の算出手段は、更に、用途別に、前記予測単位及び予測期間における供給量の分布を算出する
ことを特徴とする請求項1に記載の最大供給量予測システム。
2. The maximum supply amount prediction system according to claim 1, wherein the second calculation unit further calculates a supply amount distribution in the prediction unit and the prediction period for each application.
前記第2の算出手段は、前記予測単位及び予測期間における供給量の分布の最大値のみを記憶部に記憶し、
前記第3の算出手段は、前記予測単位及び予測期間における供給量の分布の最大値から前記最大供給量を算出する
ことを特徴とする請求項1又は請求項2に記載の最大供給量予測システム。
The second calculation unit stores only the maximum value of the distribution of the supply amount in the prediction unit and the prediction period in the storage unit,
3. The maximum supply amount prediction system according to claim 1, wherein the third calculation unit calculates the maximum supply amount from the maximum value of the distribution of the supply amount in the prediction unit and the prediction period. .
前記第1の算出手段は、前記予測単位及び予測期間における気温を、過去の実績値から取得する
ことを特徴とする請求項1乃至請求項3のいずれかに記載の最大供給量予測システム。
The said 1st calculation means acquires the temperature in the said prediction unit and prediction period from the past performance value, The maximum supply amount prediction system in any one of the Claims 1 thru | or 3 characterized by the above-mentioned.
コンピュータが実行し、エネルギーの最大供給量を予測する最大供給量予測方法であって、
前記コンピュータが、
前記最大供給量の予測を行う予測単位及び予測期間の入力を受け付ける入力受付ステップと、
前記入力受付ステップにより受け付けられた前記予測単位及び予測期間に基づく、過去の供給量実績データを取得する取得ステップと、
前記予測単位及び予測期間における気温を算出する第1の算出ステップと、
前記取得ステップにより取得された前記過去の供給量実績データおよび前記第1の算出ステップにより算出された前記気温に基づいて、前記予測単位及び予測期間における供給量の分布を算出する第2の算出ステップと、
前記第2の算出ステップにより算出された前記予測単位及び予測期間における供給量の分布から、前記最大供給量を算出する第3の算出ステップと、
を含むことを特徴とする最大供給量予測方法。
A maximum supply amount prediction method executed by a computer and predicting the maximum supply amount of energy,
The computer is
An input receiving step for receiving an input of a prediction unit for performing prediction of the maximum supply amount and a prediction period;
An acquisition step of acquiring past supply amount actual data based on the prediction unit and the prediction period received by the input reception step;
A first calculation step of calculating the temperature in the prediction unit and the prediction period;
Second calculation step for calculating distribution of supply amount in the prediction unit and prediction period based on the past supply amount actual data acquired in the acquisition step and the air temperature calculated in the first calculation step. When,
A third calculation step for calculating the maximum supply amount from the distribution of the supply amount in the prediction unit and the prediction period calculated in the second calculation step;
The maximum supply amount prediction method characterized by including.
JP2012046159A 2012-03-02 2012-03-02 Maximum supply amount prediction system, and maximum supply amount prediction method Pending JP2013182448A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012046159A JP2013182448A (en) 2012-03-02 2012-03-02 Maximum supply amount prediction system, and maximum supply amount prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012046159A JP2013182448A (en) 2012-03-02 2012-03-02 Maximum supply amount prediction system, and maximum supply amount prediction method

Publications (1)

Publication Number Publication Date
JP2013182448A true JP2013182448A (en) 2013-09-12

Family

ID=49273046

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012046159A Pending JP2013182448A (en) 2012-03-02 2012-03-02 Maximum supply amount prediction system, and maximum supply amount prediction method

Country Status (1)

Country Link
JP (1) JP2013182448A (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08110332A (en) * 1994-10-11 1996-04-30 Komatsu Ltd Method and apparatus for ultrasonic inspecting by articulated robot
JP2001318970A (en) * 2000-05-12 2001-11-16 Nippon Steel Corp Method for controlling supply-demand of power, supply- demand controller of power, supply-demand control system of power and storage medium
JP2004112869A (en) * 2002-09-13 2004-04-08 Toshiba Corp Predictive system for electric power demand
JP2007004646A (en) * 2005-06-27 2007-01-11 Hitachi Ltd Power transaction support system and power transaction support program
JP2009294969A (en) * 2008-06-06 2009-12-17 Mitsubishi Electric Corp Demand forecast method and demand forecast device
JP2012034444A (en) * 2010-07-28 2012-02-16 Toshiba Corp Power supply-demand planning device and method thereof

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08110332A (en) * 1994-10-11 1996-04-30 Komatsu Ltd Method and apparatus for ultrasonic inspecting by articulated robot
JP2001318970A (en) * 2000-05-12 2001-11-16 Nippon Steel Corp Method for controlling supply-demand of power, supply- demand controller of power, supply-demand control system of power and storage medium
JP2004112869A (en) * 2002-09-13 2004-04-08 Toshiba Corp Predictive system for electric power demand
JP2007004646A (en) * 2005-06-27 2007-01-11 Hitachi Ltd Power transaction support system and power transaction support program
JP2009294969A (en) * 2008-06-06 2009-12-17 Mitsubishi Electric Corp Demand forecast method and demand forecast device
JP2012034444A (en) * 2010-07-28 2012-02-16 Toshiba Corp Power supply-demand planning device and method thereof

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
CSNG200500464001; 冨山芳幸: '気象リスクの定量的マネジメント' 経営の科学オペレーションズ・リサーチ第49巻第5号 第49巻第5号, 20040501, 282〜287, 社団法人日本オペレーションズ・リサーチ学会 *
CSNG200500464002; 冨山芳幸: '気温の確率予報を用いたリスク・マネジメント' 天気 50.3, 200303, 43-55, 日本気象学会 *
JPN6014045121; 冨山芳幸: '気温の確率予報を用いたリスク・マネジメント' 天気 50.3, 200303, 43-55, 日本気象学会 *
JPN6014045122; 冨山芳幸: '気象リスクの定量的マネジメント' 経営の科学オペレーションズ・リサーチ第49巻第5号 第49巻第5号, 20040501, 282〜287, 社団法人日本オペレーションズ・リサーチ学会 *

Similar Documents

Publication Publication Date Title
Zennaro et al. Initial conditions for accurate N-body simulations of massive neutrino cosmologies
Munoz et al. New bounding and decomposition approaches for MILP investment problems: Multi-area transmission and generation planning under policy constraints
Ying et al. Bi-objective reentrant hybrid flowshop scheduling: an iterated Pareto greedy algorithm
Lu et al. PRADA: Prioritizing android devices for apps by mining large-scale usage data
US20170322241A1 (en) Non-intrusive fine-grained power monitoring of datacenters
US20180268283A1 (en) Predictive Modeling from Distributed Datasets
Andradóttir A review of random search methods
US9208462B2 (en) System and method for generating a marketing-mix solution
Parnianifard et al. Comparative study of metamodeling and sampling design for expensive and semi-expensive simulation models under uncertainty
Chakraborty et al. A generally weighted moving average chart for time between events
US10075510B2 (en) Analytics-based update of digital content
Song et al. Multi-fidelity sampling for efficient simulation-based decision making in manufacturing management
JP2021064049A (en) Calculator system and mathematical model generation support method
US9342606B2 (en) Intelligent software installation
Gamboa et al. Decomposition methods for Wasserstein-based data-driven distributionally robust problems
Almomani et al. Selecting a good stochastic system for the large number of alternatives
CN115270013B (en) Method and device for evaluating emission reduction measures during activity and electronic equipment
US20230004560A1 (en) Systems and methods for monitoring user-defined metrics
US10298648B1 (en) Managing a distributed system processing a publisher&#39;s streaming data
CN111309712A (en) Optimized task scheduling method, device, equipment and medium based on data warehouse
JP2013182448A (en) Maximum supply amount prediction system, and maximum supply amount prediction method
JP2013156762A (en) Graph generation program, graph generation method, and production management apparatus
Quintana et al. Optimal allocation of testing resources for statistical simulations
CN113610174A (en) Power grid host load prediction method, equipment and medium based on Phik feature selection
JP2014164336A (en) Model analysis device, model analysis method and model analysis program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140305

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20141015

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20141028

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20141219

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150616

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150727

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20160209