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

1988年3月

551．465／．508．5／．55

発達した低気圧による海上強風の乱流特性と長周期スペクトル構造
内藤玄一＊
国立防災科学技術センター平塚支所

Tl1rbu1emt

Propertiesa．ndLong−period

Spectm1Behaviours

Winds

Deve1oped

CausedbyofStrongCyc1㎝es

By

Gen

ichi

Naito

舳淋〃肋8〃〃ピ／1，Nα〃o〃／肋∫θ肌ソlC〃θ7〃1）ゴ∫α∫伽・〃20〃ゴo〃，9−2，」V批gα11α閉α，H加αお〃尾α，Kα〃αgo〃α一加〃254，1αψα〃Abstractusing

Long1astingstrongaofcoupleevaluated

with

Gust

winds

sonic

discussion

factoris

wereobservedofbehaviours

theatmosphericstationary

becomes1arge

bydeve1oped

cyc1（〕nepasses．Turbulent

different，vith
Power

by

the

the

coefficient

under

the

wind

of

largeroftheobservation

tion

tal

kind

spectra

long

are

of

condition
vectors

than

does

longitudinal

more

large
not

the

the

intensities

surface

of

the

and

tower

by

statistical1y

disturbances．

turbu1ence

about3．2whenthoseoff−shore

wereofstrongboundary

from7mcausedcoefficientunderturbu1ence，and

are

drag

1ayer

the1ayer

disturbance

disturbance．
vectorvectors

verticalexponentia1Iy

is

wind

than400min，andhorizonta1wind
of

to

of

properties

theabruptly

disturbancesheightsofhigh，butstrongtwo

turbulent

fromユ．46to1．36under

to23m
the

at

anemometers．Those

dency

of

mean

wind

speed

not

cannotexpressthespatia1correlation

in

very

wide

those

featuresdecreaseswithfrequencies

were

decay

with

increasing

c1ear．Therefore，the
of

wind

turbulences．

ホ沿岸防災第2研究室

341

f1uctuations

mean

wind

coherenceoffrequency，and

its

Davenport
1argerobtained

discussed．Correla−increasing

separation．Additiona1ly，eddy

weregeometric
than

speed

horizon−depen−

simi1arity

surface

1ayer

Ocean

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

1988年3月

1．はじめに強風による災害は家屋や樹木の倒壊，自動車の横転，海洋上における船舶の座礁などの直接的なものから，高波，高潮，塩害などの問接的なものまで広範囲にわたっている．これらの災害を防止するためには先ず強風の発生原因と強風の性質を充分解明することが不可欠である．そして強風の時問的空問的挙動に関する知識の積み重ねは，災害防止に最も必要な都市計画，耐風設計に資するものである．沿岸地域に生活圏が集中する我が国において，これから長大橋，高層建築物など増々構造物が大型化する方向にある．この条件下では地震よりも現象の時問スケールが大きい風の問題が，相対的に比重を増してくる．

我々が生活している空間は，ほとんど地表面から数100mまでの大気境界層の内である．従って大気境界層における強風の挙動が災害の要因となる．境界面すなわち地面（または海面）の特徴がその上の気層の性質を支配するが，乱流よりも大きな上空の気象擾乱も地表面近くにまで達することが少なくないと考えられる．

海上風は地形に影響されないため，大気境界層の性質を明らかにするのに最も都合がよい．

ただし海面の空気力学的粗度が陸上の草地などと比べてかなり小さいため，乱流の特性が両方の大気層でかなり異る．陸上における強風は，気象擾乱によって生じる本来の大きさに複雑な地形効果が加わるため局地性が非常に強い．他方海洋上では非常に長い吹送距離が期待できるため，気象擾乱のみによる変質されない強風が観測できる．

境界面に近い接地層の乱流構造については多くの研究者によってかなり詳しく調べられているが，ほとんど陸地面上での気流に関してである．海上の気流の構造についての実験的研究は少なく，Kondoθ≠α五（1972）は波浪との相関に関して貴重な成果を得た．境界層の気流のスペクトノレ構造は，Davenport（1961），Hino（1971），Kaima1θ左αム（1975）らがパワースペ

クトルなどの実験式を提案することによってモデル化をはかった．しかしこれらのスペクト

ル形はおよそ30分問の測定時問における風速変動に適用されるものであり，大気の長周期の挙動を無視している．長時問にわたる風速のパワースペクトルはVander

Hoven（1957）の他にほとんど報告されていない．大気境界層の乱流特性は，乱流よりも長い周期をもつ大気擾
乱と切り離して論じることが出来ない．したがってかなり広い周期帯における連続スペクト
ノレを実験的に得ることは，境界層の構造を理解するためにも非常に重要なことと考えられる．

大気乱流の空問相関に関する研究は，Davenport（1961）の幾何学的相似則の提案以降，急速に進展した．Shiotaniθ左α／．（1976），Pie1keθオα五（1970）らは，乱渦の周波数別相関を表現するコヒーレンスの性質を実験から追求した．Naito（1982）は海面上の接地層で3次元風ベク
トルのコヒーレンスを，測定データより詳細に調べてDavenportの相似則を発展させた．

今回の報告は主として，時間スケールの大きな乱流の特性を明らかにすること，及び強い

気象擾乱が接地層の乱流に干渉する程度を評価することを目的としている．そのため長時問

一342一

発達した低気圧による海上強風の乱流特性と長周期スペクトル構造一内藤にわたって強い海上風が吹く条件下で，海上の観測塔の2高度に設置された超音波風速計によって連続測定が行われた．強風を対象としているため，大気状態の安定度は中立であるとし，気流の解析には安定度は特に考慮していない．

2．測定

強風の測定は相模湾の最奥部にある当所の波浪等観測塔で行った．相模湾は南方に向って開いており，南寄りの風に対しては無限吹送距離が期待できる．また後背地は平野であり，丘などによって気流が大きく変質することがない．

風の測定には2台の3次元超音波風速温度計（海上電機，モデノレDAT−300）を用いた．各々の受感部は塔自体による気流の乱れを避ける配慮をし，更に南寄りの風を想定して取り付けた（写真1）．下方の受感部の高度は平均海面より6．85m，上方の受感部の高度は23．Omである．両方の風速計は海面よりの同一鉛直線上にある．また海面の波浪状況は容量型波高計の出力から求めた．良質な測定値を取得するために，塔に設置されている他の測器によって風写真1P1loto1

2台の超音波風速温度計による観測
Obser▽ation

by

using

a

cou−

p1e of sonic anemometer・
thermometers．

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

1988年3月

向風速，気温，水温等を常時観測して測定条件の総合的な判断を行っている．

風速変動と波高の出力は観測塔内でディジタル化され，オンライン・データ収集システム
で陸上の計算機に伝送される、3次元風ベクトルは0．06sで読みとり，解析の目的に従って，

10分，30分及び200分を一つの単位として統計処理する．スペクトル解析の場合，200分以上の連続観測データに対しては全体のスペクトル挙動が明らかになるように処理する．解析方法としてFFT法を用いた．

3．結果と考察3．13次元風ベクトルの乱流特性
海上風の3次元成分の乱流特性については，Naito（1978）［は比較的定常な場合の観測値か

ら乱流強度，海面の抵抗係数などを示した．長時間にわたって強風が吹くとき，必ずしも充分定常な状態が続かない．強風が観測される場合，その原因となる大気擾乱，例えば台風，

メイ・ストームをおこす発達した温帯低気圧などの中心が時刻とともに移動する．またそれらの大きな低気圧に伴う中小規模の擾乱が速く変動する．従って地表面で測定される風向と風速はかなりよく変化する．しかし別な角度から強風時の風向風速の変化をみれば，その発生原因となる大気擾乱を捕えていると言えよう．

今回解析の対象とした事例は日本海を通過した中程度の強さの台風と，前線を伴っている

発達した温帯低気圧による強風である．いずれの場合も12時間以上連続して10m／s以上の平均風速が観測された．

地表面近くの乱流境界層では，大気の安定度が中立であるとき平均風速σの鉛直分布は高度2の対数分布で与えられる．すなわち

一

σ二

〃＊

K

肋

2

（1）

2o

ここでK＝0．40はK直rman定数，〃。は摩擦速度，2。は空気力学的粗度である．〃。は大気乱流の運動量の鉛直輸送スケールを表わすパラメーターであり，接地層では一定と考えられる．（1）は風浪が発達するとともに2。が大きくなる海上においても成立するが，海面の2。は極めて小さい．今回の一連の観測で，高度6．85mの風速変動より直接求めた〃。は非常に合理的な値であった．他方高度23．Omの風速変動より求めた〃。は，高度6．85mで求めた値と異るかあるい
は算出不可能（運動量フラックスが海面へ向し）ていない）であった．この原因は風の鉛直シアーが高度23．0mでは小さいため，測定が困難になり相対的に誤差が大きくなるという理由だ
けでは説明できない．（1）が成立する定フラックス層すなわち接地層は海洋上では薄く，高度20mを越えると異った性質を示すようになると思われる．

一344一

発達した低気圧による海上強風の乱流特性と長周期スペクトル構造一内藤

平均風速の鉛直分布として
σz

T＝（丁）α（2）の形がかなり高い層まで成立する．ここでσユは基準高度21での平均風速，αは地表面の粗度によって決まる定数である．図1は（2）が成立すると仮定して高度6．85mと23．0mでの10分問
平均風速よりαを求め，21＝6．85mでの平均風速σ1に対して示したものである．

風速分布として近似的に（1）と（2）が同時に成立する境界層では，両式より次の関係が導
かれる．

1

〃＊

1

α一一一 ＝一πK σ1 K（3）

CDは海面の抵抗係数である．図1よりαはσ1の増大とともに，わずかに小さくなる傾向をもつ．後述するがCoは風速が大きくなると，少し増大する．従って上式よりαはσ1とともに少

O．3

0．2

㌃雑燃＿

O．1

図1

平均風速に対して示した風速鉛直分布

のベキ指数
5

10

15

Fig．1 Power of the vertica1wind
against mean wind speed．

20

σ1（π／・）

profi1e

しずつ増加する性質を示すはずであり，図1での風速依存性とは逆になっている．しかし測定値のばらつき程度を考慮すると大きな矛盾を与えるものでなく，αはσ1に対して一定と考えて良いだろう．図よりσ1＞10m／sで，α＝0．13〜0．18の値をとり，平均的には，σI＝14m／s近傍で，α＝0，145（≒1／7）となる．この値は今までに報告されているαの値とほぼ等しい．

図2に2高度における10分問の突風率Gを平均風速に対して示す．突風率は平均風速に対する最大瞬問風速の比でもって定義されるため，風速計の応答特性に依存する．今回のデー
タ処理における評価時間（読み取り時間）はO．24sである．図より上方の高度のGよりも海面に近い高度のGの方が少し大きい性質が見出されるが，二つの高度でのGの差はあまり大きくない事が分る．同図の破線に囲まれた一群の測定値及び右上方の矢印を付した測定値を除くと，あとの測定値は海洋上の定常な境界層乱流によって与えられるものである．この条件下

一345一

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

1988年3月

でGの平均風速に対する依存性は小さく，風速の増加とともに少し小さくなる．突風率の中一こ、値は，σ＝12m／sの近傍で，

2＝6．85m

で，G二1．36（4）

2＝23．Om

で，

G＝1．26（5）が得られる．

図2の破線で囲まれた測定値は発達した温帯低気圧が日本海を通過するとき，前線に挾ま

3．5

z＝6．85mz≡23．Om

3．O

へSEVERESTO舳

（∵

2．5

＼

＼

2．O

1．5

、。

1．0

051015202530
ラ（m／S）

図22高度での平均風速に対する突風率

Fig．2

Gust

factoragainst

mean

wind

speed

at

two

heights．

れた暖域での強風に対応する．このときGの値は海面近くで大きく，高度23mでは通常の海上
風とほとんど同じ値をとった．上空の擾乱が接地境界層に深く及んだ事例である．右上方の矢印を付したGは高度23．Omにおける突風を示すもので，この場合は上記の事例とは逆に海面近くでは通常のGの値を示している．矢印は時問の経過を指すもので，最大瞬間風速が20m／s程度の状態から突然50m／sまで上昇し，ほとんど同じ強さを保ちながら約50分問続いた．
この問に平均風速が単調に上昇していったため，図に見られるようなGの傾向となった．強い

中規模擾乱が通過したと考えられる．やはり発達した温帯低気圧の南側に存在した暖域で起った．

これらの例は突風率に関してこれまでに提案されている種々の実験式に従わないものである．海洋上では明瞭に現われるため観測が可能となったが，複雑な地形をもつ陸上では現象を抽出するのは困難であろう．乱流のスケーノレを越える強い大気擾乱による突風の存在は，耐風設計において充分考慮されなければならない．

一346一

発達した低気圧による海上強風の乱流特性と長周期スペクトル構造

O．20

内藤

o
o
o
o
■
o o
o
o o ooo

O．15

1b

＼

O．10

oo

o

＿11。。淡1∴一。黙紅■■■ ○■ ． ■
・8．．
・
．
■

O．05

■

■

■
○

呂

；

6．85m

．

z

＝

23．Om

10

20

15

σ1而ノS）

O．20
︒︒

o

o ■○

■
●

。。呼ぎ．

O．15

■

1ミ0．10σ／σ＝0．

．

．。

L．

σ／σ・O．

一
○．

■

O．05

．、

10

■ ■
■○■

20

15

u｛皿／S）

O．10

■

o／u二0

・主．

1ミ0．05σ／7・0

一■■■■■■

。ドー

■
．

10

．

15

20

σ（m／S）

図32高度での平均風速に対する風ベクトル（〃，o，〃）の乱流強度

Fig．3

Turbulent

intensities

of

wind

vectors（〃，ひ，〃）against

height、

一347一

mean

wind

speed

at

two

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

1988年3月

風速3成分（〃，〃，〃）の30分問における乱流強度（σ、／σ，σ、，／σ，σ、，／σ）を平均風速に対して示したのが図3である．（σ、，σ．，σ、．）は各々の風速変動のRMSを指す．図中の2高度の値は同時に測定されたものである．

強風時の風速変動は海面の粗度のみで決定される大気乱流の性質に，乱流よりもスケールの大きい擾乱の挙動が重ね合わされると考えられる．従って図でみられるように乱流強度も平均風速に対してかなりばらついている．しかし海面から上方に離れると3成分ともに乱流強度が小さくなっていく事が分る．この事は突風率の性質と対応する．

強風時の2高度における平均的な強度は図中の直線で示される．各成分とも乱れの強さは
平均風速に依存しないと言えよう．
2＝6．85m，σ≧10m／sのとき，（σ、／σ，

σ、。／σ，

σ．ノσ）＝（O．116，

0，091，

0，050）（6）

2＝23，Om，σ≧12m／sのとき，（σ｝、／σ，σ、，／σ，

σ、．、／σ）＝（O．090，

0，071，

O．036）（7）

σ＜10m／sでは乱流強度のばらつきは各成分とも大きく，見かけ上平均的には大きな乱流強度を持つ性質を示す．上記の測定値は定常状態が非常に長く保たれている海上風の場合と比
べて20％程度大きい．

突風率は水平風速の乱れの強さと密接な関係にあり，Ishizaki（1983）は次の式を提案している．

G＝1＋一1軌＿ 伽 τ，2 σ△τ

（8）△τは評価時問，τは観測時間である．評価時問は充分小さくなると必ずしもGに良い相関をもつとは考えられないが，今回の条件である△τ二〇．24s，τ、＝10minを（8）に代入してGを求めると次の値をとる．

2＝6．85mでG二1．45，
2＝23．Om

で

G＝1．35

上記の値は実測値である（4），（5）とよく一致する．従って定常状態が続く場合，（8）で突風率と乱流強度の関係が近似的に示すことができよう．

海面上の乱流境界層の構造を示すのに最も重要な尺度の一つは，運動量の鉛直フラックス
すなわち風の鉛直シアー・ストレスτ。であり，次式で定義される．
（9）

τo＝一ρ〃〃＝ρ〃柔＝ρC1）乙τ2

ここでρは空気密度である．海面へ輸送されたτ。は風浪を発生，発達させるため，海面の空気
力学的粗度2。が多少変化する．τ。は通常海面の抵抗係数CDでもって評価される．

一348一

発達した低気圧による海上強風の乱流特性と長周期スペクトル構造一内藤

図4に30分問の測達値から求めたCDを2＝6．85mの平均風速σ、に対して示す．シアー・ストレスτ。も接地層の乱流のみでなく，大きなスケールの大気擾乱の干渉を受けると考えられる．

図には台風時の測定値と，発達した温帯低気圧時の測定値を分類して示してある．明らかに台風時のCDの方が大きい．図中の2本の直線は測定値を近似するものである．温帯低気圧の場合の近似線はDavies＆Flather（1978）が地衡風から海上風を求めるために提案した式と一致する．このことは接地層を越えた広域的な見地から与える海上風は，中小規模の擾乱を含んだ「短時問では定常でない測定値」に対応すると考えてよいだろう．

3．O

o

EXTRATROPICAL．

TYPH00NCYCLONE

■

■

2．0

o
x

o

q
1．0

Eq．（1O〕／o、、．（11）／

図46，85mの高度における平均

風速に対する抵抗係数
Fig．4

5

10

15

20

σ（m／s）

Drag

coefficientagainstmean wind speed atheight of6，85m．

the

図の近似線は次式で表わされる．

台風のとき，

103CD＝0，186σユ

温帯低気圧のとき，103CD■一012＋0137σ1

（10）

（11）

但しσはm／s単位の値をとる．上式はσ1＞7m／sの条件下で適用できる．定常的な海洋上の乱流では，Naito（1978）は103CD二0．36＋0，118σ1，＝1．01＋O，036σ1，3＜σ1≦8m／s
8m／s＜ひ1を得ている．上式と（11）を比較すると強風時では（11）の方が平均風速の増大とともに急速にCDが大きくなるが，σ1＝12m／s近傍における値はほぼ等しい．
表1に海洋上の強風の接地層におけるモデル化した値を示す．平均風速σ。＝12m／sに対する台風時と温帯低気圧時の各バラメーターの差異を比較する．表中のαは，平均風速の鉛直分
布をベキ乗則で表わしたときの指数である．（3）で示されるようにαは低抗係数CDから導くことができる．表1のαの値は（10），（11）より求めたCDを（3）に代入して得たものであり，図1から直接求めた値よりも小さい．この理由として考えられることは接地層内の平均風速の鉛直分布に関して対数則とベキ乗則が同時に成立するとした点があげられる．ここで表より

一349一

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

1988年3月

台風時の海面は，他の強風時と比較してかなり粗い事が見出せる．

種々の大気擾乱による強風の10分間乱流統計量を表2に示す．各RUNは長時問強風が続い

表1

強風時の接地境界層における乱流輸送のモデル

Tab1e1

Mode1ofturbu1ent

strongwinds．

transfer

Extratropica1

cyc1oneTyphoonthesurface1ayer

under

CD

〃｝

20

（m／s）

×1000

（cm／s）

（Cm）

12．0

1．52

46

O．02

0．097

（二1／10．2）

12．0

2．23

56

0．13

O．118

（＝1／8．5）

σ1

Case

in

表2

強風時の海上風の10分問乱流統言十Table2

Turbu1ent

winds

under10mininseverestorms．

σω

〃串

CD

ξm

（cm／s）（cm／s）（cm／s）（cm／s）

×1000

（Cm）

8402−7

1，436

2，584

266

184

82

78

2．99

208

O．51

2．60

8402−13

1，46！

2，206

194

162

92

71

2．41

266

0．38

3．17

8414−A118

1，425

2，268

193

126

72

70

2．23

200

O．17

2．94

8414−B26

1，233

66

ll︶

59

2．29

：1：1︶

161

0．04

2．70

150

86

、ll：：）

：lll︶
llll︶

：1；1︶

；：︶

1：︶

1：：︶

：ll︶

226

：1：1︶

265

κ山

：ll：︶

1；；）

ll：：︶

：：︶

1：：︶

1：：︶

：：：；1︶．：lll）

llll︶

282
1ご︶

二︶

ll︶

lll︶

1：1︶

：1：：1︶

8611−129

1：1）

S加

：1：：︶

1，788

：1：；：︶

lll：；︶

8611−119

1：lll︶

8608−41

1：；ll︶

8608−39：：ll：︶

σ〃（cm／s）

RUN

llll：︶

σm（cm／s）

1：：︶σ！。

strong

ll：︶

of

一σ

characteristics

α352

た擾乱の，最も強い10分問平均値が現れた時刻を中心として解析されたものである．従って必ずしも充分定常な状態が続いた条件下で測定されたものではない．しかし風向変化は60似
内のものとし，大部分はこの値よりもはるかに小さい．表の測定値のうちσ朋は最大瞬問風速，∫とκ、は風速〃成分の3次と4次のモーメント，すなわちSkewnessとKurtosisを表わす．ξmは同時刻に容量型波高計で測定された10分問の最大波高を示す．2行にわたって数値がある
RUNでは，上段が高度6．85m，下段が23．Omにおける測定値である．RUN8414−A，Bは台風，

他のRUNは発達した温帯低気圧による強風の性質を与えている．RUN8608−39，41は突風率

一350一発達した低気圧による海上強風の乱流特性と長周期スペクトル構造

内藤の考察でとりあげた図1の矢印を付した測定値の場合を示す．同様にRUN8611−119，129は同図の破線で囲まれた測定値の場合に対応する．この表からも強い中規模擾乱が通過するときの強風は，境界層乱流とかなり異った特性をもつことが分る．定常な海洋上の乱流では風速〃成分のKurtosis

K、は，Gauss分布（K、、＝3．0の値をとる）に近い値を持つことが知られている．すなわちK、、≒2．7であるが，RUN8606−39，41の2＝23．0mでは（下段）κ、、＞5であって突風の形が全くGauss分布と異る統計を示すことが分る．更にRUN8611−119，129の2二6．85m（上段）ではK＞4であり，海面に近い気流のみが全くGauss分布の形をもたない事が見出せる．後者の事例では海面上に定フラックス層の存在を仮定することが出来ない．従って摩擦速度〃。
の意味も通常の乱流と異ったものとなろう．

最大波高ξmは観測塔の位置の水深が平均20mであるため，底面の影響をほとんど受けていないと考えてよい．なお台風通過時の海上風の乱流特性については，内藤（1985）に詳細に報告されている．

3．2海上風の長周期スペクトル
大気境界層の風速のパワースペクトノレ形に関しては，多くのモデルが提案されている．これらのほとんどは規格化された形をとっている．すなわち30分から1時間程度の観測によって得られたデータをもとにして，等方性乱流を基礎とする解析である．換言すればO．1sから10min程度のスケールの乱渦が大気境界層を支酉已しいるとし，乱流よりも大きなスケールの大気擾乱は意識的に，あるいは無意識的に排除している．しかし大気境界層の風の乱れは，境界面に起因する強制対流と大きな規模の大気擾乱が重なったものである．従って特定の周期でもって乱れを二つの集団に分けることができない．

大きなスケールの大気擾乱の挙動が強く影響する周波数帯，すなわち境界層乱流のパワースペクトルのピークよりも低周波数帯でのスペクトル特性はあまり調べられていなくて，報告されているスペクトノレの実験式もK身rman型スペクトル形に準じているのが多い．風速の長期にわたるデータからスペクトル挙動を調べたVander

Hovenの例は，風速変動の大きさに対する解析であって風向を考えていない．これは風をベクト／レ量として扱っていないため，大気擾乱の空間構造を調べる事にほとんど利用できない．

図5に強風時の風速3成分（〃，o，〃）のパワースペクトルの例を示す．高度23mにおける測定で，30分間の平均風速はσ。＝1，769cm／sであった．このスペクトルは典型的な境界層乱流の特徴を示す形を持っていて，3成分の周波数〃に対するスペクトル密度の差異はそれをよく示している．また高周波数側の減衰は，図の直線で示されるように一5／3乗則に従っている．中程度以下の平均風速のときは，高度23m近傍の風速変動は必ずしも良く境界層乱流の特性を示さない．しかしこの例で見られるように，強風時には海面上20mを越える層において充分境界層乱流の性質を備えるようになると考えられる．

一351一

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

1988隼3月

10μ

u
く

へ
（ノヘノ

103

言

・

バ（ノ

．

｛

、
＼

戸ノ
＼ソ

と
q

一、

〃、／

、
上

ユ02

RUN

図5

8611−65

おける風速3成分のバワースペ

U＝1，769cm／s

クトル

2

101

Fig．5

ユO−1

lO−2

強風時における23．Omの高度に

100

101

Power

spectra

of

3

wind

vectorsattheheightof23．0m

η（HZ）

in

strongwinds．

図6に長時問にわたって観測された〃成分のパワースペクトルを示す．解析対象とした比較的定常性が保たれた風速の測定時問Dは四つの事例のいずれも400min以上であり，読み取
り時問は60msである．図中に付した平均風速σエ，σ。はそれぞれ高度6．85mと23mにおけるもので，スペクト／レ形もそれぞれの高度に対応する．ひとユは観測時問Dの問の平均値である．

なお図5は図6の事例の一つRUN8611の測定時の一部分から求めた．RUN8402に付した滑
105

RU
D

8402
＝

・一．帆∴

400m〒n

10叫

1Oヨ

1O叫

デ1

Eq．（7ジσ。＝1・239・・／…

一・．

RUN 8414＿B
D ＝ 720m1n

RUN8608

図6・一1J1

強風が定常状態である条件下における風向方向風ベクトルのパ
ワースペクトル．σ1とσ2はそれ

1〕！

ぞれ6．85mと23，Omの高度での観測時間Dの平均風速である．

820m1nRUN8611

Fig．6、ノ『ノ、・〆〜V

、
D

103

≡

㌦、。

一一
・、、＼
σ・1，330㎝／s 一、

400而in1

口・1，613㎝／s
2
10

一』

10

一ヨ

10

■2

10

−1

10

0

10

1

η（Hz）

spectraconditionofathe1ongitu−

under

the

stationary

wind．σ1，andσ2are

the mean wind speed during
the observation 1） at tw0

heightof6．85mand23．0m
respective1y、

一352一

of

vectorstrong

＼

l02

10■5

Power

dina1wind発達した低気圧による海上強風の乱流特性と長周期スペクトル構造一内藤らかな曲線は低周波数帯の一部を除いて測定値を近似するものであり，次式で表わされる．
〃P。（〃）
σ、2

71〃

（12）

（1＋93〃）5／3

ここでσ、＝188cm／sは30分問の〃成分のRMSを400分問で求め平均したものである．

図の四つの事例をみると低周波数側のピーク近傍は大気擾乱の挙動を示すものであり，この周波数帯のスペクトル密度の大きさが接地層乱流の挙動を示す高周波数帯のスペクトノレ密度に影響を与えていることが分る．RUN8608では明瞭でないが，各スペクトルは5x1O■4Hzから5xlO■3Hzにおいて谷を持つ．この谷の周波数よりも高い周波数帯が境界層乱流が卓越する領域とされよう．表3に谷の周波数〃，及び〃よりも高い周波数帯でのスペクトノレ・ピー
クを示す周波数仰を与える．同時にそれぞれの周波数に対応する波数〃／σも示す．但しRUN

RUN

・。

〃。／σ

（Hz）

（m■1）

（Hz）

〃t

〃t〃

8402

O．O16

0．00130．00055

0．000044

8414−B

0．028

0．00260．00120．000118611

lll：1）lllll：〕llll：：〕llllll：）

（m■1）

表3

図6のスペクトルのピーク及び谷における周波気と波気．
＊印は高度23，0mにおけるもの．

Tab1e3Frequencies
numbers
andoftroughandwave

spectralpeakcorresponding

toFig．6．

8608では求めることが出来ず，RUN8611は2高度でのスペクトルに対して求めており，下段が高度23mの値を示す．表3で示されているピークの波数は，全体として今までに提案されているスペクトル密度のモデルよりもかなり小さい．これは各モデルは陸面上の風のデータを主として解析したものであるため，海面上の風では良く適合しない事を示している．海洋上の風のスペクトル・ピークは陸風に比べてかなり長い周期の方へ偏っているのである．

RUN8402及びRUN8414−Bでの高周波数帯に付した点線は，長い観測時問にこの点線に挾まれた幅を測定値が変動することを示す．2高度の同時測定がなされたRUN8608とRUN8611で見出される重要な事は，低周波数帯では上の高度のスペクトルが大きく，谷を越えた
高周波数帯では下の高度のスペクトルが大きい事である．これはある限界周波数よりも大きい領域では海面の粗度による強制対流が風の乱れの中で卓越するという事をスペクトル密度分布から示しているのである．RUN8611ではその限界周波数が谷の周波数と一致している．

3．3乱渦の鉛直相関
風速変動の空問相関についてNaito（1983）は海上での詳細な実験よりモデル化を行っている．これは接地層の乱渦に対するもので，鉛直方向に大きな高度差があるような場合，この一353一

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

1988年3月

モデルでは必ずしも充分に相関を説明できない．今回の2高度6．85mと23mで同時測定された3成分（〃，o，〃）のうち，〃成分はいかなる強風の条件においても有意な相関をもたなかった．従ってこの節で扱う風速は二つの水平ベクトル〃とoである．高度差16．15mは接地層の乱流のスケールと比べると比較的大きい．そのため小さな乱渦はこの高度差では
てしまう

通り抜け

．また乱渦のスケールは平均風速の増大とともに大きくなるために，弱風の場合も

また相関値が低い．

Shiotaniらは水平風速の相関はかなり大きな時問遅れ（Lagtime）でも有意な相関がある事を指摘している．これは大きな大気擾乱が境界層乱流に重なっているために示される特性と考えられる．

今回風ベクトルの鉛直相関は20分間と200分問の観測値を単位として解析し，長い時問における相関の特徴を主として調べた．2高度及びその高度の平均風速から，代表させる高度と平均風速は幾何平均を行って求めた．従って代表高度は三＝12．55mである．図7に水平風速〃と〃の200分問における相互相関関数Rを示す．相互相関関数の形は200分問に含まれる擾乱の1．0

7・1，016㎝／。！／

＼

／

＼
＼

■

0．5

＼v
＼

1

u＼

／

＼

＼、

1．O

7・1，146㎝／s／／

＼

／

＼

／

串0・5

＼

／

＼

＼

︑口

！

＼7
＼

／
／

＼

＼

■

1．0

7・1，440㎝／s
u

O．5

ノ

／
一

O

／

／

図7

」＼

v

！

観測時問200分での水平風ベク

トル鉛直相互相関関数
Fig．7 Vertical cross・correlation

一

functions−1・000−500

0

500

1，OOO

of

horizonta1wind。。。t。。。。。d。。th．d。。ti。。。fτ（s）

200min．

種類によって異るが，図は代表的な三つの事例を示す．図から見出せる傾向は，風速が強くなると相関はLagtimeτのすべての範囲で弱くなることである．そして中程度の風速でτが1，000sを越えても，〃，o成分ともに有意な相関がある．換言すれば測定時間が200分よりも少なくなると相関値が小さくなる事を意味する．更にほとんどのLagtimeで，o成分の相関値の方が〃成分の相関値よりも大きい．また三つの事例ともLagtimeτ＝Oの中心線に対して対称な形をもっていなくて，擾乱の空間的特徴の一端が見出される．相関関数のピークはτ＝2s前

一354一

発達した低気圧によろ海上強風の乱流特性と長周期スベクトル構造

内藤

後のときに現われるため図では詳細に示すことが出来ないが，τ＝Oでの相関値とピーク値は
かなリ異っている．

図8に観測時間200分の相関係数と20分の相関係数R、を平均風速σに対して示す．図の値
はいずれも相関関数のピーク値である．図より観測時問200分の相関の方が20分の相関よりもかなり大きい事が分る．これは図7の相関関数の形からも示唆された．測定値のばらつきもかなり大きいが，〃，o成分の両方とも平均風速の増大とともに相関が小さくなっていく傾向にある．これは20分問の相関の方が低い風速から始まる．
1．O

＼

＾
・■o＾ ・
＼
o一．
■
6ら＼・
＼ ■
o o ＼、、；、一。。、1、 の哨。。▲ε＼＼

6

O．5

＼

1：1：：：ll∵〜淳ぺ：：1：ll；
v；

1〕

＝

0

200m1n

5

10

15

20

25

7（m／s）

図8
Fig．8

平均風速に村すろ200分問と20分間の水平風ベクトルの相関係数Correlation

coefficients

ing200minofand20minhorizontal

against

mean

wind

wind

vectors

dur−

speed．

Panofsky＆Singer（1965）は鉛直相関係数を二つの高度の関数として次の形で表現した．
1〜、＝expl一（22U3−2I1／3）1，

i＝〃，

ひ

ここで2。＞21をとる．彼らは相関係数の平均風速に対する依存性がないと仮定し，した．実測よりα

（13）

を定数と

＝O．5，α、，＝0．86を得ている．図8で明らかなように今回の測定より相関係

数R、は平均風速σと観測時間Dの関数であるとするのが極めて合理的である．兄のDに対する依存性は今回の二つの条件からでは数式化できないが，σに対する依存性を（13）の形を踏襲して行う．すなわち減数定数を改めてσの関数とする．図中の近似線は〃成分に対するもので次のように表わされる．

D＝200min，

σ＞8m／sのとき

α、、＝O．076（σ一8．O）
D＝20min，（14）

σ＞5m／sのとき

α、、＝O．122（0

一5．0）

（15）

ここでσはm／S単位の値を用いる．o成分の相関は〃成分に比べてわずかに大きい事が図より
見出せる．

一355一

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

1988年3月

乱渦の鉛直相関の周波数スペクトルを表わすコヒーレンスγは次式で定義される．

G身）寺絆一ππ…／lφ（い）／（16）

ここでハ（〃），乃（〃）は高度差／の2点A，Bにおける風速変動のバワースペクトルであり，C・（〃）十クQ（〃）は同じくクロススペクトルである．またφ（1，〃）はAB問の位相の遅れを表わす．Naito（1982）は接地層乱流の3成分について3次元空間でのコヒーレンスを，無次元周波数鵬／σと無次元距離／／2の関数として表現した．このうち鉛直方向についてのみ示すと

γ（÷，害）一…／一ん（÷）一箒／・（17）（λ、、，λ、、，λ砒∫）＝（24．O，12．5，8，8）（18）

である．Davenportが最初に示した実験式は，（17）の／／2の指数1．26を1．0に置き換えたものである．従ってDavenportはコヒーレンスの高度依存性はないとしている．

空問距離1が大きい場合，小さいスケールの高周波数乱渦が相関をもたなくなるため，（17）ではコヒーレンスを充分に表現できなくなる．また接地層乱流のスケールを越える擾乱に対しても（17）で表現できる限界を調べる必要があろう．今回，2高度6．85mと23．0mでの風速変動における強風時の乱渦のコヒーレンスについて吟味し，大きなスケールの乱れの挙動を調べる．

〃成分については相関関数においても指摘したと同様に，すべての周波数について全く有意な相関をもたない．従ってコヒーレンスは水平2成分〃とoについてのみ評価する．〃，o成分のコヒーレンスの実測値を観測時問30分程度で求めると，周波数依存性はこれまでの多くの解析の前提となった指数関数的減衰を示さない．今回の観測での高度差1＝16．15mは接地
層乱流の挙動を調べるために大きすぎるのである．しかし観測時間を充分長くとれば大きなスケールの擾乱をとらえることができ，コヒーレンスが周波数の増加とともに単調に減少する形が期待できる．

図9に観測時間200分のコヒーレンスの例を示す．図中の曲線はコヒーレンスが周波数〃とともに指数関数的に減衰するとし，〃＜0，05Hzの〃成分の測定値に対して近似したものである．かなりばらついているが3例とも周波数の増加とともに単調に減少していると言えよう．

しかし〃が充分大きくなってもコヒーレンスは0に充分近づかず，かなり広い周波数帯で有意な値を持ち続ける．そしてこの領域では〃成分のコヒーレンスよりo成分のコヒーレンスの方が大きいことが見出せる．これは以前も指摘したように接地層の乱渦の特性である．
図から分るように，〃＞O．05Hzではコヒーレンスは指数関数で良く近似できない．他の多一356一発達した低気圧による海上強風の乱流特性と長周期スペクトル構造

内藤

1．O

も

O．5

策

σ

o、

．

＝

1，110cm／s

．古・・

o

o

㌣。；8θ．：1・g、・1；1：・

1．O

塾。一1…㎞／・
・※。㌔
圭O．5

8

箏・○・．・。．・

圭

。い。。

崖
掌1．O

・

・・；88。

。．

8

羨び＝1舳
O．5

o

■

O、斉舳二．．1・：
o■

O

O

o■

■o■ ooO．05

図9200分問における水平風速ベク
トルの鉛直コヒーレンス

Fig．9

Verticalcoherenceszonta1wind

0．10

vectorsofhori−during200

min．

パHz）

くの観測例から，風速が強くなるとともにこの傾向が増大することが分る．従って（17）でもってコヒーレンスを定式化することが困難であるが，長周期の大気擾乱の鉛直相関を調べるために比較的指数関数的な減衰を示すRUNについて同様な方法で解析した．図10にコヒーレ
ンスγが周波数〃とともに指数関数的な減衰をするとき，その減衰定数一z閉γ／〃を平均風速σに対して示した．図中の2本の曲線は，（17）でAに24及び31の値を与えて更に今回の測定条件を代入し減衰定数の風速依存性を示したものである．各々の曲線は〃成分及び〃成分の測定
値を（17）で近似したときの様子を示す．全件的に測定値はばらついていて，平均風速に対す
80

㌧ふ

？60
詣

■

ミ叫0

■

ト
ミ

λ≡31

20

図10

λ＝2叫

u

平均風速に対するコヒーレンスの指数関数減衰定数v

Fig．10Exponential
5

10

15

20

25

石（皿／・）

一357一

decay

eddy

coherencewind

speed．

factoragainstofmean

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

1988年3月

る依存性も明瞭でない．しかしひ成分の値の方が〃成分の値よりも小さいことが分る．すなわ
ちo成分の方が〃成分よりも相関が強い．

図中の曲線は測定値の特徴を表現するのに必ずしも適当な近似でないことが分る．しかし
境界層乱流の空間相関を表わす（17）と今回の結果を比較するため，敢えて図示した．（18）で与えられる定数λとλ王，に比べて図中の曲線で与えられる定数λはかなり大きい．従って境界層乱流の特性が，乱流よりも充分大きな大気擾乱において持続されるとして適用した（17），（18）の乱渦のコヒーレンスよりも実測値はかなり小さいと言える．ここで再び言えることは，接地層にまで大きなスケールの擾乱が及び，その擾乱は乱流に重なって境界層の諸特性にばらつきや偏りを与える．しかし大きな擾乱は乱流とは基本的に異質なものであると思われる．

4．結論

長時問続けて吹いた海洋上の強風を観測し，統計的に解析することによって接地境界層の乱流特性及び乱流のスケールを越える大気擾乱の挙動を調べた．3次元風ベクトルは海面上6．85mと23．Omの2高度で同時に達続測定した．主な結果は次のとおりである．

（1）海上風の突風率は，通常の境界層乱流では高度約7mから23mの層において平均的に1．46から1．36の値をとる．他方発達した低気圧に伴う大気擾乱が通過するとき，同層において3．2程度までの大きな値をとるようになる．

（2）強い大気擾乱が境界層乱流に重なるとき，強風の乱れの強さや海面の抵抗係数は定常な乱流の状態で測定された値よりも平均的に少し大きくなり，ばらつく．更に大気擾
乱の種類によって平均的な値が異ってくると考えられる．

（3）測定時間が400分以上にわたる広周波数領域の風のパワースペクトルを求めた．10分を越える長周期のスペクトル挙動は強風の原因となった擾乱によるものと推定され，一つのピークを形成する．

（4）大きな高度差をもつ2点での相関係数は，平均風速の増大とともに減少する．

（5）大きな高度差をもつ2点での乱渦のコヒーレンスは，長時問の測定値でも周波数の増大とともに必ずしも指数関数的に減少しないし，また平均風速に対する依存性も明瞭でない．従ってDavenportの相似則では接地層乱流のスケールを越える擾乱に対して良い近似で表現できない．

謝辞

平塚沖の海上観測塔における測器の設置等の準備に関して，海洋科学技術センター佐々木
保徳氏に助力をいただき感謝いたします．

一358一

発達した低気圧による海上強風の乱流特性と長周期スペクトル構造

参1）Davenport，A．G、（！961）：The

spectrum

考

文

内藤

献

of1ongitudina1gustiness

near

the

gromd

in

high

winds．Quart．

J．Roy．Meteorol．Soc、，87，194−221．

2）Davies，A．M，andContinental

R．A．Flather（1978）：Application

Shelf

and

the

North

Sea

to

the

of

numerica1mode1s

computation

of

the

of

storm

the

surgesNorthwestofEuropeanNo▽ember−December

1973．Dtsch．Hydrogr．Z．Eng．H．A．14，1−72，

3）Hino，M．（1971）：Spectrum

of

gusty

wind．Proc，3rd

Int．C㎝f．㎝Wind

Effectson

Bldgs．and

Structures．Tokyo，1971，19−7！．
4）Ishizaki，H、（！983）：Wind

regoins．J，windprofi1es，turbulent

E㎎．andintensities

and

gust

factorsfordesignintyphoon・prone

characteristics

ofsurface1ayer

Ind．Aerodynamics，13，55−66．

5）Kaimal，J．C．，J．C．Wy㎎aard，Y．Izumi

and0．R．Cote（1972）：Spectral

turbu1ence．Quart．J．Roy．Meteorol．Soc．，98，563−589．

6）Kondo，J．，Y．Fujinawa

andG．Naito（！972）：Wave−induced

wind

fluctuation

over

the

sea．J．Fluid

Mech．，51，751−771．

7）Naito，G．（1978）：Direct
open

measurementsofmoment㎜andsensib1eheat

fluxesatthetowerinatmosphericinthesea，J．Meteoro1．Soc，JPN、，56，25−34．

8）Naito，G．（！982）：Three−dimensional
layer

above

the

space

structureofturbu1ent

eddy

the

bomdary

ocean．J．Meteorol．Soc．JPN．，60．1299−1315．

9）内藤玄」（1985）：台風通過時の海洋風の乱流特性．国立防災科学技術センター研究報告第35号，315−333．
10）Panofsky，H．A．and

I．A．Singer（1965）：Vertica1Structure

ofturbu1ence．Quart．J．Roy．Meterol．Soc．，

91，339−344．

11）Pie1ke，R，A．and
Layer12）Shiotani，M．andstrong13）VanH．A．Panofsky（1970）：Turbu1ence

characteristics

a1ong

correlationsofseveraltowers．Boundary−

Meteoro1．，1，115一ユ30．

der

0007towind．J

Y，Iwatani（1976）1Horizonta1space

ve1ocity

f1uctuations

during

Meteoro1．Soc．，54，59−67．

Hoven，I．（1957）：Powerg00cyc1es

perspectmmofhorizontal

wind

speed

in

the

frequency

range

from0．

hour．J．Meteoro1．、14，160−164．
（1987年11月24日

359

原稿受理）