国立防災科学技術センター研究報告第41号

1988年3月

550．34：681．3

地震前兆解析システムによる自動前兆監視
堀貞喜‡・松村正三＊・島田誠一榊・大久保正＊

井元政二郎＊舳・岡田義光榊＊

国立防災科学技術センター

Automaticby

theMonitorimg
Alla1yzing

SystemSystemforfor

Earthq11ake

Precursors

PmdictionofEarth叩akes

By

Sadaki

TadashiHori，Shozo

Matsumllm，Seiichi

Ohkllbo，Masajiro

Imoto

amd

Shima曲，

Yoshimitsu

Okada

Nα肋〃α1肋∫θακ乃C励θ7〃〃∫ω伽肋舳肋〃，伽〃Abstract

opedAn

automatic

atthefunctionofmonitoring

systemrate

theprimaryandconsists

threestagesofwave

modifieddatawave．Thearepossib1ybasesofdata）and1ow
to

the

data

sets

of

also

createdatwith

some

and

this

anomaly−leve1for

each

deviation

parameter

pattemorparameters

are

of

the

present1yparametersbegraphic

on

afor

a

variety

the

data

dai1y

va1uesanddisplay

from

for

of

are

eva1uated

the

va1ue

for

the

in

of

the

be

data

examined

on

their

analyzing
ordinarytenthestate．

items：seismicofstrainandradonanomaly−1eve1judgementcan

easi1yま第2研究部地震前兆解析研究室，舳第2研究部地殻力学研究室，

101

to

basic

magnitudes，spectra1parameters

termina1．Fina1ly，the

榊第2研究部，舳舳第2研究部地震活動研究室

Q■工ofsampling−ratethethe

decay，crusta1tilt，crustal

resu1ts

hypocentral

asconsidered

byfol1owing

sampling−

ana1yses．The
of

such

from

isdifferentof

the1owparametersparameters

highmaindata，etc．）aresets

parameters

ca1culatedparameter

monitored

ratio

waves，foca1mechanisms，coda

concentration．Thesedata（crusta1ti1tare

deve1−

automaticcompi1es

processing．These

activity，リーva1ue，b−va1ue，frequency

monitored

precursors

been

Earthquakes（APE）．The

stage．Secondly，various

earthquake

primaryandhasawaveform−related

minutely，hour1yprediction

Prevention（NRCDP）asprocessing．First1y，thereadydatathefluctuating

seismic

bases

atandprecursory

data

themusualness
Theofofsampling−ratesamp1ing−ratesolutions

DisasterPrecursors

fmdamenta1data
high

re1atedpreparedforfor

processingcodaforμα6此α1earthquake

CenterSystemparameters，foca1mechanism
asystem

Ana1yzingdata（seismicprocessed

monitoring

Nationa1Researchsystemoutputstheregiona1

国立防災科学技術センター研究報告
anoma1y−level

andobtained

combining

threebyprocessings
eachsuggeststhethe

probabi1ity

individualofof

processing

scientificdeveIopment

in

the

ancanbeeasilyfie1dofearthquake

1988年3月

earthquake

anoma1y−1evels

areヵ〃一α〃o〃〃此α妙operated

procedure

第41号

bya

modifiedoccurrence

ca1cu1ated
rea1time

by

which

eachcomputer．The

corresponding

to

any

has

been

item，Thesemethodortechnica1orprediction．

1．はじめに
国立防災科学技術センター（以下「防災センター」と略す）は，昭和60年度末に「地震前兆解析システム（TheAna1yzing

SystemforPrecursorsofEarthquakes，以下APEと略す）」を導入した〔松村・他（1987）〕．APE導入の目的は，単に増大したデータを迅速に処理するためだけではなく，短期的地震予知の実用化に向けて，従来それぞれの専門家が行ってきた地震を予知するための様々な解析とその結果に対する判断を，計算機が自動的に行うシステムを実現することにある．従って，本稿で述べる地震前兆自動監視のためのサブシステム（以下「自動監視システム」と略す）は，APEの中核を成す部分になっている．この部分はAPEを構成する3台の計算機（1号機，2号機，3号機）のうち2号機によって行われる．具体的な地震前兆の検出手法や，地震発生の危険度の数値的表現，例えば地震発生確率等を算出するためのアルゴリズムについてはまだ開発途上にあるが，処理の流れ，起動タイミングやデータの格納形式等，システムとしての枠組みは，昭和61年度をもって完成した．本稿ではこの枠組みの内容について述べる．

2．自動監視システムにおける処理の種類

自動監視システムにおける処理には図1に示すように，大きく分けて3つの段階がある．
まず，入カデータの種類によって，高速採取データ（80Hzサンプリングデータ＝地震波形データ）に対する部分と，低速採取データ（1Hzサンプリングデータ＝地殻変動などのデータ）に対する部分の2つに分かれ，それぞれに対して1次処理と2次処理がある．1次処理とは，生データに対して定常処理で行うような加工をし，最も基礎的なデータを作成することである．

2次処理とは数種類の異常監視項目を設定し，1次処理による基礎的データから，それぞれのパラメータを算出し，その変動状況を監視することである．さらに最終処理として，高速採取データ，低速採取データ両者の監視結果を統合し，検出されたパラメータの異常性が地震前兆に結び付くかどうかの判定を行う総合異常判定処理がある．高速採取データに対する1次処理としては，震源計算，メカニズム決定，および波形情報の検測／松村・岡田・堀（1987）〕が行われる．また，2次処理としては地震前兆との関係が想定されている各種異常監視項目

一102一

地震前兆解析システムによる自動前兆監視一堀・松村・島田・大久保・井元・岡田

HlGH

SAMPLlNG−

RATE

LOW

SAMPLING−

RATE

DATA

DATA

PRIMARY
PROCESSING

PRIMARY
PROCESSINGSECONDARYPROCESSINGPROCESSING

SEC◎NDARY

F1肚PR㏄ESSlNG（ANOMALY
図1
Fig．1

JUDGEM日NT）

地震前兆白動監視システムにおけるデータ及び処理の流れ
Flow

chart

monitoringsorsofinthe

automaticsystemforprecur・earthquakes．

（後述）についてそのパラメータの算出及びそれらのパラメータの異常性の判定が行われる．

一方低速採取データに対しては，1次処理として，収録されたデータをチャネルごとに並べ替える編集処理が実行される．次に，2次処理として1次処理の結果から降雨・地球潮汐等，その挙動の分析が可能な「非前兆的変動」の影響を除去し，残存量の変動パターン等からその異常性が判定される．両者の2次処理結果は総合異常判定において次のように統合される．

関東・東海地域を図2a〜cに示すような約60の領域に細分し，その個々の領域ごとそれぞれの地殻活動状態の異常性が総合され，最終結果としては異常性の度合が出力される．地震予知を実現するためには，発生が予想される地震の規模，発生時刻，及び発生場所を推定しなければならないが，現在のところそれら全てを予知するための決定的手法はまだ確立されていない．従って，自動監視システムではその手法を模索するためのデータを蓄積する意味で，このような出力を」応の最終結果としている．最終処理のため監視対象とした小領域群は，主として関東・東海地域における微小地震の震源分布をみることによって分割されたものであり，個々の領域は多面体で構成されている．領域を多角柱ではなく多面体で定義することによって，各領域は3次元的に分割されており，たとえば斜めに沈み込んだプレートの形状などもそのままの形で領域分割に反映するように工夫されている．これらの領域の大部分は，震源が密集する領域を取り囲むように設定されているが，地殻上部は，監視対象領域が隙問なく割り当てられており，地震活動度が非常に低い領域も存在する．

■103一

国立防災科学技術センター研究報告
NRCDP

REGlONS｝｝｝SECTlON

16

ヨアN

35

酬㌧

38

17

NO．1

一，14

埠

19・1亭6

1988年3月

‡ウ13

ρ一．

第41号

524，．g2弓212β7−j．244o

1035刊

一3942287629

ザ

図2（a）

自動監視システムの監視対象

Fig．2aDistributions

ofmonitoredthe11・1

領域断面1（浅部）

4134■3413ア…

1眈NRCDP03233EREGlONS

1！9■E

12140・E＾｝SECT■ON

14一一E

142

E

Sectionbytheregions

system．

No．1（Sha1low）．

NO．2，榊

46

51

48

g．
50

43〈！5・㎝

44

図2（b）

自動監視システムの監視対象
領域断面2（中深部）O

Fig．2bDistributionsofmonitored

theSection

NRCDP

REGlONS■｝｝SECTlON

NO．3

by

the

regions

system．

No．2（Intermediate）

｝

45
52

54

彰庫
l54伽L一＿．．

53

Q、l

図2（c）

57

自動監視システムの監視対象
領域断面3（深部）0

Fig．2cDistributionsofmonitored

theSection

104

by

the

regions

system．

No．3（Deep）．

地震前兆解析システムによる白動前兆監視一堀・松村・島田・大久保・井元・岡田

3．自動監視システムにおける各処理の内容3．l1次処理
3．1．1高速採取データに対する処理
図3に高速採取データに対する処理の流れを示す．この処理はデータ収録の処理が終了すると直ちに起動される．まず磁気ディスク中に収録された高速採取データは，検測処理用に並べ替えられた後，磁気ディスク中に編集済み波形ファイルとして保存される．このファイルは9．6MBの容量を持っており，1〜2週間程度の波形データが保存可能であり，同データを使ってグラフィックディスプレイ（GD）上に，後に行う検測結果と共に波形を表示することができる（図4）．得られた編集済み波形に対して，PおよびS波到着時刻の自動検測が行われ／堀・松村（1987）〕，続いて震源計算，メカニズム決定／岡田（1987）〕，そして波形情報の検測が実行される．これらの計算結果は磁気ディスク中にそれぞれ別ファイルとして格納される．

震源ファイル，メカニズムファイノレ，そして検測パラメータファイノレは，2次処理におけるPROCESS！NロH

ロF

Hl〔；HS＾MPLlGHS＾MPL

NG−R＾↑E

PRl

M＾RYMECH＾N

R巨ADl

NローR＾TED＾丁＾D＾↑＾

ROCESS1NG
DETERM

HYP0［ENTERF口C＾L

Il

NAT10NSM DET巨RM

llN＾T10NWAVE−REL＾TED P＾R＾METERSNCW＾VE．

MEC卜1．

HYP0．

SE〔：口ND＾RY PROCESSl

NG＾NOM＾LY−LEVEL JuD口EMENT BY 1TEM
一．

1■一■・1一・・

SE1SM■・l

●■C

1＾CTl

VI

■

一一TY

リーV＾LUEb−V＾LuE
M＾口N1TuDE−FRE口uENCY

FOC＾L MECH＾NI SM

SPE［TRuM
CODA DEC＾YP＾R＾METERS TO BE MON1TORED

図3

高速採取データの流れ．HYPO．：震源要素のデータセット．MECH．：発震機構解のデー
タセット．WAVE．：波形検測パラメータのデータセット

Fig．3

RESuLTS 0F JUDCEMENT

F1owofF】NAL

chartofprocessing

forthehighsamp1ing−rate data．HYPO．：Dataset of
hypocentral parameters．MECH．：Dataset

PROCESSING

foca1mechanism

Dataset一105一ofsolutions．WAVE．：

wave−relatedparameters．

国立防災科学技術センター研究報告

HlGH＝1〉NMT

SAMPLlNG

RATE

第41号

1988年3月

I DN02960
GAlN
口1．OO

DATA

0ATE 87！09！21
STATlON
NMTuD

P
SC

NMT

NS

S

NMT

EW

OO＝10：28

図420468101214161820

地震波形のモニター画面．PA一：P波の到着時刻，読み取り精度のランク，初動の方

向．P：P波の計算走時．SC：S波の到着時刻，読み取り精度のランク．S：S波の計算走時
Fig．4

Examp1eofoftheaccuracy

andCa1culatedP

Calculated

S

monitorpo1arityofscreentheofaseismic

initia1motion

time，SC：Arriva1time

andofwaveP

rank

data．PA一：Arriva1time，rank

wave

of

detectedbycomputer．P：

accuracyofSwave．S：

time．

便用目的を考慮し，それぞれ2通りの形式で格納されている．まず，主としてGD上に表示する目的で発生時間順に並べたデータファイルがある．容量は，震源ファイノレが地震1500個分，メカニズムファイノレは225個分，そして検測パラメータファイルは15000トレース分であり，いずれもほぼ1〜2ヵ月分の情報が保存される．もう1つのデータファイルは，領域毎に分けたもので，容量は各領域とも震源ファイルが地震300個分，メカニズムファイルが45個分，そして検測パラメータファイルが3000トレース分となっている．震源データは図5に示す形
式で常時GD上に表示されているが，より信頼度の高い震源マップを表示するため，10日毎に3号機（FACOM

M360）で行われる定常処理の震源計算結果を2号機（FACOM

S3500）側に転送し，2号機による自動処理結果とあわせて同一画面上に色を変えて表示している．通常は両者の重複表示を避けて，定常処理結果についてはその最新データまでを水色で表示し，それ以降から現在までの期問についての自動処理結果を緑色で表示している．さらに自動処
理結果のうちの最新の震源位置は赤色で点減表示させ，より識別し易くしている．また，自動震源計算については，即時決定のルートが別に用意されており，この結果は同画面上に黄一106一

地震前兆解析システムによる自動前兆監視一堀・松村・島田・大久保・井元・岡田

1987

05

22〜

19．87

06

16

N昌

O＾TE

g47

古体

37．0N．ぜ、｝

l1

87．06．一6

TIME

1全li．43．i1

o ；（ROUTIHE）

於1・・肚舳・明物（日UT㎝PTIC）。1出。｝日附．昨閉・目・m・二・｛；ξ□外二仙1＾ 由．

、ぺ
ぺ・1〆；

申．ザ

触、1・＝、㍑，．．．．

一、．、≡；淋・由ム・．、．

グ邑・い！

二削二・：・・一・・

1饒1：：：；1；簑、㌔義・・■

LPTESTEuENT

DEPTH：15．ヨRh
直省ぺ｛

MP6トlITUDEぺ。い・箏出＝ピ垣。。．㎝oム≡
芦震源分布のモニター画面例
Examp1emonitor

screen

DEPTH

φ0

図5

the1OOoof○。．。
＾■

Fig．5

2．O工旧くヨ．0

8：1：1；：：1：（1くm）〇三Zく1010＜Zく≡≡O

口20工ZくヨO＾ ；OくZく60
古 60＜Zく90廿 90くZof

hypocentra1distribution．

色で表示されるようになっている．ただし，これに続いて実行される自動震源計算が終了した時点で，結果は自動的に置換されるようになっている．すなわち，1つの地震の震源位置は黄色から赤点滅，そして緑を経て水色に変わっていくわけである．別ルートを用いて即時的な自動震源計算を行っているのは次のような理由による．本来の自動監視システムにおける白動震源計算では，結果の信頼性を重視するため，波形の収録を完全に済ませてから計算を開始するため，特に大きい地震に関しては，震源決定結果が見られるまでにかなり時間がかかってしまう．これに対して信頼性はある程度犠牲にしても，より即時的な状況の把握を行えることが業務的な要請として提出されたからである．この即時震源計算では，扱う波形データをトリガオンとなった時刻を含む3秒問に限定し，P波の検測値のみを用いて震源計算を行っているため，地震発生から1分以内で計算結果を表示することができる．

3．1．2低速採取データに対する処理
図6に低速採取データに対する処理の流れを示す．低速採取データに対しては，1次処理として分値，時問値，そして日値の自動編集処理が行われる．時間値と日値作成の起動タイ

ミングはそれぞれ毎正時と0時であるが，分値作成処理だけは，10分毎に過去10分分のデータをまとめて処理することになっている．ここでの処理は，データをチャネノレ別に並べ替えることのみで，3号機においてオペレータが定常処理として行う補間やオフセットの補正は

一107一

国立防災科学技術センター研究報告

PROCE1SSlNG

LOW

OF

LOW

SAMPL

SAMPLlNG−R＾TEING−RATE

第4！号 1988年3月

DATA

DAT＾

PRlM＾RY PROCESS lNG
1三Dl↑1NG

OF MlNuTELY DATACREAT10N 0F HOURLY DA↑A
CRE＾丁一0N 0F DAILY DATACREA↑10N

H．

M．

D．

SECONDARY PROCESS1NG
ANOM＾LY−LEVEL JUDGEMENT BY 1↑EM
CRuSTAL T−LTCRuS↑＾L STRAIN
RADON CONCEN↑RAT lON一I．

一・一・一

一

一・．

・…

…

．

．

．一．

・

●一・・一・．1

図6P＾R＾MEl↑ERS

T0 8E MON−TORElD

ツト
Fig．6

RESuL↑S 0F JuDGEMENT

低速採取データの流れ．M：分値のデータセット。H：時問値のデータセット．D：日値のデータセ
Flowchartofprocessing

for

the

lowsampling・ratedata．M：
Datasetofminutelyva1ues．H：

Datasetofhour1yva1ues．D：

FlNAL

PROCESSING

Datasetof

daily

values．

行わない．

3．22次処理3．2．1高速採取データに対する処理
2次処理は震源やメカニズムなどの基礎データから，地震前兆を検出するための各種異常監視項目についてそのパラメータを計算する部分と，それらのパラメータの変動量，過去の履歴等から対象領域における異常性を判断する部分から成り立っている．異常を判定する部分では，各パラメータの平均値と標準偏差を用いることが多いが，これらを毎回計算するのは余り効率的でない．そこで，主にこの2つを格納しておくために，異常判定用パラメータを確保するためのファイルが磁気ディスク中に設けられている．現在高速採取データに関して設定している異常監視項目は，図3に示すように地震活動度，b値，リ値，規模別頻度，発一108一

地震前兆解析システムによる自動前兆監視一堀・松村・島田・大久保・井元・岡田

震機構解，スペクトル，そしてコーダ減衰の7つで，これらの起動のタイミングには2種類ある．地震活動度とリ値は定時起動で，1時問おきとなっている．その他の5項目はいずれも1次処理終了時，つまり地震検出時に起動がかかる．また，前者については全領域について計算が行われるが，後者については地震の発生した領域についてのみ行われる．ただし，後者のうちコーダ減衰については観測点単位で計算が行われることになっている．それぞれの項目の出カとしては最大5個のパラメータを想定しており，2次処理の前半部分の処理結果，すなわち各項目について算出されたバラメータと，後半部分の処理結果，すなわち異常性の判定結果は磁気ディスク中にそれぞれ保存される．その容量は1領域あたり，定時起動項目
一一一一一一一一一一一一一一一一一く

前兆指標モニターメニュー）一一……一一…一一一一一一一一

対象：高速採取データ

起動1地震時区分：領域別領域番号＝・）26

項目名・・）BV＾LUE

パラメタ名（赤〕・⇒BV1（緑）・・〉BV2（青）＝＝〉

データスケール
最大値（赤）・・）

2．00（緑）＝・）

2．00（青）・＝〉

O．0

最小値（赤）三⇒O．0（緑）＝・〉

O．O（青）・・）

O．O

表示期間（日＝時＝分）・・〉60

終了時刻（年＝月＝日＝時1分）三⇒

図7（a）地震前兆パラメータ（b値）のモニター画面例

メニュー画面
Fig．7aExampleofthemonitorparameter（ろ・value）．

screenMemofaprecursomy

screen．

NRCOP

APE

DA

SYsTEM〈＜〈B−VALUEREGlON NO．26
OFF E IZU PEN−8〉〉〉

TlME

87TE

87．10．15

1？．19．47■■■ EVENT TRiGGER ，■■07！18α〕＝α〕一8711011600：O0

8

BV1

FBV2

図7（b）地震前兆パラメータ（b値）のモニター画面例

モニター画面
Fig．7bExampleof

themonitor screen of aMonitor screen．

parameter（6−value）．

一109一

precursory

国立防災科学技術センター研究報告第41号

1988年3月については60日分，地震時起動項目については300地震分，そしてコーダ減衰の場合は，1観測点あたり300地震分となっている．これらの計算結果は，図7aに示すメニュー面画を経由して図7bのように，横軸を時問とした任意の時問幅でGD上に表示される．時間軸は，図7aの画面上で表示期問と表示終了時刻を入力することによって自動的に設定されるが，表示終了時刻に空白を入力した場合，図7bの画面はリフレツシュ機能を持つことになり，常に最新のデータが表示されるようになっている．

3．2．2低速採取データに対する処理

1次処理で作成された分値，時問値，日値それぞれの編集済み低速採取データは，地球潮汐，降雨，気圧等の影響による変動だけでなく，地震や停電時のショックによるオフセットなど，地震の前兆を検出するためには妨げとなるものを数多く含んでいる．2次処理ではこ
うした非前兆的変動を除去し，後に残ったものを真の地殻変動に関わるデータとして異常監視項目の一つに加え，磁気ディスク中に格納する．そしてそのパラメータの変動パターン等をモニターすることによって異常な変動をとらえ，その情報も磁気ディスク中に格納する．
非前兆的変動を除去する場合，例えば，気圧補正係数や降雨に対するタンクモデノレの定数等，固定したパラメータを用いることが多いので，各チャネル毎に124個の変動除去用パラメータを磁気ディスク中に確保できるようにしている．低速採取データの場合，それぞれの観測点に付随するデータとなるため，個々の処理における計算結果は監視対象となる領域毎にまとめられるということはない．また低速採取データのチャネノレの中には，雨量，気圧，温度等それ自体は異常地殻変動の監視対象にならないものも含まれているので，上記の処理の対象としては，データ採取時の960チャネノレの内から，最大150個の「観測点項目番号」で識別されるチャネルデータが選択される．個々のチャネノレデータは，最大8個（異常判定結果の場合は最大12個）のパラメータを持つことができる．低速採取データの場合，現在のところ，観測点数は全てを合わせて30点に満たない．また，観測している項目も地殼傾斜，地殻歪，及び地下水中のラドン濃度の3つで，将来観測点数が増えたり，全ての観測点でこれら3項目の観測を始めたとしても，設定された容量で対応できるようになっている．磁気ディスク中の各ファイノレの容量は，分値が8日分，時問値が60日分，日値が960日分となっている．処理結

果は，高速採取データの場合と同様，メニュー画面で必要なパラメータを入力した後，GD上に表示することができる．

3．3最終処理
図8に最終処理の流れを図示する．ここでは2次処理で作成されたデータをもとにして，各領域に対して総合的な異常性の度合を表すパラメータ（以下「領域別異常指数」と称す）を最大5個まで算出する．最終処理の起動のタイミングは1時問問隔である．この処理では，まず各観測点毎に計算された異常判定結果（低速採取データと高速採取データの内のコーダ減衰パラメータ）を領域毎の異常指数パラメータに換算する．そしてすでに各領域毎に言十算さ一110一

地震前兆解析システムによる自動前兆監視一堀・松村・島田・大久保・井元・岡田

F

INAL

PROCESS＝NG

RESULTS 0F ANOMALY−LEVEL JuDGEMENT BY 1TEM
Hl GH SAMPしlNG−RATEDA↑A LOW SAMPL lNG−RA↑E DATACODA DECAY

TRANSしA↑10N T0 REG10NAL DATA

（R1三⊂；10NAL〕

JUDGEMENTANOMALY−LElVELRESUL↑S 0F JuDGEMEN↑

ALARMINGh■■■■D＾＾l，＾SPLAY＾■図8

最終処理におけるデータ及び処理の流れ

Fig．8

Flowchartof

thefina1processing．

れている高速採取データに関する異常判定結果と合わせて評価することによって，領域別の総合的な異常指数を計算する．5個のパラメータの割当を，例えば，異常判定におけるデータの時間幅によって分けてやれば，それぞれの潜伏期間の異なる地震前兆現象に対応してその領域で地震の発生が予想される時刻を推定するための資料になることも考えられる．この計算結果は磁気ディスク中に1年分保存される．そして，各項目別の異常判定の場合と同様，リフレッシュ機能を持つモニター画面上に以下のような2通りの形式で表示される．一つは，

図9aのように，指定した1領域について，横軸を時問に取って任意の期問だけ表示する方法である．もう一つは，図9bのように任意の時間断面の空問的な分布として，異常指数の大きさによって各領域を色分けして表示する方法である．前者の場合，画面は上下2つに分割され，上半分にはその領域，ないしは指定した観測点における項目別の異常判定結果が最大3パラメータまで表示される．下半分には，指定した領域における総合異常指数が表示される．

4．議論

以上で述べたようにAPEにおける地震前兆自動監視については，処理の起動タイミング，

一111一

第41号

国立防災科学技術センター研究報告

nCOP
＾PE

〈〈〈REGlONAL

SYS．1■EM

RE610N

NO．26

ANOMALY0FF

FOR

EARTHQUAKE＞
←〉

50＾TE

87．06．16l

16．28．42

TE

87／0411800：00−87／06／17001◎0Uu」〈事⊃一くq〉o⊃」

り血o2o：

OCCuRENCE〉〉〉

E 1ZU PEN．

8s

1988年3月

8V2（9〉ou』

8V1

ω暑8
亨

88

AER28
亨

竃dNooCP2η阯o−o−o8

OP．1

80P3
11一■
APR1987

M＾Y1987

図9（a）

地域別異常判定結果の表示例

Fig．9aExamp1eJ∪N1987

時間変イヒ

of

themonitor

screenfor

resu1ts

R［GION＾L＾NO＾LYN＾1〕ION｝｝｝

ofregional

anoma1y−1evel

judgement．

Tempora1variation．
ト旧Cl）PnPE

くく〈SYSTEホホ｝

SECTNO．

1〉〉）
87／03■20）・乃・一31

16．
・斗5・1

3

㍗︑・35

17O＾T［8r．U3．26

T1N

lE

l．52．20

21＝00

？

，戸i、！

・1・§・、■，i9q

9一＝、1381

2幸1

129
1

千0・一、、！O

1享・3・＝帖ヂ、；．、、

f一」

l1

躰．

冷・ヅ■

8

戸G

P1＊lO

l27

一津6・・

3，0

む3士り．

二．2．．

34一／

33
32■■一

図9（b）地域別異常判定結果の表示例
空間分布

Fig．9hExamp1e
Spatialof

themonitor

screenfordistribution．

112

resultsofregional

anoma1y−1eve1judgement．

地震前兆解析システムによる自動前兆監視一堀・松村・島田・大久保・井元・岡田

資源の割当，モニターの形式等，システムの枠組みに相当するプログラムが完成されている．

一方，地殻活動の異常検出のために設定した各種の項目については，個々のプログラムが独立して開発され，昭和62年10月現在では，予定した10項目のうち地震活動度，ひ値，b値，規模別頻度，地殻傾斜，そしてラドン濃度の6項目が組み込まれている．これらについてはまだテスト的な運用が始まったばかりであるが，今後はこれまでの研究によって得られた知見を生かして，実際に異常を判定するためのプログラムの開発を推し進めていかなければならない．今回完成したシステムの枠組みでは，こうしたプログラムの開発や改良が簡単にできるような環境が設けられている．例えば，各項目の処理は，バッチジョブとして次々に起動されていくのであるが，おのおのに対するイニシエータは，各起動タイミング毎に1つになっているため，一般のリアルタイムプログラミングに際してのように，各ジョブ問の処理の交錯に注意を払うことはそれほど大きな問題にはならない．また，汎用言十算機である3号機の方に，シミュレーションのできる環境が用意されているため，プログラム作成・デバッグを3号機で行ってから完成したプログラムをリアルタイム機である2号機に移行することができる．新しいプログラムについて自動監視システムでの運用を始める場合には，起動タイミング別に用意された各ジョブの起動を行うプロシジャ中に，新しいプログラムを起動する部分を挿入するだけである．また，すでに運用されているプログラムを修正する場合も，スクリーンエディタの機能を用いてディスプレイ上で修正した後，翻訳・結合を行うのみで，次回の起動時からは修正後のプログラムが走り始める．

現在の地震予知体制は，小数の専門家の知識や経験に基づく判断に頼っているところが少なくない．しかし，観測網の整備などが進むにつれ処理すべきデータ，目を通すべき資料も膨大なものになりつつあり，現体制では対応できない局面が生ずることが予想される．将来的には，これまでの地震学的知見を反映させた計算機システムに，判断を含めた処理の大部分を移行していく方向になっていくであろう．その意味でも自動監視システムはきわめて先駆的なものであり，その成否は新しい地震予知体制を探求する上で重要な意味を持ってくるといえる．

謝辞

地震前兆自動監視システムの開発は小西秀之氏をはじめ富士通株式会社のスタッフの精力的な作業に負うところが少なくない．ここに厚く謝意を表する．また，最終処理における領域分けは，国立防災科学技術センター第2研究部の野口伸一主任研究官によっている．

参考文献
1）堀貞喜・松村正三（1988）：地震前兆解析システムにおける自動震源決定，国立防災科学技術センター研

一113一

国立防災科学技術センター研究報告第41号

1988年3月

究報告，41，89−100．

2）松村正三・岡田義光・堀貞喜（1988）：地震前兆解析システムにおける地震データ（高速採取データ）の処

理，国立防災科学技術センター研究報告，41，45−64．
3）松村正三・岡田義光・井元政二郎・島田誠一・堀 貞喜・大久保正・大竹政和・浜田和郎（1988）：地震前兆解析システムの機能と構成，国立防災科学技術センター研究報告，41，35−44．
4）岡田義光（1988）：震源計算・発震機構解計算プログラムの改良，国立防災科学技術センター研究報告，

41，153−162．

（1987年11月16日

114一

原稿受理）