JPH0820438B2 - Nondestructive measurement method of physical properties of formation using acoustic wave - Google Patents
Nondestructive measurement method of physical properties of formation using acoustic waveInfo
- Publication number
- JPH0820438B2 JPH0820438B2 JP21891992A JP21891992A JPH0820438B2 JP H0820438 B2 JPH0820438 B2 JP H0820438B2 JP 21891992 A JP21891992 A JP 21891992A JP 21891992 A JP21891992 A JP 21891992A JP H0820438 B2 JPH0820438 B2 JP H0820438B2
- Authority
- JP
- Japan
- Prior art keywords
- receiver
- pseudo
- source
- measurement method
- random code
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Description
【発明の詳細な説明】Detailed Description of the Invention
【0001】[0001]
【産業上の利用分野】本発明は、地層の物理特性の音響
波を用いた非破壊測定方法に係り、特に、土木工事や基
礎工事の計画、設計のための土質調査に用いるのに好適
な、海底堆積層における、間隙率、透水率、堆積物の変
化や、剪断変形率、剪断力等の断面分布状況を得るため
のものである。より詳細に述べると、本発明は、液体が
浸透した堆積層における、擬似ランダムバイナリシーク
エンスコードで駆動された音響伝達の使用に関するもの
である。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a nondestructive measurement method using acoustic waves of physical properties of a stratum, and is particularly suitable for use in soil investigation for planning and designing civil works and foundation works. The purpose is to obtain cross-sectional distributions such as porosity, water permeability, sediment changes, and shear deformation rate and shear force in the seabed sedimentary layer. More specifically, the present invention relates to the use of pseudo-random binary sequence code driven acoustic transfer in a liquid infiltrated sedimentary layer.
【0002】[0002]
【従来の技術】“Synthetic Seismograms for Mar
ine Sediments and Determinationof Porosity an
d Permerbility (海底堆積層の合成震動記録及び間隙
率と透水率の決定)”(Geophysics ,Volume 53,
No.8(1988年8月),1056−1057頁)と
称する発明者の論文で、発明者は、海底堆積層の垂直方
向の地震プロファイルの数値計算シミュレーションを発
表している。この論文中で引用した形の垂直方向の震動
プロファイルは、沖積鉱床の探査や、地質工学上の問題
解決や、水中構造物プロジェクトの建設に非常に有効で
ある。典型的な垂直方向の震動プロファイルにおいて
は、高出力の電気的震動源を用いて海底堆積層の研究が
行われていた。そのような研究において使用可能な周波
数範囲は、ほぼ100Hz から200Hz であることが
発見されており、これによって、海底堆積層中の透過距
離は、1000m 以上にも達することがあった。2. Description of the Related Art "Synthetic Seismograms for Mar"
ine Sediments and Determination of Porosity an
d Permerbility (synthetic seismic record of seafloor sediments and determination of porosity and permeability) "(Geophysics, Volume 53,
No. 8 (August 1988, pp. 1056-1057), the inventor has published a numerical simulation of vertical seismic profiles of submarine sediments. The vertical seismic profile of the form cited in this paper is very useful for exploration of alluvial deposits, geotechnical problem solving, and construction of underwater structure projects. In a typical vertical seismic profile, high power electrical seismic sources have been used to study submarine sediments. The frequency range that can be used in such studies has been found to be approximately 100 Hz to 200 Hz, which can lead to transmission distances in seabed sediments of up to 1000 m or more.
【0003】前記論文中で記載した一連の試験において
は、単一の探査孔を用いて、この同じ探査孔の中に地震
エネルギの発生源と受信器を共に配置して、垂直方向の
震動プロファイルを得るようしていた。前述の論文で詳
細に報告しているように、このような配置条件の元で
も、非常に興味ある結果が得られた。In the series of tests described in the above article, a single exploration hole was used in which the source of seismic energy and the receiver were both placed in the same exploration hole to obtain a vertical seismic profile. Was trying to get. As reported in detail in the above paper, very interesting results were obtained even under such arrangement conditions.
【0004】前記論文で記載したように、垂直方向の震
動プロファイルを用いることによって、理論的にも実際
的にも有効な結果を得ることができる。Biot の理論と
垂直方向の震動プロファイル方法によって得られた結果
を用いて、合成例についてスペクトル比を計算すれば、
海底堆積層の間隙率や透水率を決定することができる。
しかしながら、堆積層の物理特性をより正確に決定する
こと、特に、間隙率や透水率や堆積物の変化の正確な断
面分布状況を得ることが、待ち望まれていた。As described in the above paper, by using the vertical seismic profile, theoretically and practically valid results can be obtained. Using the Biot's theory and the results obtained by the vertical seismic profile method to calculate the spectral ratio for the synthetic example,
It is possible to determine the porosity and permeability of the seabed sedimentary layer.
However, it has been desired to more accurately determine the physical properties of the sedimentary layer, and in particular to obtain an accurate cross-sectional distribution of changes in porosity, water permeability and sediment.
【0005】“Measurements of Acoustic Wave
Velocities and Attenuation in Marine
Sediments(海底堆積層における音響波速度と減衰の測
定)”(Turgut and Yamamoto 、J.Acoust .
Soc. Am .´87(6)、1990年6月)にお
いて、堆積層の物理特性を測定するための非破壊方法が
論じられている(特願平2−333477号参照)。こ
の方法は、特に海底堆積層に適用され、間隙率や透水率
の値の断面分布を得るのに有用である。この方法は、
又、堆積物質の変化や、剪断変形率及び剪断力も明らか
にする。この方法において、速度の測定結果は、震動エ
ネルギ信号を、発振源から受信器に、直接探査孔間トモ
グラフィ(断層写真法)によって直接伝播させることに
よって得られる。この方法は、従来の探査孔間測定技術
を用いることなく、これ迄よりも、非常に優れた精度及
び明確さをもたらす。震動エネルギは、2つ以上の、互
いに離れて配置された垂直探査孔間に照射される。"Measurements of Acoustic Wave
Velocities and Attenuation in Marine
Sediments (Measurement of acoustic wave velocity and attenuation in seafloor sediments) "(Turgut and Yamamoto, J. Acoust.
Soc. Am. 87 (6), June 1990), a nondestructive method for measuring the physical properties of a deposited layer is discussed (see Japanese Patent Application No. 2-333477). This method is especially applied to seabed sedimentary layers and is useful for obtaining cross-sectional distributions of porosity and permeability values. This method
In addition, the change of the deposited material, the shear deformation rate and the shear force will be clarified. In this method, velocity measurements are obtained by propagating the seismic energy signal directly from the source to the receiver by direct interhole tomography. This method provides much greater accuracy and clarity than ever without the use of conventional inter-hole measurement techniques. Seismic energy is applied between two or more vertically spaced wells that are spaced apart from each other.
【0006】[0006]
【発明が解決しようとする課題】しかしながら、鉄鋳造
所のような雑音の多い環境下で、そのような技術を用い
る際には困難な点があり、1つの断面探査孔から他の断
面探査孔まで離すことのできる距離について、限度があ
った。更に、そのような測定技術を、高い周波数で用い
る際にも問題点があった。However, there is a difficulty in using such a technique in a noisy environment such as an iron foundry. There was a limit to the distance that I could separate. Furthermore, there are problems when using such measurement techniques at high frequencies.
【0007】即ち、従来の発振源は、爆薬、重垂落下、
板たたき等であり、規則波でなかったため、従来の音響
波では、周囲の雑音、騒音との識別が難しかった。特に
土中における発振音響波の減衰は著しく、10m 以上の
長距離間の測定は不可能であった。That is, the conventional oscillation sources are explosives, heavy drops,
Since it is a tapping, etc., and it is not a regular wave, it is difficult to distinguish it from ambient noise and noise with the conventional acoustic wave. In particular, the attenuation of oscillating acoustic waves in soil was remarkable, and it was impossible to measure over a long distance of 10 m or more.
【0008】従って、本発明は、特に、雑音の多い環境
下で実施された場合であっても、間隙率や透水率の値や
それらにの変化の、正確な断面分布を得て、正確で明確
な測定結果を得ることが可能な、新規で、非常に正確
な、地層の物理特性の音響波を用いた非破壊測定方法を
提供することを目的とする。Therefore, the present invention can obtain an accurate cross-sectional distribution of the values of porosity and water permeability and changes thereof, even when it is carried out in a noisy environment, and thus, it is possible to obtain accurate values. It is an object of the present invention to provide a new, highly accurate, non-destructive measurement method using acoustic waves of physical properties of the formation, which can obtain clear measurement results.
【0009】本発明の他の目的は、この方法を、1マイ
ル(1.6Km )以上の距離のように、非常に離れた断
面探査孔間で実施することにある。Another object of the invention is to implement the method between cross-section surveying holes that are very far apart, such as for distances of one mile (1.6 km) or more.
【0010】更に、本発明の他の目的は、これ迄は不可
能であった高い周波数での、新規で高精度の非破壊測定
方法を提供することにある。Still another object of the present invention is to provide a new and highly accurate nondestructive measuring method at a high frequency which has been impossible up to now.
【0011】本発明の他の目的は、高周波数の音響トモ
グラフィによって、間隙率や透水率を正確に評価するこ
とにある。Another object of the present invention is to accurately evaluate porosity and water permeability by high-frequency acoustic tomography.
【0012】本発明の更に他の目的は、地層の透水率分
布を、迅速、且つ、正確に決定するための新規な方法を
提供することにある。Still another object of the present invention is to provide a new method for determining the permeability distribution of a formation quickly and accurately.
【0013】本発明の他の重要な目的及び効果は、以下
の明細書や図面によって更に明らかになるであろう。Other important objects and effects of the present invention will become more apparent by the following specification and drawings.
【0014】[0014]
【課題を解決するための手段】本発明は、地層の物理特
性の音響波を用いた非破壊測定方法において、地層の表
面から少なくとも所定の深さまで伸び、入口を有する、
互いに所定距離だけ離れた少なくとも一対の探査孔を形
成し、前記探査孔の既知の位置に受信器を配置し、前記
探査孔の他の1つの中の、前記受信器に関して既知の位
置に、擬似ランダムコード震動エネルギの発振源を配置
し、前記擬似ランダムコードに従って該発振源を励起し
て、該発振源から前記受信器に震動エネルギを伝播さ
せ、前記発振源の位置から前記受信器の位置に向かって
伸びる経路中の前記震動エネルギの擬似ランダムコード
地震波特性を測定することにより、前記目的を達成した
ものである。The present invention provides a nondestructive measurement method using acoustic waves of physical properties of a formation, which extends from the surface of the formation to at least a predetermined depth and has an inlet.
Forming at least a pair of probe holes separated from each other by a predetermined distance, arranging a receiver at a known position of the probe hole, and simulating a position of the receiver in another one of the probe holes at a position known with respect to the receiver. Arranging an oscillation source of random code vibration energy, exciting the oscillation source according to the pseudo-random code to propagate vibration energy from the oscillation source to the receiver, and from the position of the oscillation source to the position of the receiver. The object has been achieved by measuring the pseudo-random code seismic wave characteristics of the vibration energy in a path extending toward.
【0015】更に、前記発振源を前記受信器に対して移
動し、該発振源を再び励起して、前記受信器に向かって
伸びる新しい経路中の擬似ランダムコード地震波特性を
更に測定するようにしたものである。Further, the oscillator source is moved with respect to the receiver and the oscillator source is re-excited to further measure the pseudo-random code seismic properties in a new path extending towards the receiver. It is a thing.
【0016】又、互いに所定距離だけ離れた少なくとも
一対の探査孔を形成し、少なくとも1つの前記探査孔内
の既知位置に、多数の受信器を互いに間隔を空けて配置
し、前記探査孔の他の1つの中の、前記受信器に関して
既知の位置に、擬似ランダムコード震動エネルギの発振
源を配置し、前記擬似ランダムコードに従って該発振源
を励起して、該発振源から前記受信器に震動エネルギを
伝播させ、前記発振源の位置から前記受信器の位置に向
かって伸びる多数の経路中の前記震動エネルギの擬似ラ
ンダムコード地震波特性を測定して記憶することによ
り、前記目的を達成したものである。Further, at least a pair of search holes separated from each other by a predetermined distance are formed, and a large number of receivers are arranged at known positions in at least one of the search holes. A source of quasi-random code seismic energy at a known location with respect to the receiver in one of, and exciting the oscillating source according to the pseudo-random code to oscillate energy from the source to the receiver. Is achieved by measuring and storing pseudo random code seismic wave characteristics of the vibration energy in a number of paths extending from the position of the oscillation source toward the position of the receiver. .
【0017】又、地層の表面から少なくとも所定の深さ
まで伸び、入口を有する、互いに所定距離だけ離れた少
なくとも一対の探査孔を形成し、前記探査孔間の地層表
面に実質的に沿う所定の既知の位置に、複数の受信器を
互いに間隔を空けて配置し、少なくとも1つの前記探査
孔内の既知位置に、多数の受信器を互いに間隔を空けて
配置し、前記探査孔の他の1つの中の、前記受信器に関
して既知の位置に、擬似ランダムコード震動エネルギの
発振源を配置し、前記擬似ランダムコードに従って該発
振源を励起して、該発振源から前記受信器に震動エネル
ギを伝播させ、前記発振源の位置から前記受信器の位置
に向かって伸びる多数の経路中の前記震動エネルギの擬
似ランダムコード地震波特性を測定し、前記発振源を前
記受信器に対して移動し、該発振源を再び励起して、前
記受信器に向かって伸びる新しい多数の経路中の擬似ラ
ンダムコード地震波特性を更に測定することにより、前
記目的を達成したものである。Further, at least a pair of exploration holes extending from the surface of the formation to at least a predetermined depth and having inlets and separated from each other by a predetermined distance are formed, and a predetermined known shape substantially along the formation surface between the exploration holes is formed. , A plurality of receivers are spaced apart from each other, and a plurality of receivers are spaced apart from each other at a known location within at least one of the exploration holes An oscillating source of pseudo-random code vibration energy is disposed at a known position in the inside of the device, and the oscillating source is excited according to the pseudo-random code to propagate the oscillating energy from the oscillating source to the receiver. , Measuring a pseudo-random code seismic wave characteristic of the seismic energy in a number of paths extending from the position of the oscillation source toward the position of the receiver, and setting the oscillation source to the receiver. Dynamic and, again exciting the emitting vibration source, by further measuring the pseudo-random code seismic properties in a new multiple paths extending toward the receiver, is obtained by achieving the above object.
【0018】又、前記震動エネルギの発振源を、圧電セ
ラミック変換器としたものである。Further, the oscillation source of the vibration energy is a piezoelectric ceramic converter.
【0019】又、前記探査孔を、互いに実質的に平行と
したものである。Further, the exploration holes are substantially parallel to each other.
【0020】又、前記受信器を、実質的に等間隔で配置
としたものである。Further, the receivers are arranged at substantially equal intervals.
【0021】又、前記地震波特性の測定に際して、伝播
した震動エネルギの固有減衰及び速度を測定するように
したものである。Further, in measuring the seismic wave characteristic, the characteristic attenuation and velocity of the propagating vibration energy are measured.
【0022】又、前記測定された地震波特性が、縦波速
度の周波数特性や固有減衰の周波数特性を含むようにし
たものである。Further, the measured seismic wave characteristics include a frequency characteristic of longitudinal wave velocity and a frequency characteristic of intrinsic damping.
【0023】又、前記擬似ランダム信号の周波数が変化
するようにしたものである。Further, the frequency of the pseudo random signal is changed.
【0024】又、前記擬似ランダムコードのキャリア周
波数が、ほぼ1、2、4、8及び10KHz の所で変化
するようにしたものである。Further, the carrier frequency of the pseudo random code is changed at about 1, 2, 4, 8, and 10 kHz.
【0025】又、前記各探査孔を、ほぼ70−200m
離したものである。Also, each of the above-mentioned exploration holes is approximately 70-200 m
It's separated.
【0026】[0026]
【作用】発明者の研究の結果、擬似ランダムバイナリシ
ークエンス(Pseude −Random Binary Sequence ;
以下PRBSと称する)の形の規則波の擬似雑音源を使
った電磁音響システムを使い、直接探査孔間トモグラフ
ィによって、震動エネルギ信号を発振源から受信器に直
接伝播させることにより、堆積層の特性を測定すること
が非常に重要であることがわかった。擬似ランダムバイ
ナリシークエンスは、例えば図1に示すように、所定の
不規則なパターンを有する繰返しシークエンスを持つ、
不規則な形状のマルチレベル信号である。この繰返しシ
ークエンスは、キャリア信号とも呼ばれる。As a result of the inventor's research, a pseudo random binary sequence (Pseude-Random Binary Sequence;
(Hereinafter referred to as PRBS), using an electromagnetic acoustic system using a regular-wave pseudo-noise source, direct propagation hole tomography propagates the seismic energy signal directly from the source to the receiver, It has been found that measuring the properties is very important. A pseudo-random binary sequence has a repeating sequence with a predetermined irregular pattern, for example, as shown in FIG.
It is an irregularly shaped multilevel signal. This repetitive sequence is also called a carrier signal.
【0027】地層の探査孔間配置を動作する際に、この
ような擬似ランダムバイナリシーエンスを使用すれば、
エンジンノイズや他の外来の干渉が存在する、雑音の多
い環境下であっても、効率の良い正確な動作が可能とな
り、実験者が、600フィート(960m )又は1マイ
ル以上のような非常に離れた距離を通して彼の断面探査
孔間信号を送っても、素晴らしい結果を得ることができ
るようになる。Using such a pseudo-random binary sequence when operating the inter-hole placement in the formation,
Efficient and accurate operation is possible even in noisy environments where engine noise and other extraneous interference are present, allowing the experimenter to be very sensitive, such as 600 feet (960 m) or more than a mile. Sending his cross-hole signal over large distances will also give excellent results.
【0028】即ち、図2に製鉄所における測定例を示す
如く、受信した音響波に、杭打ちや重機の走行等による
雑音や騒音と、PRBSコードによる信号が入り混ざっ
ている場合でも、PRBSコードとの相関をとることに
より、他の雑音と識別でき、発振から受信迄の時間測定
が可能となる。That is, as shown in a measurement example in a steel mill in FIG. 2, even if noises and noises due to pile driving, running of heavy machinery, etc. and PRBS code signals are mixed in the received acoustic waves, the PRBS code is mixed. By correlating with, it can be distinguished from other noises, and the time from oscillation to reception can be measured.
【0029】PRBSの使用は、非常に高い周波数であ
っても、素晴らしい動作を可能とする。PRBSエネル
ギの発振源は、これまでに使われていた圧電セラミック
変換器、又は、受信器としてのハイドロホン(Hydroph
one ;水中聴音器)を含む探査孔から離れた、もう1つ
の探査孔中の既知の場所に置かれたときに単一レベルの
規則波を発生する、スパーカ(Sparker)と呼ばれるも
のとは、大いに異なり、優れている。The use of PRBS allows excellent operation even at very high frequencies. The oscillation source of the PRBS energy is a piezoelectric ceramic converter that has been used up to now, or a hydrophone (Hydroph) as a receiver.
One, which produces a single-level regular wave when placed at a known location in another well, away from the well containing the hydrophone, is called a Sparker. Much different and better.
【0030】従って、本発明によれば、入口を有し、互
いに所定の距離だけ大きく離れて配置された、少なくと
も一対の探査孔を形成することによって、断面分布分析
結果を得るための地層の物理特性の音響波を用いた非破
壊測定方法を得ることができる。前記探査孔は、地層表
面から少なくとも既知の深さまで伸び、該探査孔内の所
定位置と、恐らく地層の表面に沿って互いに間隔を空け
て配置された、ハイドロホンやジオホン(Geohone;地
中聴音器)のような、複数の受信器を含む。擬似ランダ
ムバイナシリシークエンス発振器の形をとる雑音源は、
他の探査孔内の、前記受信器に関して既知の位置に配置
され、励起されて、該エネルギ源から全受信器に地震波
エネルギを伝播させる。発振源の位置から受信器の位置
に伸びる多数の経路について、発振源から受信器に直接
伝播された地震波の特性と経路を測定するに際して、重
要な測定は、地層の特性に関するものである。更に、発
振器探査孔内で発振器を別の既知の位置に動かして再び
発振器を励起することによって、受信器に向かって伸び
る多数の新しい経路についての測定を行うことができ
る。このような一連の測定の結果を組合せることによっ
て、堆積物の正確で且つ信頼性の高い断面の値や変化を
得ることが可能である。PRBS信号を直接伝播させ受
信することによって得られた結果は、従来利用可能なも
のに比べて、はるかに信頼性の高い結果をもたらす。Therefore, according to the present invention, by forming at least a pair of exploration holes having inlets and arranged at a large distance from each other, the physical properties of the formation for obtaining the cross-sectional distribution analysis result are formed. A nondestructive measurement method using a characteristic acoustic wave can be obtained. The exploration hole extends from the surface of the stratum to at least a known depth, and is located at a predetermined position in the exploration hole and possibly along the surface of the stratum at a distance from each other such as a hydrophone or geophone (geohone). Multiple receivers, such as a receiver). A noise source in the form of a pseudo-random binarity sequence oscillator is
Located in a known location with respect to the receiver in another probe hole and excited to propagate seismic energy from the energy source to all receivers. When measuring the characteristics and paths of seismic waves directly propagated from the oscillation source to the receiver for a number of paths extending from the position of the oscillation source to the position of the receiver, an important measurement relates to the characteristics of the stratum. Furthermore, by moving the oscillator to another known position in the oscillator probe hole and exciting the oscillator again, measurements can be made on a number of new paths extending towards the receiver. By combining the results of such a series of measurements, it is possible to obtain accurate and reliable cross-sectional values and changes of the deposit. The results obtained by directly propagating and receiving the PRBS signal give much more reliable results than previously available.
【0031】本発明によって、PRBS音響波の減衰を
用いて直接速度を測定することにより、長距離にわたっ
て信頼性の高い間隙率や透水率が得られ、より正確な透
水率の値を得ることができるということは、非常に重要
なことである。これは、石油の探鉱等の地下水関連問題
を解決するときに、非常に重要である。本発明による間
隙率の測定は、石油の探鉱の場合に、非序に有益であ
る。何故ならば、大部分の石油は堆積層の孔隙に吸収さ
れており、通常は、埋蔵されている石油のほんの少しの
割合しか発見できないからである。本発明の方法を利用
すれば、孔隙率が最も小さな場所、即ち、堆積層の孔隙
に吸収されていない自由な石油が存在する確率の高い場
所を正確に予測することができ、従って、透水率分布に
基づいて、石油の探査及び採取に好適な位置を示すこと
ができる。According to the present invention, by directly measuring the velocity using the attenuation of the PRBS acoustic wave, it is possible to obtain a highly reliable porosity and water permeability over a long distance, and to obtain a more accurate water permeability value. Being able to do is very important. This is very important when solving groundwater related problems such as oil exploration. The porosity measurement according to the invention is intriguingly useful in the case of petroleum exploration. This is because most of the oil is absorbed in the pores of the sedimentary layer and normally only a small percentage of the oil reserves can be found. The method of the present invention can be used to accurately predict where the porosity is the smallest, i.e. where there is a high probability that there is unabsorbed free oil in the porosity of the sedimentary layer, and therefore the permeability. Based on the distribution, suitable locations for oil exploration and extraction can be indicated.
【0032】本発明においては、エネルギ源と受信器を
同じ探査孔内に配置した単一の探査孔のみを有する垂直
方向の震動プロファイルを利用する代わりに、既に説明
したように、互いに離れた探査孔を利用して、1つの探
査孔から他の探査孔への伝播を測定することが重要であ
る。ここでは、この新しいPRBS技術を、縮めて、
「PRBS探査孔間トモグラフィ」と称する。この方法
は、地面の透水率プロファイルをチェックする際に非常
に有効であることが判った。In the present invention, instead of utilizing a vertical seismic profile having only a single exploration hole with the energy source and receiver located in the same exploration hole, as already explained, It is important to utilize holes to measure the propagation from one exploration hole to another. Here, this new PRBS technology is shortened,
It is called "PRBS inter-hole tomography". This method has been found to be very effective in checking the permeability profile of the ground.
【0033】本発明の一実施例においては、発振源とし
て、圧電音響素子を用いて、PRBSキャリア周波数を
1、2、4、8及び10KHz で変化させることによ
り、圧電発振源からの音響伝播を利用し、受信器とし
て、堆積層の表面に沿って、又、前記圧電発振源を含む
探査孔から離れた場所の独立した探査孔内を下に向かっ
て伸びるように配置したハイドロホンを用いて実験を行
うことができる。探査孔内の発振源の場所を変えること
によって、同じ方法で追加の実験が行われる。In one embodiment of the present invention, a piezoelectric acoustic element is used as an oscillation source, and the PRBS carrier frequency is changed by 1, 2, 4, 8 and 10 KHz, so that the acoustic propagation from the piezoelectric oscillation source is caused. Utilizing, as a receiver, a hydrophone arranged to extend downward along the surface of the deposited layer and in an independent exploration hole apart from the exploration hole including the piezoelectric oscillation source. You can carry out experiments. Additional experiments are performed in the same way by changing the location of the source within the borehole.
【0034】一例として、前記離隔した探査孔は、堆積
物中にほぼ鉛直方向に伸びることができる。これらの探
査孔は、互いに、例えば約70−200m だけ離して配
置することができ、発振源から全てのハイドロホンに対
して、直接的な伝播が行われるようにして、探査孔間の
縦波(圧縮波)速度及び固有減衰を測定するように、各
ハイドロホンからの信号を、アンプ及び/又はレコーダ
を介して測定することができる。固有減衰の測定は、ス
ペクトル比による方法を用いて行うことができ、縦波速
度及び固有減衰は、位相スペクトル及び伝達関数から容
易に計算することができる。As an example, the spaced apart exploration holes may extend substantially vertically into the deposit. These exploration holes can be arranged at a distance of, for example, about 70-200 m from each other, and direct propagation is performed from the oscillation source to all hydrophones, so that longitudinal waves between the exploration holes are provided. The signal from each hydrophone can be measured via an amplifier and / or a recorder, so as to measure the (compressed wave) velocity and the intrinsic damping. Intrinsic attenuation measurements can be made using the spectral ratio method, and longitudinal wave velocity and intrinsic attenuation can be easily calculated from the phase spectrum and transfer function.
【0035】即ち、Biot の理論より、縦波速度V及び
固有減衰Q-1は、次式で表わされる。That is, according to the theory of Biot, the longitudinal wave velocity V and the intrinsic damping Q -1 are expressed by the following equations.
【0036】 V=V0 +[(R2 −R1 )f ×360°]/θ …(1) Q-1=[V/{πf (R2 −R1 )}]・{ln (R1 /R2 ) +ln (A1 /A2 )} …(2)V = V0 + [(R2-R1) f × 360 °] / θ (1) Q- 1 = [V / {πf (R2-R1)}] {ln (R1 / R2) + ln ( A1 / A2)} (2)
【0037】ここで、V0 は基準周波数ω0 における縦
波速度、R1 、R2 は、発振源から2つの受信器までの
距離、f は発振周波数、θは位相角、A1 /A2 は、2
つの信号の振幅比である。Here, V0 is the longitudinal wave velocity at the reference frequency ω0, R1 and R2 are the distances from the oscillation source to the two receivers, f is the oscillation frequency, θ is the phase angle, and A1 / A2 is 2
It is the amplitude ratio of two signals.
【0038】このようにして、本発明による探査孔間ト
モグラフィは、空間的な異質部分をチェックするのに使
用される。In this way, inter-hole tomography according to the present invention is used to check for spatial disparities.
【0039】従って、本発明においては、堆積物質中で
の音響パルス伝播が行われ、特に固有減衰が最大である
特定の周波数範囲を含んで伝播性のある縦波速度を示す
ことが判った。間隙率及び透水率は、前記Biot の理論
を用いることによって、縦波速度及び固有減衰から容易
に評価される。なお、その詳細は、全て、我々の“Mea
surements of acoustic wave velocities and atten
uation in marine sediments(海底堆積層中の音響波速
度及び減衰の測定)”J.Acoust .Soc.Am .,V
ol.87,No.6(1990年6月),2376−23
83頁に詳細に記載されているので、詳しい説明は省略
する。Therefore, it has been found that in the present invention, acoustic pulse propagation is carried out in the deposited material, and in particular, it exhibits a propagating longitudinal wave velocity including a specific frequency range where the intrinsic attenuation is maximum. Porosity and permeability are easily estimated from longitudinal wave velocity and intrinsic damping by using the Biot theory. The details are all in our "Mea
surements of acoustic wave velocities and atten
uation in marine sediments (Measurement of acoustic wave velocity and attenuation in seafloor sediments) "J. Acoust. Soc. Am., V
ol. 87, No. 6 (June 1990), 2376-23.
The detailed description is omitted because it is described in detail on page 83.
【0040】前記探査孔間トモグラフィに、PRBS発
生器を備えた圧電音響変換器を用いることによって、縦
波だけでなく横波も発生する。従って、横波速度の知識
から、堆積形状内の横波速度分布も正確に測定でき、従
って、堆積物の剪断変形率及び剪断力も決定することが
できる。これらは、構造物を建設する目的や、石油やガ
スのような炭化水素を正確に決定したり、存在を知る場
合に重要である。By using a piezoelectric acoustic transducer equipped with a PRBS generator for the inter-probe tomography, not only longitudinal waves but also transverse waves are generated. Therefore, from the knowledge of the shear wave velocity, the shear wave velocity distribution within the deposit shape can also be accurately measured, and thus the shear deformation rate and shear force of the deposit can also be determined. These are important for the purpose of constructing structures and for accurately determining and knowing the presence of hydrocarbons such as oil and gas.
【0041】横波速度の値は、直接、堆積物の剪断変形
率や剪断強度の値を決定する。縦波速度の値に加えて横
波速度の値を知ることにより、堆積形状内における、石
油やガスのような炭化水素の存在や非存在に関する重要
な情報を得ることができる。The shear wave velocity value directly determines the shear deformation rate and the shear strength value of the deposit. Knowing the value of the transverse wave velocity in addition to the value of the longitudinal wave velocity can provide important information about the presence or absence of hydrocarbons such as oil and gas within the deposit geometry.
【0042】[0042]
【実施例】以下図面を参照して、本発明の実施例を詳細
に説明する。Embodiments of the present invention will now be described in detail with reference to the drawings.
【0043】図3は、音響トモグラフィ測定のエネルギ
伝播源として使われる、探査孔11内に配置されたPR
BS発振源10を示すものである。探査孔11は、例え
ば直径5インチ(12.7cm)のプラスチックケーシン
グで作ることができる。符号12は、堆積層の表面に沿
って間隔を空けて配置された多数のハイドロホンを示
し、符号13は、沿直方向に伸びる探査孔14内で互い
に沿直方向に間隔を空けて配置された多数のハイドロホ
ンを示している。探査孔14と11は、図示の如く、互
いに、例えば70−200m 又はそれ以上の距離だけ離
れて配置されており、各探査孔自体も、例えば3−10
m の深さとされ、互いに平行とされている。従って、測
定は、好ましくは、発振源によって伝播される波動エネ
ルギに関して、異なる角度の広い範囲に亘って同時に行
われる。そのような測定は全て、堆積層内のエコー又は
他の現象を必要とすることなく、PRBS発振源10か
ら受信器であるハイドロホン12、13への直接伝播を
基礎として行われる。FIG. 3 shows a PR arranged in the probe hole 11 used as an energy propagation source for acoustic tomography measurement.
1 shows a BS oscillation source 10. The exploration hole 11 can be made of, for example, a plastic casing having a diameter of 5 inches (12.7 cm). Reference numeral 12 denotes a number of hydrophones arranged at intervals along the surface of the sedimentary layer, and reference numeral 13 is arranged at intervals in the longitudinal direction in the exploration holes 14 extending in the longitudinal direction. It shows a number of hydrophones. The exploration holes 14 and 11 are, as shown, separated from each other by a distance of, for example, 70-200 m or more.
It has a depth of m and is parallel to each other. Therefore, the measurements are preferably made simultaneously over a wide range of different angles with respect to the wave energy propagated by the oscillator source. All such measurements are based on direct propagation from the PRBS source 10 to the receiver hydrophones 12, 13 without the need for echoes or other phenomena in the deposited layer.
【0044】発振源は、PRBSコード発生器21とパ
ワーアンプ20によってパワーを与えられ、受信信号
は、アンプ23で増幅されて、レコーダ22で記録され
る。The oscillation source is powered by the PRBS code generator 21 and the power amplifier 20, and the received signal is amplified by the amplifier 23 and recorded by the recorder 22.
【0045】図の位置10′に発振源10を置いて多数
の測定を行った後、発振源は、他の位置10′′に下げ
られ、測定プロセスが繰返される。ハイドロホン12、
13のいずれも移動させることなく、発振源探査孔内で
の他の様々な位置で更に他の測定が行われる。After placing the oscillator 10 at the position 10 'in the figure and making a number of measurements, the oscillator is lowered to another position 10''and the measurement process is repeated. Hydrophone 12,
Still other measurements are made at various other locations within the source probe hole without moving any of the 13.
【0046】本発明の1実施例においては、PRBS信
号を発生し測定するために用いられる主な装置は、発振
源、ハイドロホンアレイ、信号アンプ及びフィルタ、P
RBS発生器、周波数シンセサイザ、パワーアンプ、関
数発生器、データ獲得ボックス、A/Bスイッチボック
ス、直流電源及びコンピュータである。この装置は、図
4に示され、以下、更に詳細に説明する。In one embodiment of the present invention, the main equipment used to generate and measure the PRBS signal is the oscillator source, hydrophone array, signal amplifier and filter, P
RBS generator, frequency synthesizer, power amplifier, function generator, data acquisition box, A / B switch box, DC power supply and computer. This device is shown in FIG. 4 and will be described in more detail below.
【0047】一般的に30で示される周波数シンセサイ
ザ(Frequency Synthesizer)は、変調可能な、一定
振幅、周波数及び位相のキャリア波を発生する。変調
は、ここで使用されたように、PRBS発生器によって
引き起され、キャリア波の周波数、位相及び/又は振幅
の制御された変化を生じる。これは、信号波を伝播する
ために、任意の周波数を取ることができる。周波数シン
セサイザ30は、ライン54から電力を受取る。A frequency synthesizer, generally designated 30, generates a carrier wave of constant amplitude, frequency and phase that can be modulated. Modulation, as used herein, is caused by a PRBS generator, resulting in a controlled change in frequency, phase and / or amplitude of the carrier wave. It can take any frequency to propagate the signal wave. Frequency synthesizer 30 receives power from line 54.
【0048】周波数シンセサイザは、ライン57を通し
て、一般的に32で示されている擬似ランダムバイナリ
シークエンス(PRBS)コード発生器(PRBS C
odeGenerator)に接続されている。PRBSコード発
生器32は、ライン56を通して電力を受取る。このP
RBSコード発生器32は、例えば、212又は4096
ビットまでの長さの擬似ランダムバイナリシークエンス
を発生する。この発生器からの出力信号は、振幅、位相
又はパルス変調器によって変調され得る。PRBSコー
ド発生器32の出力コードの周波数は、例えば周波数シ
ンセサイザ30によって発生される周波数の1/64と
されている。The frequency synthesizer is connected through line 57 to a pseudo random binary sequence (PRBS) code generator (PRBS C), generally indicated at 32.
odeGenerator). The PRBS code generator 32 receives power over line 56. This P
The RBS code generator 32 may be, for example, 2 12 or 4096.
Generate a pseudo-random binary sequence up to bit in length. The output signal from this generator can be modulated by an amplitude, phase or pulse modulator. The frequency of the output code of the PRBS code generator 32 is, for example, 1/64 of the frequency generated by the frequency synthesizer 30.
【0049】PRBSコード発生器32は、ライン58
を通して、一般的に34で示されている可変出力スイッ
チングパワーアンプ(Power Amplifier)に接続され
ている。パワーアンプ34は、任意の最大電力出力及び
最大電圧出力を持つことができる。都合のいい例では、
それぞれ4k VAと1400ボルトである。出力で最大
電圧を得るために、パワーアンプ34に対する入力は、
ライン66及び68を介して、例えば+−1.0ボルト
に制限されるべきである。パワーアンプ34は、ライン
64を通してグラウンドに接続されている。パワーアン
プ34は、全ての種類の短絡や過大入力に耐える。パワ
ーアンプ34の周波数応答は、好ましくは700Hz か
ら20k Hz の間で実質的に平坦とされる。The PRBS code generator 32 has a line 58.
Through a variable output switching power amplifier (Power Amplifier), generally designated 34. The power amplifier 34 can have any maximum power output and maximum voltage output. In a convenient example,
They are 4k VA and 1400 volts respectively. To get the maximum voltage at the output, the input to the power amplifier 34 is
It should be limited to + -1.0 volts, for example, via lines 66 and 68. The power amplifier 34 is connected to the ground through the line 64. The power amplifier 34 withstands all types of short circuits and excessive inputs. The frequency response of power amplifier 34 is preferably substantially flat between 700 Hz and 20 kHz.
【0050】一般的に36で示される圧電音響変換器
(Source )は、ライン70を通して、一般的に34で
示されるパワーアンプに接続される。圧電音響変換器3
6は、電気信号である電圧を発生する。圧電音響変換器
は、音響波を発生する。圧電音響変換器36は、例えば
ビルトインされた4倍の電圧トランスを有している。従
って、圧電音響変換器36に送られる電圧は、好ましく
は該圧電音響変換器36の最大電圧定格の約1/4に限
定されるべきである。A piezoelectric acoustic transducer (Source), shown generally at 36, is connected through line 70 to a power amplifier, shown generally at 34. Piezoelectric transducer 3
6 generates a voltage which is an electric signal. The piezoelectric acoustic transducer produces an acoustic wave. The piezoelectric acoustic transducer 36 has, for example, a built-in quadruple voltage transformer. Therefore, the voltage delivered to the piezo-acoustic transducer 36 should preferably be limited to about 1/4 of the maximum voltage rating of the piezo-acoustic transducer 36.
【0051】一般的に38で示されるハイドロホンアレ
イ(Hydrophone Array)を作り出すために、例えば2
4個のハイドロホンが、所定の固定されたパターンで配
列されている。ハイドロホンアレイ38は、例えば、油
によって封入された24素子のアレイである。このハイ
ドロホンは、水に媒介される音波に応答する電気音響変
換器であり、本質的に等価な電気波を発生する。各ハイ
ドロホンは、周波数が1.0Hz 又はそれ以下の低い信
号を感知することができる。一般的に40で示される1
2ボルト直流電源(12 Volt DC Supply )は、
ライン72を通してハイドロホンアレイ38に電力を与
える。To produce a Hydrophone Array, generally designated 38, for example, 2
Four hydrophones are arranged in a fixed pattern. The hydrophone array 38 is, for example, a 24-element array filled with oil. The hydrophone is an electroacoustic transducer that responds to water-borne sound waves and produces essentially equivalent electrical waves. Each hydrophone is capable of sensing low signals with a frequency of 1.0 Hz or less. 1 generally designated by 40
2 volt DC power supply (12 Volt DC Supply)
The hydrophone array 38 is powered through line 72.
【0052】一般的に74と76で示されるケーブル
は、ハイドロホンアレイ38を、一般的に42で示され
るA/Bスイッチボックス(A/B Switch Box)に
接続する。例えば12個のハイドロホンが、端子Aに接
続され、残りの12個のハイドロホンが、端子Bに接続
される。Cables, generally designated 74 and 76, connect the hydrophone array 38 to an A / B switch box, generally designated 42. For example, 12 hydrophones are connected to the terminal A, and the remaining 12 hydrophones are connected to the terminal B.
【0053】一般的に84で示される、好ましくはRS
232ケーブル等であるケーブルが、A/Bスイッチボ
ックスを、一般的に44で示されるカスタム信号増幅フ
ィルタ(Signal Amplifier Filter )に接続する。
信号増幅フィルタ44は、ライン90を通してグラウン
ドに接続されている。信号増幅フィルタ44は、ライン
86を通して15ボルト、ライン88を通して−15ボ
ルトのような低い電圧に接続される。信号増幅フィルタ
44は、好ましくは、可変ゲイン、差動入力、回路の一
方がアースされた出力信号の増幅器等である。入力差動
信号は、好ましくは、500Hz の近傍で、30 dB高
域通過される。この信号は、それから、好ましくは、例
えば、1、2、5、10、20、50、100、20
0、500、1000、2000のゲインのいずれかに
よって増幅される。この特別のアンプ44は、例えば1
2チャンネルを、上記ゲインによって、同時に増幅する
ことができる。しかしながら、12チャンネル以上又は
以下の増幅が可能な他のアンプも使用できる。一般的に
52で示される手動スイッチが、信号増幅フィルタ44
に接続されている。Numeral 84, preferably RS
A cable, such as a 232 cable, connects the A / B switch box to a custom signal amplification filter, indicated generally at 44 (Signal Amplifier Filter).
The signal amplification filter 44 is connected to the ground through the line 90. The signal amplification filter 44 is connected to a low voltage such as 15 volts through line 86 and -15 volts through line 88. The signal amplification filter 44 is preferably a variable gain, a differential input, an output signal amplifier in which one of the circuits is grounded, or the like. The input differential signal is preferably 30 dB high pass in the vicinity of 500 Hz. This signal is then preferably, for example, 1, 2, 5, 10, 20, 50, 100, 20.
It is amplified by any of the gains of 0, 500, 1000 and 2000. This special amplifier 44 is, for example,
Two channels can be simultaneously amplified by the above gain. However, other amplifiers capable of amplifying more than 12 channels or less can also be used. A manual switch, generally indicated at 52, is a signal amplification filter 44.
It is connected to the.
【0054】一般的に92で示される、好ましくはRS
232Cケーブルであるケーブルが、信号増幅フィルタ
44を、一般的に46で示されるデータ獲得ボックス
(Data Acquisition Box)に接続する。PRBSコ
ード発生器32は、又、ライン60を通してデータ獲得
ボックス46に接続されている。更に、PRBSコード
発生器32は、ライン62を通して、一般的に50で示
されるファンクションゼネレータ(Function Generat
or)に接続され、これは次に、ライン63を通してデー
タ獲得ボックス46に接続されている。データ獲得ボッ
クス46は、好ましくは接続ボックスとして作用し、P
RBS発生器32及び信号増幅フィルタ44を、一般的
に96で示される、好ましくはフラットケーブルである
ケーブルへ接続し、このケーブル96は次に、386/
33MHz IBMコンパチブルコンピュータのような、
膨大な記憶容量を備えた、一般的に48で示されるコン
ピュータ(Computer )中に内臓されたマルチチャンネ
ルアナログ/デジタル(A/D)ボードに接続される。
このA/Dボードは、16チャンネルを有することがで
きるが、それ以上又は以下であってもよい。このA/D
ボードは、好ましくは、ハイドロホンからの感知信号を
チャンネル1−12を通して受取り、PRBS発生器の
出力信号をチャンネル13を通して受取り、外部クロッ
クをチャンネル15を通して受取る。Numeral 92, preferably RS
A cable, which is a 232C cable, connects the signal amplification filter 44 to a data acquisition box, generally indicated at 46. The PRBS code generator 32 is also connected to the data acquisition box 46 via line 60. In addition, the PRBS code generator 32 is connected through a line 62 to a function generator, generally designated 50.
or)), which in turn is connected to the data acquisition box 46 via line 63. The data acquisition box 46 preferably acts as a connection box and P
The RBS generator 32 and the signal amplification filter 44 are connected to a cable, generally designated 96, which is preferably a flat cable, which cable 96 then
Like a 33MHz IBM compatible computer,
It is connected to a multi-channel analog / digital (A / D) board built into a computer (Computer), generally designated 48, with enormous storage capacity.
The A / D board can have 16 channels, but can have more or less. This A / D
The board preferably receives the sense signal from the hydrophone on channels 1-12, the output signal of the PRBS generator on channel 13 and the external clock on channel 15.
【0055】周波数シンセサイザ30は、PRBSコー
ド発生器32中のデジタル動作に必要なキャリア周波数
を供給する。このPRBSコード発生器は、キャリア波
を変調して、擬似ランダムバイナリシークエンス信号波
を伝達する。パワーアンプ34は、この信号波を増幅す
る。圧電変換器36は、この増幅された信号波を伝達す
る。The frequency synthesizer 30 supplies the carrier frequency required for digital operation in the PRBS code generator 32. The PRBS code generator modulates a carrier wave and transmits a pseudo random binary sequence signal wave. The power amplifier 34 amplifies this signal wave. The piezoelectric converter 36 transmits the amplified signal wave.
【0056】ハイドロホンアレイ38は、擬似ランダム
バイナリシークエンス信号波を感知する。ハイドロホン
アレイ38は、24個の各ハイドロホンを含むことがで
きる。12個のハイドロホンによって感知された信号波
は、好ましくは、A/Bスイッチボックス42の端子A
に伝えられる。他の12個のハイドロホンによって感知
された信号波は、好ましくは、A/Bスイッチボックス
42の端子Bに伝えられる。信号増幅フィルタ44に取
付けられたスイッチ52は、信号増幅フィルタ44が、
該スイッチ52の位置に依存して、端子A又は端子Bか
ら、感知された信号波を受取るようにする。信号増幅フ
ィルタ44は、各ハイドロホンによって感知された信号
波を増幅してフィルタを通し、この情報をデータ獲得ボ
ックス46に出力する。データ獲得ボックス46は、外
部クロック信号及びPRBS発生器32からの元の擬似
ランダムバイナリシークエンス信号も受取っている。デ
ータ獲得ボックス46は、更に、ファンクションゼネレ
ータ50からの入力も受取っている。データ獲得ボック
ス46は、この情報を、386/33MHz IBMコン
パチブルコンピュータのようなコンピュータ48中に内
蔵されたマルチチャンネルアナログ/デジタル(A/
D)ボードに伝える。The hydrophone array 38 senses a pseudo-random binary sequence signal wave. The hydrophone array 38 can include 24 hydrophones each. The signal wave sensed by the twelve hydrophones is preferably the terminal A of the A / B switch box 42.
Conveyed to. Signal waves sensed by the other twelve hydrophones are preferably transmitted to the terminal B of the A / B switch box 42. The switch 52 attached to the signal amplification filter 44 is
Depending on the position of the switch 52, the sensed signal wave is received from the terminal A or the terminal B. The signal amplification filter 44 amplifies and filters the signal wave sensed by each hydrophone and outputs this information to the data acquisition box 46. The data acquisition box 46 also receives the external clock signal and the original pseudo-random binary sequence signal from the PRBS generator 32. The data acquisition box 46 also receives input from the function generator 50. The data acquisition box 46 sends this information to a multi-channel analog / digital (A /
D) Tell the board.
【0057】この装置を、最も好適に動作させるため
に、PRBSコード発生器32によって発生される周波
数レートは、好ましくは、発振源周波数の4倍の範囲で
ある。A/Dボードのコンピュータサンプリングレート
が高くなるほど、同時に測定可能なハイドロホンの数が
少なくなる。何故ならば、A/Dボードでの入力レート
は、チャンネルとサンプリングレートを掛けた数に等し
く走査しなければならないからである。好ましくは、1
つのハイドロホンにつき1チャンネルある。これは、も
し4個のハイドロホンのみが一時にサンプリング可能で
あれば、コンピュータは好ましくは、始めの4個のハイ
ドロホン、それから、2番目の4個のハイドロホン、そ
の後、3番目の4個のハイドロホンからの感知データを
サンプリングし、又は記録する。コンピュータは、この
サンプリングを常に繰返すことができる。For the best operation of this device, the frequency rate generated by the PRBS code generator 32 is preferably in the range of four times the oscillator source frequency. The higher the computer sampling rate of the A / D board, the fewer hydrophones can be measured simultaneously. This is because the input rate on the A / D board must scan equal to the number of channels times the sampling rate. Preferably 1
There is one channel for each hydrophone. This means that if only 4 hydrophones can be sampled at a time, the computer preferably has the first 4 hydrophones, then the second 4 hydrophones, and then the third 4 hydrophones. Sampling or recording the sensing data from the hydrophone. The computer can always repeat this sampling.
【0058】平均化は、PRBS信号のような多重シー
クエンス長を取入れて、互いの上に積重ねていき、雑音
を軽減して信号を強調する処理である。平均化がいつも
必要でない場合には、コンピュータの内部クロックは、
感知データを記録するのに使用される。Averaging is a process of taking in multiple sequence lengths such as PRBS signals and stacking them on top of each other to reduce noise and emphasize signals. If averaging is not always necessary, the computer's internal clock is
Used to record sensing data.
【0059】コンピュータ48、ファンクションゼネレ
ータ50及びPRBSコード発生器32からの外部クロ
ックは、共同して平均化を実行するように作用する。P
RBSコード発生器からの外部クロックは、感知信号を
PRBS信号と同期させる。平均化に際して、PRBS
発生器からの外部クロックは、発振源信号とハイドロホ
ンによって測定された信号の間でドリフトが発生しない
ように使われるべきである。ファンクションゼネレータ
50は、いつ擬似ランダムバイナリシークエンスが終了
したかをコンピュータに伝える。バイナリシークエンス
の指定された数が感知されると、これらのシークエンス
は、違いの上に積重ねられて雑音を軽減し、信号を強調
する。発振源周波数の4倍であるPRBS発生器の周波
数レートは、入力周波数に適合するように、ファンクシ
ョンゼネレータを使って分周される。例えば、1k Hz
では、外部クロックの周波数は4k Hz であり、8k H
zのサンプリングレートで1チャンネルを測定するため
に、この信号は、ファンクションゼネレータによって1
4サイクルに分周される。従って、A/Dボードの外部
クロックチャンネルに適用される入力レートは56k H
z になる。A/Dボードのチャンネル15からくる方形
パルスは、外部クロックが使用中であることを示す。External clocks from the computer 48, the function generator 50 and the PRBS code generator 32 work together to perform the averaging. P
An external clock from the RBS code generator synchronizes the sense signal with the PRBS signal. PRBS for averaging
The external clock from the generator should be used so that there is no drift between the source signal and the signal measured by the hydrophone. The function generator 50 tells the computer when the pseudo-random binary sequence is complete. When a specified number of binary sequences are sensed, these sequences are stacked on top of the differences to reduce noise and enhance the signal. The frequency rate of the PRBS generator, which is four times the source frequency, is divided using the function generator to match the input frequency. For example, 1 kHz
Then, the frequency of the external clock is 4 kHz, 8 kHz
To measure one channel at a sampling rate of z, this signal is
Divided into 4 cycles. Therefore, the input rate applied to the external clock channel of the A / D board is 56 kHz.
becomes z. A square pulse coming from channel 15 of the A / D board indicates that the external clock is in use.
【0060】コンピュータのメモリに蓄えられたデータ
獲得プログラムは、ハイドロホンからの感知信号を読み
取り、感知信号をメモリに蓄え、平均化のような計算を
実行するために、外部又は内部クロックを使用する。The data acquisition program stored in the memory of the computer uses an external or internal clock to read the sensing signal from the hydrophone, store the sensing signal in the memory, and perform calculations such as averaging. .
【0061】ファンクションゼネレータ50、PRBS
コード発生器32及び適当なデータ獲得プログラムを備
えたコンピュータ48が、発振源信号と測定信号の相関
を取り、平均化処理を実行するために、共同で作用す
る。この相関の出力は、1つの特定発振源位置で各ハイ
ドロホンの1つに対して一連の時間で見ることができ
る。これから、横波速度、縦波速度、剪断力、間隙率、
透水性及び標準貫通力が決定され得る。標準貫通力は、
貫通させるためにパイプを打撃しなければならない回数
によって定義される。Function generator 50, PRBS
A computer 48 equipped with a code generator 32 and a suitable data acquisition program work together to correlate the oscillating source signal with the measured signal and to perform the averaging process. The output of this correlation can be seen in a series of times for one of each hydrophone at one particular source position. From this, transverse wave velocity, longitudinal wave velocity, shear force, porosity,
Water permeability and standard penetration can be determined. The standard penetration is
It is defined by the number of times the pipe must be struck to penetrate.
【0062】本発明の一実施例に従い、60m の深さの
2本の深い探査孔を使って、探査孔間測定を行った。圧
電音響変換器36は、一方の探査孔中に配置され、ハイ
ドロホンアレイ38は、他方の探査孔中に配置された。
探査孔間の距離は66.21m であった。シークエンス
長さ4095を有するPRBSコードが、170 dBの
発振源を、3k Hz のキャリア周波数で励起するために
使用された。4回の測定の平均が取られた。コンピュー
タのサンプリングレートは12k Hz であった。最初に
感知されたデータは、図5及び図6の時間経過グラフに
示す如くであった。これらのグラフは、信号が最初にハ
イドロホンアレイ38に着いた時の時刻を示している。
伝播時間倒置法のコンピュータプログラムを使って、縦
波速度のトモグラフを作成することができた。一般的な
合同ストレス割合を仮定して、縦波速度プロファイルか
ら、間隙率、横波速度、剪断力及びSPT(標準貫通
力)打撃回数を決定し、それからトモグラフを作ること
ができる。最後に、透水性及び粒度が、2つの異なるP
RBSキャリア周波数で測定した縦波速度トモグラフか
ら評価できた。In accordance with one embodiment of the present invention, inter-probe measurements were made using two deep exploration holes 60m deep. The piezo-acoustic transducer 36 was placed in one probe hole and the hydrophone array 38 was placed in the other probe hole.
The distance between the exploration holes was 66.21 m. A PRBS code with a sequence length of 4095 was used to excite a 170 dB source with a carrier frequency of 3 kHz. The average of 4 measurements was taken. The computer sampling rate was 12 kHz. The first sensed data was as shown in the time course graphs of FIGS. These graphs show the time when the signal first arrived at the hydrophone array 38.
Using a computer program of the time-of-flight inversion method, we were able to create a longitudinal wave velocity tomograph. Assuming a general joint stress rate, the porosity, shear wave velocity, shear force and SPT (standard penetration force) hit count can be determined from the longitudinal wave velocity profile and a tomograph can be made from it. Finally, the water permeability and particle size are two different P
It could be evaluated from the longitudinal wave velocity tomography measured at the RBS carrier frequency.
【0063】堆積物データの処理及び解釈は、定性的及
び定量的性質を共に有する。堆積物データ処理の品質
は、処理データをいかに上手く解釈するかによって測定
される。図5及び図6に示した時間グラフは、非常に優
れた性質の結果を受取っている。 この例においては、
図4に示される本発明の実施例を用い、コンピュータ4
8は、12個のハイドロホンからのデータを受取ること
が可能な、A/Dボードを含んでいた。もしA/Dボー
ドが、24チャンネルを受信可能であったならば、A/
Bスイッチボックス42及びスイッチ52は不要だった
であろう。The processing and interpretation of sediment data has both qualitative and quantitative properties. The quality of sediment data processing is measured by how well the processing data is interpreted. The time graphs shown in Figures 5 and 6 receive results of very good nature. In this example,
Using the embodiment of the invention shown in FIG.
The 8 included an A / D board capable of receiving data from 12 hydrophones. If the A / D board was able to receive 24 channels, A / D
B switch box 42 and switch 52 would not have been needed.
【0064】この例においては、図4に示される実施例
を用いて、170 dBの発振源が使用されていた。もし
発振源の dBが大きければ、1.0マイル離れた探査孔
間の品質の高い測定結果を得ることができるであろう。In this example, using the embodiment shown in FIG. 4, a 170 dB oscillation source was used. If the source dB is large, it will be possible to obtain high quality measurements between the probe holes 1.0 miles apart.
【0065】更に、実験に際して、4回の平均が行われ
ていた。しかしながら、本発明を利用する際には、任意
の回数の平均を行うことができる。Further, in the experiment, four times of averaging were performed. However, any number of averagings can be performed when utilizing the present invention.
【0066】本発明は、ある特定の実施例を引用して説
明されていたが、本発明の精神及び範囲から離れること
なく、様々な変形を行うことができることは明らかであ
る。例えば、本発明による探査孔間トモグラフィのため
に、より多い探査孔を用いることができる。又、それぞ
れの位置が判っていれば、PRBS発振源と受信器間の
位置の多くの相対的な位置関係や変形を用いることがで
きる。更に、様々な理論や関係式を用いることによっ
て、空隙率、透水性、透水性プロファイル、剪断変形率
や剪断力のような特性の決定及び計算が行われていた
が、PRBS発振源がある探査孔内に位置され、該PR
BS発振源から離れた別の探査孔内に単一又は複数の受
信器が位置され、PRBS発振源が励起されて受信器に
向けて地震エネルギが伝播され、PRBS発振源と受信
器間で測定が直接行われるような、本発明に従う方法の
結果に基づくものである限り、他の分析的な方法を利用
することも可能である。Although the present invention has been described with reference to certain specific embodiments, it will be apparent that various modifications can be made without departing from the spirit and scope of the invention. For example, more holes can be used for inter-hole tomography according to the present invention. Also, if the respective positions are known, many relative positional relationships and variations of the positions between the PRBS oscillator source and the receiver can be used. Furthermore, various theories and relational expressions have been used to determine and calculate properties such as porosity, water permeability, water permeability profile, shear deformation rate and shear force. Located in the hole, the PR
Single or multiple receivers are located in a separate exploration hole away from the BS oscillator, the PRBS oscillator is excited and seismic energy is propagated towards the receiver, measured between the PRBS oscillator and the receiver It is also possible to use other analytical methods, as long as they are based on the results of the method according to the invention, such as is carried out directly.
【0067】[0067]
【発明の効果】以上説明した通り、本発明によれば、擬
似ランダムバイナリシークエンスコードを発振音響波に
用いて、相関により他の雑音と識別するようにしたの
で、杭打ちや重機の走行等の雑音がある場所でも、地層
の非破壊測定が可能となる上に、一回の測定可能距離が
数百メートルに及び、精度も著しく向上することができ
るという優れた効果を有する。As described above, according to the present invention, the pseudo random binary sequence code is used for the oscillating acoustic wave so as to distinguish it from other noises by the correlation. Even in a noisy place, the non-destructive measurement of the formation is possible, and the measurable distance per time is several hundred meters, and the accuracy is remarkably improved.
【0068】発明者の実験によれば、製鉄所内の、杭打
ち機や、ブルドーザが走行する雑音の多い場所におい
て、220m スパンまで測定可能であることが実際に確
認できた。According to the experiment conducted by the inventor, it was actually confirmed that it is possible to measure up to 220 m span in a noisy place where a pile driver or a bulldozer runs in a steel mill.
【図1】本発明で用いられる擬似ランダムバイナリシー
クエンスコード信号の例を示すタイムチャートFIG. 1 is a time chart showing an example of a pseudo random binary sequence code signal used in the present invention.
【図2】同じく受信した音響波と相関をとった後の音響
波の例を示す線図FIG. 2 is a diagram showing an example of an acoustic wave after having been correlated with an acoustic wave that is also received.
【図3】音響トモグラフィのための発振源と受信器の配
置を含む、本発明の一実施例を示すための、堆積層の垂
直断面を示す断面図FIG. 3 is a cross-sectional view showing a vertical cross-section of a deposited layer to illustrate one embodiment of the present invention, including an oscillator source and receiver arrangement for acoustic tomography.
【図4】本発明の一実施例で使用された装置を示すブロ
ック線図FIG. 4 is a block diagram showing an apparatus used in an embodiment of the present invention.
【図5】あるハイドロホンによって感知された初期デー
タを示す線図FIG. 5 is a diagram showing initial data sensed by a hydrophone.
【図6】他のハイドロホンによって感知された初期デー
タの例を示す線図FIG. 6 is a diagram showing an example of initial data detected by another hydrophone.
10…PRBS発振器 11、14…探査孔 12、13…ハイドロホン 21、32…PRBSコード発生器(PRBS Code
Generator) 30…周波数シンセサイザ(Frequency Synthesize
r) 36…圧電音響変換器(Source ) 38…ハイドロホンアレイ(Hydrophone Array) 44…信号増幅フィルタ(Signal Amplifier Filte
r ) 46…データ獲得ボックス(Data Acquisition Bo
x) 48…コンピュータ(Computer ) 50…ファンクションゼネレータ(Function Generat
or)10 ... PRBS oscillator 11, 14 ... Exploration hole 12, 13 ... Hydrophone 21, 32 ... PRBS code generator (PRBS Code)
(Generator) 30 ... Frequency Synthesizer (Frequency Synthesize)
r) 36 ... Piezoelectric acoustic transducer (Source) 38 ... Hydrophone Array 44 ... Signal amplification filter (Signal Amplifier Filter)
r) 46 ... Data acquisition box (Data Acquisition Bo)
x) 48 ... Computer (Computer) 50 ... Function Generator (Function Generator)
or)
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 FI 技術表示箇所 G01V 1/36 9406−2G 1/40 9406−2G (56)参考文献 特開 平4−198794(JP,A) 特開 昭60−108750(JP,A) 特開 昭57−173776(JP,A)─────────────────────────────────────────────────── ─── Continuation of the front page (51) Int.Cl. 6 Identification number Office reference number FI technical display location G01V 1/36 9406-2G 1/40 9406-2G (56) Reference JP-A-4-198794 ( JP, A) JP 60-108750 (JP, A) JP 57-173776 (JP, A)
Claims (12)
伸び、入口を有する、互いに所定距離だけ離れた少なく
とも一対の探査孔を形成し、 前記探査孔の既知の位置に受信器を配置し、 前記探査孔の他の1つの中の、前記受信器に関して既知
の位置に、擬似ランダムコード震動エネルギの発振源を
配置し、 前記擬似ランダムコードに従って該発振源を励起して、
該発振源から前記受信器に震動エネルギを伝播させ、 前記発振源の位置から前記受信器の位置に向かって伸び
る経路中の前記震動エネルギの擬似ランダムコード地震
波特性を測定することを特徴とする地層の物理特性の音
響波を用いた非破壊測定方法。1. At least a pair of exploration holes extending from the surface of the stratum to at least a predetermined depth and having inlets and separated from each other by a predetermined distance, and arranging a receiver at a known position of the exploration holes, Arranging a source of pseudo-random code seismic energy at a known location with respect to the receiver in another one of the exploration holes and exciting the source according to the pseudo-random code;
A stratum characterized by propagating seismic energy from the oscillation source to the receiver, and measuring a pseudo-random code seismic wave characteristic of the seismic energy in a path extending from the position of the oscillation source to the position of the receiver. Nondestructive measurement method using acoustic waves of physical properties of.
記受信器に対して移動し、該発振源を再び励起して、前
記受信器に向かって伸びる新しい経路中の擬似ランダム
コード地震波特性を更に測定することを特徴とする地層
の物理特性の音響波を用いた非破壊測定方法。2. The pseudo random code seismic wave characteristic according to claim 1, further comprising moving the oscillation source with respect to the receiver, re-exciting the oscillation source, and extending in a new path toward the receiver. A non-destructive measurement method using acoustic waves of physical properties of the formation, which is characterized by further measuring.
の探査孔を形成し、 少なくとも1つの前記探査孔内の既知位置に、多数の受
信器を互いに間隔を空けて配置し、 前記探査孔の他の1つの中の、前記受信器に関して既知
の位置に、擬似ランダムコード震動エネルギの発振源を
配置し、 前記擬似ランダムコードに従って該発振源を励起して、
該発振源から前記受信器に震動エネルギを伝播させ、 前記発振源の位置から前記受信器の位置に向かって伸び
る多数の経路中の前記震動エネルギの擬似ランダムコー
ド地震波特性を測定して記憶することを特徴とする地層
の物理特性の音響波を用いた非破壊測定方法。3. At least a pair of probe holes which are separated from each other by a predetermined distance, and a plurality of receivers are spaced apart from each other at known positions in at least one probe hole. A source of pseudo-random code seismic energy at a known location with respect to the receiver in one of the
Propagating seismic energy from the oscillation source to the receiver, measuring and storing pseudo random code seismic wave characteristics of the seismic energy in a number of paths extending from the position of the oscillation source toward the position of the receiver. Nondestructive measurement method using acoustic waves for physical properties of formations characterized by.
伸び、入口を有する、互いに所定距離だけ離れた少なく
とも一対の探査孔を形成し、 前記探査孔間の地層表面に実質的に沿う所定の既知の位
置に、複数の受信器を互いに間隔を空けて配置し、 少なくとも1つの前記探査孔内の既知位置に、多数の受
信器を互いに間隔を空けて配置し、 前記探査孔の他の1つの中の、前記受信器に関して既知
の位置に、擬似ランダムコード震動エネルギの発振源を
配置し、 前記擬似ランダムコードに従って該発振源を励起して、
該発振源から前記受信器に震動エネルギを伝播させ、 前記発振源の位置から前記受信器の位置に向かって伸び
る多数の経路中の前記震動エネルギの擬似ランダムコー
ド地震波特性を測定し、 前記発振源を前記受信器に対して移動し、 該発振源を再び励起して、前記受信器に向かって伸びる
新しい多数の経路中の擬似ランダムコード地震波特性を
更に測定することを特徴とする地層の物理特性の音響波
を用いた非破壊測定方法。4. A predetermined known structure which extends from the surface of the formation to at least a predetermined depth and which has at least a pair of exploration holes having inlets and separated from each other by a predetermined distance, and which is substantially along the stratum surface between the exploration holes. , A plurality of receivers spaced apart from each other, and a plurality of receivers spaced apart from each other at a known location within at least one of the probe holes, Inside, at a known position with respect to the receiver, placing an oscillating source of pseudo-random code seismic energy, exciting the oscillating source according to the pseudo-random code,
Propagating seismic energy from the oscillation source to the receiver, measuring pseudo random code seismic wave characteristics of the seismic energy in a number of paths extending from the position of the oscillation source to the position of the receiver, the oscillation source To the receiver and re-exciting the source to further measure the pseudo-random code seismic properties in a number of new paths extending towards the receiver. Non-destructive measurement method using acoustic waves of.
振源が圧電セラミック変換器であることを特徴とする地
層の物理特性の音響波を用いた非破壊測定方法。5. The non-destructive measurement method according to claim 4, wherein the oscillation source of the vibration energy is a piezoelectric ceramic converter.
実質的に平行であることを特徴とする地層の物理特性の
音響波を用いた非破壊測定方法。6. The non-destructive measurement method according to claim 4, wherein the exploration holes are substantially parallel to each other, using acoustic waves of physical properties of the formation.
に等間隔で配置されていることを特徴とする地層の物理
特性の音響波を用いた非破壊測定方法。7. The nondestructive measurement method according to claim 4, wherein the receivers are arranged at substantially equal intervals, using acoustic waves of physical properties of the formation.
に際して、伝播した震動エネルギの固有減衰及び速度を
測定することを特徴とする地層の物理特性の音響波を用
いた非破壊測定方法。8. The nondestructive measurement method using acoustic waves of physical characteristics of a stratum according to claim 4, wherein the characteristic attenuation and velocity of propagating vibration energy are measured when the seismic wave characteristics are measured.
特性が、縦波速度の周波数特性や固有減衰の周期数特性
を含むことを特徴とする地層の物理特性の音響波を用い
た非破壊測定方法。9. The non-destructive use of acoustic waves of physical properties of a stratum according to claim 4, wherein the measured seismic wave characteristics include frequency characteristics of longitudinal wave velocity and period number characteristics of intrinsic damping. Measuring method.
周波数が変化することを特徴とする地層の物理特性の音
響波を用いた非破壊測定方法。10. The non-destructive measurement method according to claim 4, wherein the frequency of the pseudo-random signal is changed, and the physical characteristics of the formation are acoustic waves.
コードのキャリア周波数が、ほぼ1、2、4、8及び1
0KHz の所で変化することを特徴とする地層の物理特
性の音響波を用いた非破壊測定方法。11. The carrier frequency of the pseudo-random code according to claim 10, wherein the carrier frequencies are approximately 1, 2, 4, 8 and 1.
A non-destructive measurement method using acoustic waves of the physical properties of the formation, which is characterized by changing at 0 kHz.
ぼ70−200m 離れていることを特徴とする地層の物
理特性の音響波を用いた非破壊測定方法。12. The non-destructive measurement method according to claim 4, wherein the exploration holes are separated from each other by approximately 70 to 200 m.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP21891992A JPH0820438B2 (en) | 1992-08-18 | 1992-08-18 | Nondestructive measurement method of physical properties of formation using acoustic wave |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP21891992A JPH0820438B2 (en) | 1992-08-18 | 1992-08-18 | Nondestructive measurement method of physical properties of formation using acoustic wave |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP9126157A Division JP2862171B2 (en) | 1997-05-16 | 1997-05-16 | Nondestructive method for measuring physical properties of formation using acoustic waves |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH06294793A JPH06294793A (en) | 1994-10-21 |
JPH0820438B2 true JPH0820438B2 (en) | 1996-03-04 |
Family
ID=16727379
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP21891992A Expired - Fee Related JPH0820438B2 (en) | 1992-08-18 | 1992-08-18 | Nondestructive measurement method of physical properties of formation using acoustic wave |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH0820438B2 (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10186048A (en) * | 1996-04-04 | 1998-07-14 | Yamamoto Tokuo | Measuring method and device for buried object, structure of stratum, and nature of deposit |
CA2255719C (en) * | 1997-03-17 | 2001-02-20 | Junichi Sakakibara | Device and method for transmitting acoustic wave into underground, for receiving the acoustic wave, and underground exploration method using above mentioned device |
US8031555B2 (en) | 2008-01-07 | 2011-10-04 | University Of Utah Research Foundation | Seismic location and communication system |
CN105607123B (en) * | 2016-01-15 | 2017-11-10 | 中国石油天然气股份有限公司 | Method and device for calculating seismic wave characteristic information of random pore medium model |
JP6815757B2 (en) * | 2016-06-09 | 2021-01-20 | 株式会社不動テトラ | Ground liquefaction evaluation method |
-
1992
- 1992-08-18 JP JP21891992A patent/JPH0820438B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH06294793A (en) | 1994-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Paillet et al. | Acoustic modes of propagation in the borehole and their relationship to rock properties | |
US5406530A (en) | Pseudo-random binary sequence measurement method | |
Kearey et al. | An introduction to geophysical exploration | |
EP0526501B1 (en) | Discrete-frequency multipole sonic logging methods and apparatus | |
Tang et al. | Fast inversion of formation permeability from Stoneley wave logs using a simplified Biot-Rosenbaum model | |
US7602669B2 (en) | Tube-wave seismic imaging | |
OA13110A (en) | Method and apparatus for determining the nature ofsubmarine reservoirs. | |
CN100456045C (en) | Method and system for measuring sea bed deposite sediment property by shallow stratigraphic section instrument | |
US6813566B2 (en) | Method of producing continuous, orthogonal signals and method of their use for examining and for detecting changes in a body | |
US5142500A (en) | Non-destructive method of measuring physical characteristics of sediments | |
JP2862171B2 (en) | Nondestructive method for measuring physical properties of formation using acoustic waves | |
Park et al. | Seismic characterization of geotechnical sites by multichannel analysis of surface waves (MASW) method | |
McCann et al. | Application of cross-hole seismic measurements in site investigation surveys | |
JPH0820438B2 (en) | Nondestructive measurement method of physical properties of formation using acoustic wave | |
Liu et al. | In situ measurement of seismic shear-wave absorption in the San Francisco Holocene Bay Mud by the pulse-broadening method | |
EP2577359B1 (en) | Method and system for determining the rigidity of a geological layer | |
Duke | Techniques for field measurement of shear wave velocity in soils | |
Beeston et al. | Shear wave velocities from down‐hole measurements | |
Lebedev et al. | Accuracy in Estimating the Parameters of a Layered Medium Using Coherent Vector Reception of a Rayleigh Surface Wave | |
Malovichko et al. | Active-passive array surface wave inversion and comparison to borehole logs in southeast Missouri | |
Muir et al. | Comparison of techniques for shear wave velocity and attenuation measurements | |
JP2000186319A (en) | Ground investigation method | |
Potylitsyn et al. | The Results of Observing the Seismoelectric Effect in a Gas Condensate Hydrocarbon Field Using a Source of Seismic Field | |
Lebedev et al. | Coherent seismoacoustics | |
Stoll et al. | Generating interface waves using a freely falling, instrumented source |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090304 Year of fee payment: 13 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100304 Year of fee payment: 14 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100304 Year of fee payment: 14 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110304 Year of fee payment: 15 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110304 Year of fee payment: 15 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120304 Year of fee payment: 16 |
|
LAPS | Cancellation because of no payment of annual fees |