JPS59501174A - 水中地震探査における音源特性のスケ−リング - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるため要約のデータは記録されません。

Description

【発明の詳細な説明】 水中地震探査における音源特性のスケーリング本発明は水中（例えば海洋）地震 探査での音源特性のスケーリングに関する。本発明は、アントニ　マーシャン　 ジオルコラスキー名義の英国特許用＠８０１３４３８　（公開番号２０４８４８ １）　及Ｕ　８０１３４３７　（公開番ｌ’ｔ　２０４８４８０　）に記載及び クレームされた発明を敷桁改良するものである。本願は特に、本発明の利用に必 須の空気銃等の初期発破圧が制御可能である特定の特徴を有する音源に関する。

参照用に内容をここにとり入れた係属中の出願公開番号ＧＢ　２０４８４８１て は、それぞれの場所で地震探査法を２回、発破及び受信器の位置を変えず（ご行 なう。但し２回目の発破は１回目の発破をスケーリングしたものを行なう。２つ の受信された震動記録ＸＩ（ｊ）及びｘ２（ｔ）、大地のインパルス応答ｇ＜ｔ ）は発破及び受信器の位置が変わらないがら同一である。しかし、２つの源の特 性Ｓ、（ｔ）及び５ｚ（ｊ）は異なり、源のスケーリング則により互いに関係す ることが示される。よって、次の３つの式による関係がある３つの未知の１１　 Ｓｔ　（ｔ　）、Ｓ２　（［）及び＋７（ｔ）が存在する。

Ｘ＋　（ｊ　）−３ｔ　（ｔ　）＊　（］（ｔ　）’　（１）Ｘ２　（ｔ　）＝ 　５２　（［）＊　ｃａ（ｔ　）　（２）Ｘ２（ｔ）＝αＳ＋（ｔ／α）　■ ここで星印はただみこみを示しスケール定数αは次の関係で決められる。

浄書（内容に変更なし） 式（１）及び■では雑音を無視している。これらの式は、上記出願に記載された 雑音があっても安定な粗い方法によりＸ＋（ｉ）及びｘ２　（Ｈの測定からＳｔ 　（ｔ　）、Ｓ２　（ｔ　）及び（１（ｔ）を得るように解きうる。

上記出願の発明の要点は、以前の方法では式（１）の如くλ２Δ未知量を含む１ つの式を与えうるのみであったということである。Ｓｔ（ｊ）についての正確な 情゛報なしにｇ（ｔ）を正確に決定するという問題が常にあった。出願公開ＧＢ 　２０４８４８１における如（実験を繰り返しても、スケーリング開式■により 付は加えられる情報がなければ意味をなさない。

このスケーリング則を適用するのに変更しうるただ１つのパラメータは点源の友 」劃だけである。例えば、爆薬を海で使用する場合、化学成分は源１と源２では 同一でなければならず、海面下の深さも同一でなｃブればならず、大きさのみが 異なりαによる関係を有する（第２の爆薬の質ｌは第１の爆薬の質量のα３倍で ある）。空気銃が使用される場合、それらの容積は係数だけ変化するが深さ及び 発破圧力は同一でなければならない。

第２の出願公開ＧＢ　２０４８４８０では、例えば空気銃の列等の点源の２次元 的な列に思想が拡張されている。１つのスケーリングされた列の各点源素子は、 他の列では対応する素子を単にスケーリングしたものであることか提案されてい る。つまり、列２と列１との対応する素子のエネルギの比はα３　（式４による ）である。２つの列の大きさは対応して係数αによりスケーリングしている。よ って、列２の水平方向の寸法は列１のそれのα倍である。他のパラメータは変更 されない。深さは同一であり、（空気銃の）圧力又は（爆薬の）化学成分は周一 である。これらの条件下では式０は、２つのスケーリングされた列から所与の距 離と所与の方向において発生した波Ｓ２　（ｔ）及び５ｔ（ｔ）に関係する。こ れは参照のためここに内容を取り入れた出願公開ＧＢ　２０４８４８０に充分に 説明されている。

２つの先願の問題点 先願の解法には海面からの［ゴース１へ反ｍｌに関係する２つの問題点がある。

源の場合、海面は反射係数か略正確に−１である鏡の如く動く。それには、第１ 図に示す如く実際の源とは逆の極性の音波の仮想地震源を創り出す効果Ｄ１ある 。点Ｑに達する音波は、音源からの直接の波動と仮想源に由来するように見える 反射波又は［ゴーストｊの両方の波動の和である。

如何なる源から発生される波動も、それ自身の直接の波ｉ０ζそのゴースト反則 との和として考えるのが便利である。

式（１）及び■て表わされる震動記録についての桟々の概念では、このゴースト 反則を源をの小波動Ｓ＋（［）又はＳ２　（ｉ）とは別のものとして扱っていた 。我々はそねをｇ（ｔ）に属するものとし、て扱ってきた。

従って、式（１）からＳｔ（ｊ）をとり出すとこのゴース［−を含む９（ｔ）が 得られる。同様に式のから３２（［）をとり出すとこのゴーストを含む別の評価 の（１（ｔ）が得られる。従って、９（ｔ）は我々が定義した如き大地のインパ ルス応答ではなく２つの異なる場所での２つの〜でンパルスに対する応答になっ てしまう。第１のインパルスは時刻１＝０に源の位置で発生する。第２のインパ ルスは同時刻に発生が、仮想源から逆の３性で発生するように見える。本当のイ ンパルス応答を得るには、この第２の「ゴースト応答」を除かねばならず、これ は非常に不便である。

リングしないため不可能であった。

このことを餉１．２図を参照してみる。第２図は点Ｑに達する直接波動ｓＢ）及 びその「ゴースト」〜５（ｔ−ｔ、）を示す。２つの波動の和は Ｐｏ　（ｔ　）＝　Ｓ（ｊ　）　−Ｓ（ｊ　−ｔ　＋　）　（５）である。こ１ て時間遅延は ｔ＋＝２Ｄｃｏｓθ　（６） であり、負の符号は高面の反射係数か略正確に−１であるからである。

２つの異なる源がＦＡ　−（７）深さにある場合、直接波動５ｔ（ｊ）及び５２ （１）は式■の如くスケーリングするが、ゴーストの遅延１．は変化しないため 直接波動とゴーストとの和はスケーリングしない。このことは第３図に示しであ る。

問題点２ 源とその仮想像又はゴーストは、それらがあったかも略独立しているかの如゛く 説明されてきた。勿論ゴーストは知られている如く源及び表面に全面的に依存す る。しかし、源はそれ自身のゴーストから完全には独立しておらず、完全には無 視しえない一種の帰還効果がある。

例えば空気銃の場合、音波は銃の作る気泡の振動により発生される。

第４図を参照するに、気泡の内側と外圧つまり静水圧との間には圧力差があるこ とがわかる。初期では内圧は典型的には約１３５気圧であるが静水圧は約２気圧 である。従って気泡は急速に水を押し拡げ、内外圧が等しい平衡点を越えて減速 し次いて停止する。その時気泡は最大の寸法に達しており内圧は静水圧より小さ い。今度は圧力差は内方に働く。気泡は縮小し、縮小は気泡内の圧力が急速に上 昇することによってのみ停まる。この点から振動が再開する。音波の源となるの は周囲の水を押すこれらの気泡の振動である（ラウドスピーカが周囲の空気を押 すことで有効な音波源となるのと略同−の仕方である）。

ゴースト反射は、振動する気泡を越えて伝わる圧力波である。それは周囲圧力に 重畳される圧力の乱れとなる。つまり、それは水中の静止圧Ｐ＋に動態的な項ｐ （ｔ）を付は加える。よって気泡における水圧は一！ 定ではなく時間とともに次のように変化する。

Ｐｗ（ｔ　）　＝ＰＨ＋Ｄ、、　（口　（ｎここでＩ”ｖ（を目は時刻ｔでの全 水圧、ＰＨ（は静水圧及びｐ、（ｔ）はゴースト反則による動態項である。通常 の動作ではｐ　（Ｉｌ）はＰＨと比べ！ で小さいからＰｗ（ｊ）は相当に一定である。しかし効果は存在し無視し得ない 。

点源にしろ列にしろ空気銃をスケーリングする場合ＰＨは両方に対し同一である 。しかし、０（１）は源から発生した波動の反射であるからスケーリングした音 源に対しては異なり、従ってＰｗ　（ｔ　）はスケーリングした音源に対しては 異なる。このため気泡振動の動態はスケーリングした源に対しては異なり、スケ ーリング開式０は正確には成り立たない。

本願の目的はこれら２つの問題点を部分的に又は完全に解決することにある。

五艮肚 本発明の要部は、αを特性のスケーリングに所望なスケーリング係数としく式０ 参照）、スケーリングされた深さにおけるスケーリングさ机た源の動態的性質を 略一定に維持しつつ源２の深さを源１の深さのα倍に等しくなるようスケーリン グすることにある。

本発明によれば、 ａ）　所望の場所の所定の水深で所定のエネルギの弾性放射を有する第１の音波 を所定の初期圧力で発生する第１の音源を使用し、次いて上記第１の音源と同種 であり同一の場所だが第２の異なる所定の水深で第２の異なる所定のエネルギの 弾性放射を有する第２の音波を第２の異なる所定の初期圧力で発生する第２の音 源を使用し、深さ、エネルギ及び初期圧力は。

第２の音源から出される弾性放射のエネルギは第１の源から出されるそれのα３ β３４倍であり、それぞれの音波の周波数スペクトラムは重畳第２の源の深さＤ ２は第１の源の深さＤｌのα倍であり。

第２の源から出される音波の初期圧力Ｐ２は第１の源から出される音波の初期圧 力Ｐ１と式 による関係かあるよう選定されており、ここでＰＨ＋は第１の源での静水圧を表 わし、ＰＨ２は第２の源の静浄よ、内゛Ｇに変更なし） 水圧を表わし、また ＩＩ）　１又は複数の受信器により第１及び第２の地震信号をそれぞれ発生する よう音波の大地内からの反射を検出し。

Ｃ）　これらの信号を、遠距離の源の小波動及び雑音が略なく水面からの音波の 反射による誤差を略補償して大地のインパルス応答を導くよう分析比較して。

表面上境界の大地中の位置及び／又は大地の表面下の特徴の音響学的特性を決定 する方法か提供される。

第１及び第２の音源は、それぞれ点音源でも又はそれぞれ互いに関係する点音源 のケ］でもよい。後者の場合、第２の音源を形成する列は、第１の音源と同一の 数と型の点音源を含み、同一の空間的関係を有するが第２の源の列の各点音源間 の間隔は第１の源の列の対応する間隔をα倍にスケーリングしたものでなければ ならない。列の点音源は全て共通の深さになければならない。

α及びβの値は１．１から３が適当であり、１．５から３であるのが好ましい。

一般的に上述した圧力関係が満たされるよう銃の初期発破圧力は調整可能でなけ ればならないから点音源として空気銃を使用するのは好ましいものであるが、初 期圧力が調整可能な他の源を使用してもよい。

空気銃源の場合深さが固定されると可変しうるパラメータは初期内圧つまり発破 圧力及び銃の容積の２つである。銃の容積は動態にはあまり影響しない。実際上 述の如（容積が変わっても深さ及び圧力が一定に保たれるなら、αが容積比の立 方根に等しい式３に示されたスケーリングを別にして音波の形状は同一に保たれ る。このスケーリングは撮動する気泡の一定の動態の表われにすぎず、スケーリ ング係数は気泡の大きさにのみ関係して動態には係らない。上述（問題点２）の 如くゴーストによる小さな容量依存性の帰還効果は包ヲが、動態には僅かにしか 影響しない。しかし他の点では容積が気泡の動態に影響しないというのは真実で ある。

帰還効果を別にすれば気泡の動態は深さ及び内圧にのみ依存する。深さが固定さ れると、変動する状態を補償するよう変更できるパラメータは内圧だけである。

上述の如く気泡を振動させて音波を放射せしめるのは気泡璧の両側の圧力１であ る。従って、動態を維持する場合には気泡の内圧は静水圧の変化に比例してスケ ーリングされねばならない。

撮動する気泡の動態の綿密な分析により、次の仮定がされるならば銃間及び銃と ゴーストとの間に相互作用のある空気銃列に対してもスケーリングが可能である ことが示される。

（１）空気銃の作る各気泡は、発生する音の波長に比べて小さい。これは通常の 条件下で動作する空気銃についての経験的事実である。

［Ｆ］　地震波に変換される全有効エネルギの割合は２つのスケーリングされた 源に対し同一である。こ机は出願分間ＧＢ　２０４８４８１における仮定と全く 同一である。

源のスケ−１ルグを達成するために満たされねばならない気圧は、ゴーストの相 互作用もスケーリングしなければならないということである。

これらの条件下では、点源及び点源の列に対し正確なスケーリングが得られる。

まず点源のスケーリングを説明し、次いで列のスケーリングを７゛書／１″′’ ３１：ａ＊１：　Ｌゝ　慈ｏａ５９−５ｏｎ７４（４）説明する。

点源のスケーリング（空気銃） 容積Ｖ＋、初期発破圧力Ｐ１及び深さＤ＋の空気銃について考える。

容積Ｖ２．深さ０２．初期圧力Ｐ２で第２の銃により発生される音波ｌそのゴー ストが第１の銃により発生される音波及び１左ゴーストをスケーリングしたもの である如き銃を発破することが望まれる。

時間のスケーリング係数は深さの比と同じでなければならない。よって ２つの銃における静水圧は である。ここで、Ｐｏは大気圧、ρは水の密度９ｇは重力加速度である。

静水圧の比をβと置くことにする。

であるが、Ｄ２及びＤｌは式（８）による関係があるがらβはαがら独立してい ない。容積及び初期発破圧の両方の比が正しく選定された場合に正確なスケーリ ングが起こる。

１れ・（内容に変更なし） スケーシングされた源は、それぞれ深さＤｌ及びαＤ１でそれぞれ静水圧ＰＨ１ 及びβＰＨｚでありそれぞれ内外圧がＰ＋　（ｊ　）、　Ｐ２　（αｔ）及びＰ ｗｚ　（ｔ　）、ＰＷ２　（α（）である仮想上は球形の振動する気泡と考えら れるものとする。これら２つの仮想上の気泡は同一の無次元の運動方程式に従う とすると、それらの半径の関係はである。よってｔ＝Ｑとすると である。気泡の体積Ｖと半径との関係はＶ＝４／３πＲ３（＋４） であるから、２つの仮想的気泡の初期容積の比はスケーリングされた気泡の壁の 両側の圧力差には式（１６）の関係があるａ− マー？＋ττ（内容ｊこ変更なし） 得られた式を、線型の仮想的な気泡から非線型の実際の気泡へ単に移すことがで きる。つまり深さ及び静水圧の比は線型の気泡に関し残るが、初期気泡容積の比 １式（１５）及び圧力２式り１７）もそうである。

式（１５）及び（１７）の関係の下では第２の銃により発生される波の場は、第 １の銃により発生される波の場をスケーリングしたものである。

第５図を参照するに、第１の銃に関して座標位置（ｒ、θ、φ）の点Ｑ１で受信 される信号ＰＱ１　くｔ）は、第２の銃に関して座標位＠（αｒ、θ、φ）の対 応点Ｑ２で受信される信号Ｐｏｚ（ｊ）に関係する。

Ｐｏ２　（ｔ）とＰＱＩ　（［）との関係はである。距離ｒか源とそのゴースト との距離に比べて大きい場合、つまりｒが２Ｄ＋に比べて大きい場合 ｒ＞＞２ｏ　（２０） 音波の形は距離が増えても変わらない。しかし、波面領域は波動が源から伝播し ていくに従い増大するので、強度は距離ｒに反比例して減少する。この基準が成 り立つ領域は源の「遠距離用」として知られており、そこでは「源」はこの場合 １まとまりとして考えられる空気銃叉工そのゴーストからなる。この遠距離用条 件が満されている場合には及び 浄書（内容二こ］五【７．−ノ とすることができる。ここでＣは音速でありθ及びφは第５図に示した信号の方 向を特定する。遅延時間を とすると、深さがスケーリングされた点源の遠距離場放射に対するスケーリング 則 Ｓ２　（θ、φ、ατ）＝αβ　Ｓｔ　（θ、φ、ｒ）　（２４）が得られる。

あるいは、より簡単には角θ及びφで定められた源からの所与の方向では所与の 距離での遠距離場特性には次の式による関係がある。

このスケーリング則は点源のスケーリングに使用されていたものに非常によく似 ている。異なっているのは、気泡の振動の動態を異なる静水圧で同一に維持する 必要から出て来た係数Ｂだけである。係数βは静水圧の比である。

座標（η、、”ｉ、ｚｉ）を有するｎ個の空気銃の列を考える。深さの座標は水 面に関し、スケーリングが可能であるよう列中の全ての銃につき同一でなければ ならない。従って 待表昭５９−５０１１７４　（５） Ｚ＋　＝　Ｚ２　＝　Ｚ３　＝−ＺＴＩ　＝Ｚ　（２６）列の銃は異なる所定の 時刻 ｔｌ、ｔ２．ｔ３・・・、ｔＴｌ＜２７＞に発破される（勿論銃が同時に発破さ れるよう設計された場合は発破時刻は等しく １、　＝　１２＝・・・−ｔη であるが、これはスケーリングのために必須では虹）。

列久更そのゴーストにより発生された彼の場が第１の列及びｉΔゴーストにより 発生された波の場をスケーリングしたものであるようなスケーリングされた空気 銃列を創ることが望まれる。単一の空気銃とそれ自身のゴーストとの間の相互作 用つまり帰還と全く同じに列中の続開には相互作用が起こる。相互作用のスケー リングを確実にするためには列切寸法をスケーリングしなければならない。スケ ーリングされた列を表示するのに添え字のＳを用い、時間伸張係数をやはりαと すると、スケーリングされた列のｊ番目の銃の座標は（ｘ　、ｙ　、ｚ　）であ る。こＳＩ　ＳＩ　Ｓｌ こで スケーリングされた列の６銃における静水圧は、浄書（内容に変更なし） ＰＨ＝　ＰＯ＋ＤＯＺ　（３０） である。列の素子の単一の空気銃のスケーリング則から得られる。

スケーリングされた列のｎ個の銃の容積は銃の発破時刻は 第７図は実施の陣列が如何に見えるかを示す。出願公開ＧＢ　２０４８４８０に おける如く列は、第２のスケーリングされた列に対する受信器位置が第１の列と 同一となるように、受信器群間隔の倍数に等しい距離りだけ前後に離れて牽引さ れてもよい。列が相互作用をするほど近づかないならば複数の第１の列及び複数 の第２の列を使用することもできる。かがる第１の列のどれか２つの間に相互作 用が生じたならば、それらは単一の列として上記の基準により設計された第２の スケーリングされた列として処理されるへきある。

複数の干渉しない列が使用されるならば、第２及び第１の音源に対し同数の列が 、第２の音源の列のそれぞれは前）ホの如く第１の音源の対応する列に対しスケ ーリングされているようにして使用されるべきである。

列は所望の音波を発生するよう同時に発破されてもよいし、順次発破されて発生 した地震信号を所望の地震信号を得るよう加算されてもよい。

継時的発破及び同時発破を組合せて使用することもできる。

列が出来あがったならば、使用法は深さに依存する余分なスケーリング係数βを 含むスケーリング則 Ｓ２　（τ）−αβＳｔ　（τ／α）　（３５）を除けば、２つの先願公開ＧＢ 　２０４８４８１及びＧＢ　２０４８４８０に記載されたのと同じである。

この全体的スケーリング開法は２つの点で重要である。第１に、５１（１）又は Ｓ２　（ｔ）の位相スペクトラムについて何らかの仮定をする必要又は大地のイ ンパルス応答ｇ（ｔ）の統計的性質についての仮定をする必要さえない。第２の 、源が正確にスケーリングされているならＳｌ　（ｔ）又はＳ２　（ｊ＞を測定 する必要さえない。２つの独立した発破記録から収集されたデータ及びスケーリ ング則から任意の方向の大地のインパルス応答を直接計算することができる。

同時係属出願における如く、弾性エネルギは通常それぞ机何らかの適当なジオホ ン又はハイドロホンを使用して粒子速度又は音圧として検出及び測定さ机る。

手続補正書は） 昭和５９年４月１１日 特許庁長官　若　杉　和　夫　殿 １、事件の表示 ＰＣＴ、／ＧＢ　８３１００１５９ ２、発明の名称 水中地震探査における音源特性のスケーリング３、補正をする者 特許出願人 住所（居所）　イギリス国、ロンドン　ダブリュー１ワイ８ビーエヌ、ポルトン 　ストリート　２９番地名称　プリトイル　ビーエルシー 代表者　ケリー、パトリツク　マイケル　トノパン（国籍　イギリス国） ４、代理人 住所　〒１０２　東京都千代田区麹町５丁目７番地６、補正の対象 特許法第１８４条の５第１項の規定による書面。

明１８１の翻訳文及び請求の範囲の翻訳文、及び代理権を証するもの。

７、　補正の内容 （１）特許法第１８４条の５第１項の規定による書面の記載を別紙のとおり補正 する。

■　明細書の翻訳文中、第２頁、第７頁、第９頁乃至１２頁、第１４頁の多層（ 内容に変更なし）を別紙のとおり補充する。

（３）請求の範囲の翻訳文中、第１６頁の浄害（内容に変更なし）を別紙のとお り補充する。

（４）委任状及びその訳文各１通を別紙のとおり補充する。

田捺調査報告

Claims (1)

【特許請求の範囲】
1. １．８）所望の場所の所定の水深で所定のエネルギの弾性放射を有する第１の音 波を所定の初期圧力で発生する第１の音源を使用し、次いで該第１の音源と同種 であり同一の場所だが第２の異なる所定の水深で第２の異なる所定のエネルギの 弾性放射を有する第２の音波を第２の異なる所定の初期圧力で発生する第２の音 源を使用し、深さ、エネルギ及び初期圧力は。 第２の音源から出される弾性放射のエネルギは第１の源から出されるそ机のα３ βＶ町＠てあり、それぞれの音波の周波数スペクトラムは重畳しており。 第２の源の深さＤ２は第１の源の深さＤＩのα倍であり。 第２の源から出される音波の初期圧力Ｐ２は第１の源から出される音波の初期圧 力Ｐ１と式 による関係があるよう選定されており、ここでＰ＋−＋ｔは第１の源での静水圧 を表わし、ＰＨ２は第２の源の静水圧を表わし、また ｂ）　１又は複数の受信器により第１及び第２の地震信号をそれぞれ発生するよ う音波の大地内からの反射を検出し。 Ｃ）　これらの信号を、遠距離の源の小波動及び雑音が略なく水面からの音波の 反射による誤差を略補償して大地のインパルス応答を導くよう分析比較して。 表面上境界の大地中の位置及び、′又は大地の表面下の特徴の音響学的特性を決 定する方法。 ２　該第１及び第２の音波源はそれぞれ１つの点音源からなることを特徴とする 請求の範囲第１項記載の方法。 ３、該第１及び第２の源はそれぞれ相互影響する点音源の列からなり、該第２の 源の列は第１の源の列と同一の数と型の点音源を含み同一の空間的関係を有し、 第２の源の列の各点音源間の間隔は第１の源の列の対応する点音源に関し係数α によりスケーリングされており、第１の源の点音源の発破時刻は第２の源の対応 する点音源に関し係数αでスケーリングされており、該第の点音源は全て同じ深 さにあることを特徴とする請求の範囲第１項記載の方法。 ４、第１の複数の非干渉的な列が第１の源として使用され、各列による地震信号 は該第１の地震信号を発生するよう加算され、それぞれが該第１の複数列の対応 する列に関しスケーリングされている同数の非干渉的な列が第２の源として使用 され、各列による地震信号は該第２の地震信号を発生するよう加算されることを 特徴とする請求の範囲第３項記載の方法。 ５、第１の複数の非干渉的な列が第１の源として使用されて該第１の音波を発生 するよう同時に発破され、それぞれが該第１の複数列の対応する列に関しスケー リングされている同数の非干渉列が第２の源として使用され、該第２の列は該第 ２の音波を発生するよう同時に発破されることを特徴とする請求の範囲第３項記 載方法。 ６、各点音源は空気銃からなることを特徴とする請求の範囲第１項記載の方法。 ７、各点音源は空気銃からなることを特徴とする請求の範囲第３項記載の方法。 ８、複数の離間した受信器が使用され、該列は、第２のスケーリングされた列に 対する受信器位置が第１の列に対するのと同じになるよう受信器群間隔の倍数に 等しい間隔で前後して水中を牽引されることを特徴とする請求の範囲第３項記載 の方法。 ９、　αは１．１から３までの値を有することを特徴とする請求の範囲第１項記 載の方法。 １０、αは１．５から３までの値を有することを特徴とする請求の範囲第９項記 載の方法。 １１、βは１．１から３までの値を有することを特徴とする請求の範囲第１項記 載の方法。 １２、　βは１．５から３までの値を有することを特徴とする請求の範囲第１１ 項記載の方法。 １３、大地のインパルス応答の震動記録を作成することを特徴とする請求の範囲 第１項記載の方法。 １４、請求の範囲第１３項記載の方法によって得られた大地のインパルス応答の 震動記録。