JP7629376B2 - Exploration device, exploration method, and exploration system - Google Patents
Exploration device, exploration method, and exploration system Download PDFInfo
- Publication number
- JP7629376B2 JP7629376B2 JP2021146336A JP2021146336A JP7629376B2 JP 7629376 B2 JP7629376 B2 JP 7629376B2 JP 2021146336 A JP2021146336 A JP 2021146336A JP 2021146336 A JP2021146336 A JP 2021146336A JP 7629376 B2 JP7629376 B2 JP 7629376B2
- Authority
- JP
- Japan
- Prior art keywords
- points
- point
- amplitude
- waveform data
- vibration
- 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.)
- Active
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Description
本発明は、建設工事における地下埋設物の探査装置、探査方法、および、探査システムに関する。 The present invention relates to an exploration device, exploration method, and exploration system for underground buried objects in construction work.
都市土木や建築において、工事を始める前に、地下埋設物の探査が行われている。例えば、特許文献1には、地表面に設置し地中に対して圧電素子または磁歪素子による振動を生成する発振器と、地中の埋設物もしくは地層境界で反射し地表に達した反射波を検知する検知器と、検知した反射波によって埋設物もしくは地層境界の位置を推定する位置推定装置と、を備え、位置推定装置は、発振器の振動周波数を高周波から低周波へ変化させたスイープ波を印加し、位置を推定する反射法探査システムが開示されている。
In urban civil engineering and construction, underground buried objects are surveyed before construction begins. For example,
しかしながら、上記技術では、埋設物を特定するために必要である埋設物の形状まで十分把握することが難しいという問題があった。 However, the above technology has the problem that it is difficult to fully grasp the shape of buried objects, which is necessary to identify the buried objects.
そこで、本発明は上記の問題点等に鑑みて為されたもので、その課題の一例は、地中における埋設物の形状を把握できる探査装置等を提供することを目的とする。 The present invention was made in consideration of the above problems, and one example of the objective of the invention is to provide an exploration device that can grasp the shape of buried objects underground.
上記の課題を解決するために、請求項1に記載の発明は、起振装置によって地面が起振され、埋設物の大きさより小さい間隔である前記地面における複数の地点からの波形データを取得する取得手段と、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の合計走時を算出する合計走時算出手段と、前記合計走時に応じて、各前記1の受振点での前記波形データの振幅を重合する振幅重合手段と、前記重合された振幅に応じて、前記埋設物での反射点を特定する反射点特定手段と、を備え、異なる複数の位置の前記起振装置に起振されて得られた前記波形データによって算出される複数の前記反射点から前記埋設物の形状を特定し、前記合計走時算出手段が、前記埋設物の大きさより小さい間隔で設定された地下の散乱点を経由して、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の前記合計走時を算出し、前記振幅重合手段が、前記合計走時における前記1の受振点での前記波形データの振幅を、前記散乱点の振幅として、前記散乱点毎に、各受振点からの前記散乱点の振幅を重合し、前記反射点特定手段が、前記重合された散乱点の振幅に応じて、前記散乱点から前記反射点を特定することを特徴とする。 In order to solve the above problem, the present invention provides a method for obtaining waveform data from a plurality of points on the ground, the distance between which is smaller than the size of an embedded object, by exciting the ground with an excitation device; a total travel time calculation means for calculating a total travel time of elastic waves from an excitation point of the excitation device to one of the receiving points; an amplitude superimposition means for superimposing the amplitude of the waveform data at each of the receiving points according to the total travel time; and a reflection point identification means for identifying a reflection point at the embedded object according to the superimposed amplitude, and the waveform data obtained by exciting the ground with the excitation device at a plurality of different positions is stored in the storage means. the shape of the buried object is identified from the multiple reflection points calculated by a data set ; the total travel time calculation means calculates the total travel time of the elastic waves from the excitation point of the excitation device to one of the receiving points among the points, via underground scattering points set at intervals smaller than the size of the buried object; the amplitude stacking means sets the amplitude of the waveform data at the one receiving point during the total travel time as the amplitude of the scattering point, and stacks the amplitudes of the scattering points from each receiving point for each scattering point; and the reflection point identification means identifies the reflection point from the scattering points in accordance with the amplitude of the stacked scattering points .
また、請求項2に記載の発明は、請求項1に記載の探査装置において、前記反射点特定手段は、前記重合された振幅が最大である点を、前記反射点として特定することを特徴とする。 In a second aspect of the present invention, in the exploration device according to the first aspect, the reflection point identifying means identifies a point at which the superimposed amplitude is maximum as the reflection point.
また、請求項3に記載の発明は、請求項1または請求項2に記載の探査装置において、前記振幅重合手段が重合する波形データは、前記取得した波形データと前記起振装置による起振波形データとの相関データを計算し、前記相関データから所定の周波数範囲の成分の波形を抽出し、前記抽出した波形から前記起振装置による立ち上がり部分の波形をミュートして得られる波形データであることを特徴とする。
The invention described in
また、請求項4に記載の発明は、請求項1から請求項3のいずれか1項に記載の探査装置において、前記起振装置により起振され振動によって、前記埋設物の大きさ以上の間隔である複数地点からの広域波形データを取得する広域波形データ取得手段と、前記広域波形データから、前記埋設物の深さおよび地面上の位置を推定する埋設物位置推定手段と、を更に備えたことを特徴とする。
The invention described in
また、請求項5に記載の発明は、請求項4に記載の探査装置において、前記起振装置の前記異なる複数の位置において、互いに最も離れた前記起振装置の位置が、前記推定された前記埋設物の深さの2倍以上であることを特徴とする。
The invention described in
また、請求項6に記載の発明は、請求項1から5のいずれか1項に記載の探査装置において、前記取得手段が、前記地面に設置された受振センサから前記波形データを取得することを特徴とする。
The invention described in claim 6 is characterized in that in the exploration device described in any one of
また、請求項7に記載の発明は、請求項1から5のいずれか1項に記載の探査装置において、前記取得手段が、非接触の光学式の振動計から前記波形データを取得することを特徴とする。
The invention described in
また、請求項8に記載の発明は、取得手段が、起振装置によって地面が起振され、埋設物の大きさより小さい間隔である前記地面における複数の地点からの波形データを取得する取得ステップと、合計走時算出手段が、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の合計走時を算出する合計走時算出ステップと、振幅重合手段が、前記合計走時に応じて、各前記1の受振点での前記波形データの振幅を重合する振幅重合ステップと、反射点特定手段が、前記重合された振幅に応じて、前記埋設物での反射点を特定する反射点特定ステップと、を含み、異なる複数の位置の前記起振装置に起振されて得られた前記波形データによって算出される複数の前記反射点から前記埋設物の形状を特定し、前記合計走時算出手段が、前記埋設物の大きさより小さい間隔で設定された地下の散乱点を経由して、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の前記合計走時を算出し、前記振幅重合手段が、前記合計走時における前記1の受振点での前記波形データの振幅を、前記散乱点の振幅として、前記散乱点毎に、各受振点からの前記散乱点の振幅を重合し、前記反射点特定手段が、前記重合された散乱点の振幅に応じて、前記散乱点から前記反射点を特定することを特徴とする。 The invention described in claim 8 includes an acquisition step in which an acquisition means acquires waveform data from a plurality of points on the ground at intervals smaller than the size of an embedded object by exciting the ground with an excitation device, a total traveltime calculation step in which a total traveltime calculation means calculates a total traveltime of elastic waves from an excitation point of the excitation device to one of the points, an amplitude stacking step in which an amplitude stacking means stacks the amplitude of the waveform data at each of the one receiving points according to the total traveltime, and a reflection point identification step in which a reflection point identification means identifies a reflection point at the embedded object according to the stacked amplitude, the shape of the buried object is identified from the multiple reflection points calculated from the waveform data obtained by exciting the buried object at a predetermined frequency , the total travel time calculation means calculates the total travel time of the elastic waves from the excitation point of the excitation device to one of the receiving points among the points, via underground scattering points set at intervals smaller than the size of the buried object, the amplitude stacking means sets the amplitude of the waveform data at the one receiving point during the total travel time as the amplitude of the scattering point, and stacks the amplitudes of the scattering points from each receiving point for each scattering point, and the reflection point identification means identifies the reflection point from the scattering points in accordance with the amplitude of the stacked scattering points .
また、請求項9に記載の発明は、地面を起振させる起振装置と、前記起振装置によって地面が起振され振動を受振する受振装置と、探査装置と、を備える探査システムにおいて、前記探査装置は、前記起振装置によって地面が起振され、埋設物の大きさより小さい間隔である前記地面における複数の地点からの波形データを、前記受振装置から取得する取得手段と、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の合計走時を算出する合計走時算出手段と、前記合計走時に応じて、各前記1の受振点での前記波形データの振幅を重合する振幅重合手段と、前記重合された振幅に応じて、前記埋設物での反射点を特定する反射点特定手段と、を有し、異なる複数の位置の前記起振装置に起振されて得られた前記波形データによって算出される複数の前記反射点から前記埋設物の形状を特定し、前記合計走時算出手段が、前記埋設物の大きさより小さい間隔で設定された地下の散乱点を経由して、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の前記合計走時を算出し、前記振幅重合手段が、前記合計走時における前記1の受振点での前記波形データの振幅を、前記散乱点の振幅として、前記散乱点毎に、各受振点からの前記散乱点の振幅を重合し、前記反射点特定手段が、前記重合された散乱点の振幅に応じて、前記散乱点から前記反射点を特定することを特徴とする。 The invention described in claim 9 is an exploration system including an excitation device for vibrating the ground, a receiving device for receiving the vibrations generated by the excitation device, and an exploration device, the exploration device including: an acquisition means for acquiring waveform data from the receiving device when the ground is excited by the excitation device and the waveform data is from a plurality of points on the ground that are spaced apart from each other by an interval smaller than the size of an underground object; a total travel time calculation means for calculating a total travel time of elastic waves from the excitation point of the excitation device to one receiving point among the points; an amplitude superimposition means for superimposing the amplitude of the waveform data at each receiving point according to the total travel time; and a reflection point identification means for identifying a reflection point at the underground object according to the superimposed amplitude. and identifying a shape of the buried object from a plurality of reflection points calculated from the waveform data obtained by excitation by the excitation device at a plurality of different positions , the total travel time calculation means calculates the total travel time of the elastic wave from the excitation point of the excitation device to one of the receiving points among the positions, via underground scattering points set at intervals smaller than the size of the buried object, the amplitude stacking means sets the amplitude of the waveform data at the one receiving point during the total travel time as the amplitude of the scattering point, and stacks the amplitudes of the scattering points from each receiving point for each scattering point, and the reflection point identification means identifies the reflection point from the scattering points in accordance with the amplitude of the stacked scattering points .
本発明によれば、起振装置によって地面が起振され、埋設物の大きさより小さい間隔である地面における複数の地点からの波形データを取得し、起振点から地点うちの1の受振点までの弾性波の合計走時を算出し、合計走時に応じて、各1の受振点での波形データの振幅を重合し、重合された振幅に応じて、埋設物での反射点を特定し、異なる複数の位置の起振装置に起振されて得られた波形データによって算出される複数の反射点によって、地中における埋設物の形状を把握できる。 According to the present invention, the ground is excited by an excitation device, waveform data is obtained from multiple points on the ground spaced apart from each other by a distance smaller than the size of the buried object, the total travel time of the elastic waves from the excitation point to one of the receiving points is calculated, the amplitude of the waveform data at each receiving point is superimposed according to the total travel time, and a reflection point at the buried object is identified according to the superimposed amplitude, and the shape of the buried object underground can be grasped from the multiple reflection points calculated from the waveform data obtained by excitation by excitation devices at multiple different positions.
以下、図面を参照して本発明の実施形態について説明する。なお、以下に説明する実施の形態は、探査システム等に対して本発明を適用した場合の実施形態である。 Embodiments of the present invention will be described below with reference to the drawings. Note that the embodiments described below are embodiments in which the present invention is applied to an exploration system or the like.
[1.探査システムおよび各装置の構成および機能概要]
(1.1 探査システムの構成および機能概要)
[1. Overview of the configuration and functions of the exploration system and each device]
(1.1 Overview of the exploration system configuration and functions)
まず、本発明の一実施形態に係る探査システムの構成および概要機能について、図1を用いて説明する。 First, the configuration and general functions of an exploration system according to one embodiment of the present invention will be explained using FIG. 1.
図1は、本実施形態に係る探査システムの概略構成の一例を示す模式図である。 Figure 1 is a schematic diagram showing an example of the general configuration of an exploration system according to this embodiment.
図1に示すように、本実施形態に係る探査システム1は、地面Gに振動を与える起振装置10と、地面Gの振動を受振する受振装置20と、埋設物5を探査するめ受振装置20から波形データを解析する探査装置30と、を備える。
As shown in FIG. 1, the
建設工事における埋設物5は、1[m]以上の対象物で、例えば、下水管、共同溝、地下鉄等のトンネルである。また、埋設物5は、工事後に残置された基礎、スラブ、擁壁、耐圧盤等の躯体である。また、埋設物5は、工残置された杭、H鋼、鋼矢板等である。埋設物5は、例えば、20[m]程度までの深度に存在する。埋設物の深度は、2[m]以上が好ましい。
Buried
起振装置10は、振動子を有する。起振装置10は、数10[Hz]から数百[Hz]の周波数帯域の変調信号で地表から地中へ弾性波を伝搬させる。起振装置10は、埋設物5の直上だけでなく、十分に離れた位置(例えば、埋設物5の想定深度の2倍以上)にも、移動等して配置される。
The
起振装置10の振動子は、例えば、電動モータである。加振力は、450[N]である。加える振動は、縦波、横波でもよいが、縦波が好ましい。なお、振動子は、圧電素子、磁歪素子でもよい。
The vibrator of the
受振装置20は、複数の受振センサ20aから構成され、地面Gの振動を波形データに変換して探査装置30に出力する。各受振センサ20aは、格子状または測線状に、所定間隔に地面G上に配置される。図1に示すように、1方向に、所定間隔に配置された受振センサ20aが表示されている。格子状のセンサの代わりに、測線状のセンサを、受振センサ20aの並びの方向と平行でない方向にずらしながら、測定されてもよい。受振センサ20aが、地面に設置された受振センサの一例である。
The vibration receiving
受振センサ20aは、ばね式やMEMS等の機械式のセンサ、コイル式、ジオフォン等の電磁式のセンサ、圧電素子式、加速度ピックアップ等の圧電式のセンサ等の設置型のセンサである。 The vibration receiving sensor 20a is an installed sensor such as a mechanical sensor such as a spring type or MEMS, an electromagnetic sensor such as a coil type or geophone, or a piezoelectric sensor such as a piezoelectric element or acceleration pickup.
広域測定の場合、受振センサ20aの格子間隔は、例えば、50cmから数mの間隔で、予想される埋設物の大きさの半分ほどが好ましい。 For wide-area measurements, the grid spacing of the vibration receiving sensor 20a is preferably, for example, between 50 cm and several meters, approximately half the size of the expected buried object.
詳細測定の場合、受振センサ20aの格子間隔は、埋設物5の大きさより小さい間隔である。例えば、円筒状の埋設物5の直径より小さく、例えば、数cmから数10cmの高密度な間隔である。
For detailed measurements, the lattice spacing of the vibration receiving sensor 20a is smaller than the size of the buried
探査装置30は、起振装置10を所定の周波数で起振させ、受振装置20から波形データを取得する。
The
建設工事として、開削工事、シールド工事等の土木工事が挙げられる。また、建設工事として、山留工事、地盤改良工事、構造物解体工事、地中障害物撤去工事等が挙げられる。 Construction work includes civil engineering works such as open cut construction and shield construction. Other construction works include retaining walls, ground improvement work, structure demolition work, and underground obstacle removal work.
(1.2 探査装置30の構成および機能)
次に、探査装置30の構成および機能について、図2を用いて説明する。
(1.2 Configuration and Function of the Exploration Device 30)
Next, the configuration and functions of the
図2は、探査装置30の概要構成の一例を示すブロック図である。
Figure 2 is a block diagram showing an example of the general configuration of the
図2に示すように、コンピュータとして機能する探査装置30は、通信部31と、表示部32と、記憶部33と、操作部34と、入出力インターフェース部35と、制御部36とを有する。そして、制御部36と入出力インターフェース部35とは、システムバス37を介して接続されている。例えば、探査装置30は、パーソナルコンピュータやスマートフォンを含む携帯型無線電話機やPDA等の携帯端末である。
As shown in FIG. 2, the
通信部31は、起振装置10と接続し、起振装置10に、起振のパターンの変調信号を送信する。通信部31は、受振装置20の各受振センサ20aと接続し、受振センサ20aから、検出した波形データを取得する。
The
表示部32は、例えば、液晶表示素子または有機EL(Electro Luminescence)素子等によって構成されている。表示部32には、イメージング処理された結果が表示される。
The
記憶部33は、例えば、シリコンディスクドライブやハードディスクドライブ等からなる。記憶部33は、探査装置30を制御するための各種プログラム等を記憶したりする。各種プログラムは、オペレーティングシステム等が挙げられる。なお、各種プログラムは、例えば、無線通信網等のネットワークを介して取得されるようにしてもよいし、CD、DVD等の記録媒体に記録されてドライブ装置を介して読み込まれるようにしてもよい。
The
記憶部33は、受振装置20から送信されてきた波形データを、各受振センサ20aのチャンネル番号やセンサIDに関連付けて記憶する。記憶部33は、各受振センサ20aの配置情報を、受振センサ20aのチャンネル番号に関連付けて記憶する。
The
操作部34は、例えば、キーボードおよびマウス等によって構成されている。操作部34は、表示部32がタッチパネルのようなタッチスイッチ方式の表示パネルでもよい。
The
入出力インターフェース部35は、通信部31等と制御部36とのインターフェースである。
The input/output interface unit 35 is an interface between the
制御部36は、例えば、CPUと、ROMと、RAMと、を有する。制御部36は、CPUが、ROMや、RAMや、記憶部23に記憶された各種プログラムを読み出して実行する。
The
[2.探査システムの動作例]
次に、探査システムの動作例について、図3から図12Bを用いて説明する。
[2. Example of exploration system operation]
Next, an example of the operation of the exploration system will be described with reference to FIG. 3 to FIG. 12B.
図3に示すように、探査システム1は、広域測定を行う(ステップS1)。掘削工場等が行われる場所を中心付近に、起振装置10が設置される。起振装置10が設置された場所の脇または周りに、数10[cm]から数[m]の間隔の格子状または測線状に、受振装置20の各受振センサ20aが配置される。受振センサ20aの間隔は、想定される埋設物より大きさより大きいことが好ましい。
As shown in FIG. 3, the
次に、探査装置30の制御部36は、各受振センサ20aの配置情報を、受振センサ20aのチャンネル番号等に関連付けて記憶部33に記憶する。
Next, the
次に、探査装置30の制御部36は、起振装置10を起振させる。起振装置10が地面Gを起振する。制御部36は、受振装置20の各受振センサ20aから、広域波形データの一例である波形データを取得する。制御部36は、取得した受振の波形データを、受振センサ20aのチャンネル番号等に関連付けて記憶部33に記憶する。制御部36は、起振装置10を起震させた起震の波形データも記憶部33に記憶する。
Next, the
このように、探査装置30が、起振装置により起振され振動によって、前記埋設物の大きさ以上の間隔である複数地点からの広域波形データを取得する広域波形データ取得手段の一例として機能する。
In this way, the
次に、探査システム1は、埋設物5の大まかな位置を特定する(ステップS2)。具体的には、探査装置30の制御部36は、取得した波形データから、埋設物5の大まかな位置を特定する。例えば、広い範囲に複数の測線が設定され、図4に示すように、測定された波形データの散乱波が凸状に並ぶ波形のピークp(最も散乱波が速く到達した測定点)の測定点の位置が埋設物の位置とされる。深度は、散乱波の凸状分布と整合する地盤の速度を推定して埋設物のおおよその深度が求められる。地盤の速度は、図4に示すように、散乱波線とおおよそ整合する散乱波の時間から推定される。
Next, the
このように、探査装置30が、前記広域波形データから、前記埋設物の深さおよび地面上の位置を推定する埋設物位置推定手段の一例として機能する。
In this way, the
次に、探査システム1は、詳細測定を行う(ステップS3)。図1に示すように、特定された埋設物5の大まかな位置付近に、起振装置10が設置される。起振装置10が設置された場所の近傍である脇または周りに、数[cm]から数10[cm]の間隔の格子状または測線状に、受振装置20の各受振センサ20aが配置される。受振センサ20aの間隔は、想定される埋設物より大きさより小さい。
Next, the
次に、探査装置30の制御部36は、起振装置10の設置情報、および、各受振センサ20aの配置情報を、受振センサ20aのチャンネル番号等に関連付けて記憶部33に記憶する。なお、起振装置10および受振装置20が、GPSセンサを有してもよい。GPSセンサが取得した緯度経度情報が、起振装置10の設置情報および受振センサ20aの配置情報に含まれてもよい。
Next, the
次に、探査装置30の制御部36は、起振装置10を起振させる。起振装置10が地面Gを起振する。制御部36は、受振装置20の各受振センサ20aから波形データを取得する。制御部36は、取得した波形データを、受振センサ20aのチャンネル番号等に関連付けて記憶部33に記憶する。制御部36は、起振装置10を起震させた起震の波形データも記憶部33に記憶する。記憶部33において、受振の波形データと受振センサ20aの配置情報とは、関連付けられている。記憶部33において、起震の波形データと起振装置10の設置情報とは、関連付けられている。
Next, the
このように、探査装置30が、起振装置によって地面が起振され、埋設物の大きさより小さい間隔である前記地面における複数の地点からの波形データを取得する取得手段の一例として機能する。
In this way, the
次に、図5に示すように、起振装置10を移動させ、起振装置10が設置された場所の近傍である脇または周りに、数[cm]から数10[cm]の間隔の格子状または測線状に、受振装置20の各受振センサ20aが配置される。
Next, as shown in FIG. 5, the
次に、探査装置30の制御部36は、記憶部33に、起振装置10の設置情報、および、各受振センサ20aの配置情報を、受振センサ20aのチャンネル番号等に関連付けて記憶する。
Next, the
次に、探査装置30の制御部36は、起振装置10を起振させる。起振装置10が地面Gを起振する。制御部36は、受振装置20の各受振センサ20aから波形データを取得する。制御部36は、取得した受振の波形データを、受振センサ20aのチャンネル番号等に関連付けて記憶する。
Next, the
このように、起振装置10を数回移動させて、詳細測定を繰り返す。起振装置10の移動のさせ方は、埋設物5を横断するような配置が好ましい。さらに、起振装置10の前記異なる複数の位置において、互いに最も離れた起振装置10の位置が、推定された埋設物5の深さの2倍以上である。
In this manner, the
次に、探査システム1は、散乱波のイメージング処理する(ステップS4)。具体的には、探査装置30の制御部36は、記憶部33から、受振の波形データおよび起震の波形データを読み出し、後述の散乱波のイメージング処理のサブルーチンにて、イメージング処理を行い、図6Aおよび図6Bに示すように、埋設物5の散乱波のイメージング処理の画像40を求める。画像40において、横軸が地面Gでの距離、縦軸が深度である。画像40は、散乱波をイメージングした図であり、重合された振幅の一例である。なお、図6Aおよび図6Bは、起振装置10を移動させた異なる起振点での画像40である。
Next, the
次に、探査システム1は、反射点を抽出する(ステップS5)。具体的には、探査装置30の制御部36は、図7Aおよび図7Bに示すように、起振装置10を移動させた異なる起振点での各画像40において、振幅が最大である反射点41をそれぞれ抽出する。このように、制御部36は、イメージングされた散乱波のピークを反射点として抽出する。なお、画像40の各点が、深度マイグレーションの計算時に設定された散乱点である。
Next, the
このように、探査装置30が、重合された振幅に応じて、前記埋設物での反射点を特定する反射点特定手段の一例として機能する。探査装置30が、前記重合された散乱点の振幅に応じて、前記散乱点から前記反射点を特定する反射点特定手段の一例として機能する。探査装置30が、前記重合された振幅が最大である点を、前記反射点として特定する反射点特定手段の一例として機能する。
In this way, the
次に、探査システム1は、反射点の点群データを作成する(ステップS6)。具体的には、探査装置30の制御部36は、図8に示すように、抽出した反射点41を集めて、反射点41の点群データ42を作成する。埋設物5の上面の形状が点群データ42により表される。このように、起振点分の散乱点を2次元断面上にプロットして、図8のような埋設管探査結果を作成する。
Next, the
このように、探査装置30が、異なる複数の位置の前記起振装置に起振されて得られた前記波形データによって算出される複数の前記反射点から前記埋設物の形状を特定する。
In this way, the
次に、探査システム1は、埋設物探査の結果を表示する(ステップS7)。具体的には、探査装置30の制御部36は、表示部32に、埋設物探査の結果として、点群データ42を表示する。
Next, the
(2.1 散乱波のイメージング処理のサブルーチン)
次に、散乱波のイメージング処理のサブルーチンについて、図9から図10Dを用いて説明する。
(2.1 Subroutine for Scattered Wave Imaging Processing)
Next, the subroutine for the imaging process of the scattered wave will be described with reference to FIG. 9 to FIG. 10D.
図9に示すように、探査システム1は、起振の波形データと受振の波形データとの相互相関の波形を求める(ステップS10)。探査装置30の制御部36は、例えば、図10Aに示すような受振の波形データと、起振波形データとの相互相関を計算して、図10Bに示すように、相互相関の波形データを求める。相互相関の波形データが、取得した波形データと起振装置による起振波形データとの相関データの一例である。起振装置の変調信号と測定したデータの相互相関を計算することで、ノイズ軽減と距離分解能が向上し、地下の埋設物5からの散乱波とそれ以外の波とを分離する。
As shown in FIG. 9, the
次に、探査システム1は、バンドパスフィルタ処理を行う(ステップS11)。具体的には、制御部36は、相互相関の波形データから、所定値以下の低周波成分と、所定値以の高周波成分を除去し、図10Cに示すような波形データを求める。バンドパスフィルタ処理が、相関データから所定の周波数範囲の成分の波形を抽出する一例である。
Next, the
次に、探査システム1は、ミュート処理を行う(ステップS12)。具体的には、制御部36は、図10Dに示すように、起振装置10からの直接波に対応する立ち上がり部分の波形を取り除くため、ミュート関数を掛けることにより、ミュート処理を行う。なお、図10Dにおいて、ミュートされた箇所が示されている。このように、制御部36は、抽出した波形から起振装置による立ち上がり部分の波形をミュートする。
Next, the
ミュート関数として、シグモイド関数、ステップ関数等が挙げられる。ミュート関数は、起振装置10からの直接波の影響を受けている時間部分のところの振幅が低くなるような関数ならばよい。以上、相互相関からミュート処理までが、前処理である。
Examples of the mute function include a sigmoid function and a step function. Any mute function can be used as long as it reduces the amplitude during the time period affected by the direct wave from the
次に、探査システム1は、深度マイグレーション処理を行う(ステップS13)。具体的には、探査装置30の制御部36は、ミュート処理された各波形データに対して、後述の深度マイグレーション処理のサブルーチンにて、深度マイグレーション処理を行い、図6Aおよび図6Bに示すように、埋設物5の散乱波のイメージング処理の画像40を計算する。
Next, the
次に、探査システム1は、イメージング処理結果を出力する(ステップS14)。具体的には、探査装置30の制御部36は、深度マイグレーション処理結果の画像40を出力する。
Next, the
(2.2 深度マイグレーション処理のサブルーチン)
次に、深度マイグレーション処理のサブルーチンについて、図11から図12Bを用いて説明する。
(2.2 Depth Migration Processing Subroutine)
Next, the subroutine of the depth migration process will be described with reference to FIG. 11 to FIG. 12B.
図11に示すように、探査システム1は、地盤の弾性波速度を決定する(ステップS20)。具体的には、探査装置30の制御部36は、受振の波形データ、または、上記のように前処理された波形データ対して、残差ムービングアウト法や、コヒーレンシーインバージョン法等を適用して、速度モデルを構築し、地面Gの地盤の弾性波速度を決定する。このように、取得データを用いた弾性波速度解析を実施して、地下の弾性波速度を得ることで、時系列データから空間データに変換する。
As shown in FIG. 11, the
次に、探査システム1は、起振点から地下の散乱点への走時と、散乱点から受振点への走時とを計算する(ステップS21)。具体的には、探査装置30の制御部36は、図12Aに示すように、決定された弾性波速度に基づき、起振点Aから地下の散乱点Pへの走時T1と、散乱点Pから受振点Bへの走時T2とを計算する。
Next, the
次に、探査システム1は、合計走時を算出する(ステップS22)。具体的には、探査装置30の制御部36は、走時T1と走時T2との合計走時T1+T2を算出する。
Next, the
このように、探査装置30が、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の合計走時を算出する合計走時算出手段の一例として機能する。探査装置30が、前記埋設物の大きさより小さい間隔で設定された地下の散乱点を経由して、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の前記合計走時を算出る合計走時算出手段の一例として機能する。
In this way, the
次に、探査システム1は、起振点で起振して受振点で受振した波形の合計走時の振幅を、散乱点の振幅とする(ステップS23)。具体的には、探査装置30の制御部36は、受信点Bで測定された図10Dに示すようなミュート処理された波形データに対して、合計走時T1+T2での振幅を、散乱点Pの振幅とする。
Next, the
次に、探査システム1は、全ての地下の点と、起振点に対する全ての受振点に対して、ステップS20からステップS23の処理が行われたか判定する(ステップS24)。具体的には、探査装置30の制御部36は、図12Bのように、起振装置10が設置された起振点に対して、受振センサ20aが設置された全ての受振点と、所定の深度および所定の地面Gの距離の範囲にある全ての散乱点P1、P2等に対して、ステップS20からステップS23の処理が行われたか判定する。
Next, the
所定の深度および所定の地面Gの距離の範囲は、図6Aおよび図6Bの画像40に示すように、所定の深度、所定の距離の間である。所定の深度および所定の地面Gの距離の範囲に設定される、散乱点Pの計算上のメッシュは、散乱点Pの間隔は埋設物5の大きさより小さいことが好ましく、例えば、数cmから数10cmが好ましい。
The range of the predetermined depth and the distance of the ground G is between the predetermined depth and the predetermined distance, as shown in the image 40 in FIG. 6A and FIG. 6B. The calculation mesh of the scattering points P set in the range of the predetermined depth and the distance of the ground G is preferably such that the spacing between the scattering points P is smaller than the size of the buried
次に、探査システム1は、全ての振幅を足し合わせる(ステップS25)。具体的には、探査装置30の制御部36は、全ての地下の散乱点での振幅を足し合わせることで重合させて、図6Aおよび図6Bのような、散乱波のイメージング処理の画像40を求める。
Next, the
このように、探査装置30が、前記合計走時に応じて、各前記1の受振点での前記波形データの振幅を重合する振幅重合手段の一例として機能する。探査装置30が、前記合計走時における前記1の受振点での前記波形データの振幅を、前記散乱点の振幅として、前記散乱点毎に、各受振点からの前記散乱点の振幅を重合する振幅重合手段の一例として機能する。
In this way, the
上記、重合前深度マイグレーションの他に、重合前時間マイグレーション、重合後深度マイグレーション、重合後時間マイグレーションのアルゴリズムでもよい。 In addition to the pre-stack depth migration mentioned above, pre-stack time migration, post-stack depth migration, and post-stack time migration algorithms may also be used.
実施形態に係る探査装置30によれば、起振装置10によって地面Gが起振され、埋設物5の大きさより小さい間隔である地面Gにおける複数の地点からの波形データを取得し、起振点から地点うちの1の受振点までの弾性波の合計走時を算出し、合計走時に応じて、各1の受振点での波形データの振幅を重合し、重合された振幅に応じて、埋設物5での反射点を特定し、異なる複数の位置の起振装置10に起振されて得られた波形データによって算出される複数の反射点41によって、地中における埋設物の形状を把握できる。
According to the
抽出した反射点41の点群データ42の分布から埋設物の形状を把握できる。また、深度2[m]から20[m]の埋設物5を数10[cm]の分解能で位置と深度、形状を把握することにより、工事費用の削減や工期短縮が期待できる。
The shape of the buried object can be determined from the distribution of point cloud data 42 of the extracted reflection points 41. In addition, by determining the position, depth, and shape of buried
埋設物5の大きさより小さい間隔で設定された地下の散乱点Pを経由して、起振装置10の起振点Aから地点うちの1の受振点Bまでの弾性波の合計走時T1+T2を算出し、合計走時における1の受振点での波形データの振幅を、散乱点Pの振幅として、散乱点毎に、各受振点からの散乱点の振幅を重合し、重合された散乱点の振幅に応じて、散乱点から反射点41を特定する場合、起振装置10によって地面が起振され、埋設物5の大きさより小さい間隔である地面Gにおける複数の地点からの波形データを取得し、埋設物5の大きさより小さい間隔で設定された地下の散乱点Pを経由して、起振点Aから地点うちの1の受振点Bまでの弾性波の合計走時T1+T2を算出し、合計走時T1+T2における1の受振点Bでの波形データの振幅を散乱点Pの振幅として、散乱点P毎に各受振点からの散乱点Pの振幅を重合し、重合された散乱点Pの振幅に応じて散乱点Pから埋設物5の反射点41を特定し、異なる複数の位置の起振装置10から得られた波形データによって算出される複数の反射点41によって、地中における埋設物の形状を把握できる。
The total travel time T1 + T2 of the elastic wave from the excitation point A of the
重合された振幅が最大である点を、反射点41として特定する場合、反射点41が精度良く特定できる。 When the point where the superimposed amplitude is maximum is identified as reflection point 41, reflection point 41 can be identified with high accuracy.
重合する波形データは、取得した波形データと起振装置10による起振波形データとの相関データを計算し、相関データから所定の周波数範囲の成分の波形を抽出し、抽出した波形から起振装置10による立ち上がり部分の波形をミュートして得られる波形データである場合、起振装置10による直接の波形が除かれるため、反射点41が精度良く特定できる。
When the waveform data to be superimposed is obtained by calculating correlation data between the acquired waveform data and excitation waveform data generated by the
起振装置10により起振され振動によって、埋設物5の大きさ以上の間隔である複数地点からの広域波形データを取得し、広域波形データから、埋設物5の深さおよび地面上の位置を推定する場合、詳細測定の際、高密度に測定する場所を絞れるので、効率よく測定ができて、測定の時間短縮を図ることができる。
By using the vibration generated by the
起振装置10の異なる複数の位置において、互いに最も離れた起振装置10の位置が、推定された埋設物5の深さの2倍以上である場合、埋設物5の位置に応じて、地下の局所的に分布する埋設物5を対象に絞って、埋設物5から散乱する弾性波を起振点付近で高密度に測定できる。そのため、効率よく測定ができて、測定の時間短縮を図ることができる。埋設物5から散乱する弾性波を起振点付近で高密度に測定したデータにイメージング処理を適用することにより、深度2[m]から20[m]の埋設物5を数10[cm]の分解能で位置と深度、形状を把握することができる。
When the position of the
地面Gに設置された受振センサ20aから波形データを取得する場合、受振センサ20aを起振点周辺に高密度に配置することで、埋設物5の形状をイメージするための散乱波データが取得できる。
When acquiring waveform data from a vibration receiving sensor 20a installed on the ground G, scattering wave data for imaging the shape of the buried
(第2実施形態)
次に、第2実施形態について図13を用いて説明する。なお、前記第1実施形態と同一または対応する部分には、同一の符号を用いて異なる構成および動作のみを説明する。その他の実施形態および変形例も同様とする。
Second Embodiment
Next, a second embodiment will be described with reference to Fig. 13. Note that the same reference numerals are used for the same or corresponding parts as those in the first embodiment, and only different configurations and operations will be described. The same applies to the other embodiments and modifications.
本実施形態では、図13に示すように、設置型の受振装置20の代わりに、非接触の受振装置21を用いる。
In this embodiment, as shown in FIG. 13, a non-contact receiving device 21 is used instead of a
受振装置21は、レーザ式、レーザドップラ等の光学式、または、レーダ式等の電磁波式の非接触の振動計である。 The vibration receiving device 21 is a non-contact vibration meter that is optical, such as a laser type or a laser Doppler type, or electromagnetic wave type, such as a radar type.
例えば、受振装置21をクレーン車のクレーンの先に設置して高所から、広域測定が行われる。 For example, the receiver 21 can be installed at the tip of a crane on a crane truck, allowing wide-area measurements to be taken from a high position.
詳細測定では、受振装置21の高さを、広域測定を低くして、詳細測定が行われる。第1実施形態のように、起振装置10を移動させて、詳細測定が行われる。受振装置21の画像の解像度に応じて波形データが得られるが、数cmから数10cmの間隔で、波形データが取得されてもよい。
In the detailed measurement, the height of the vibration receiving device 21 is lowered from the wide-area measurement, and the detailed measurement is performed. As in the first embodiment, the
このように、探査装置30が、非接触の光学式の振動計から前記波形データを取得する前記波形データを取得する取得手段の一例として機能する。
In this way, the
非接触の光学式の振動計から波形データを取得する場合、受振センサ20aの設置作業を省くことができる。 When waveform data is acquired from a non-contact optical vibrometer, the installation work of the vibration receiving sensor 20a can be omitted.
なお、測定レイアウトを鉛直から水平にして、トンネルの前方探査に利用することができる。例えば、トンネルの掘削の掘削面に起振装置10を設置し、受振装置20または受振装置21から波形データを取得することで、掘削先の埋設物を探査する。
The measurement layout can be changed from vertical to horizontal and used for tunnel forward exploration. For example,
さらに、本発明は、上記各実施形態に限定されるものではない。上記各実施形態は、例示であり、本発明の特許請求の範囲に記載された技術的思想と実質的に同一な構成を有し、同様な作用効果を奏するものは、いかなるものであっても本発明の技術的範囲に包含される。 Furthermore, the present invention is not limited to the above-mentioned embodiments. The above-mentioned embodiments are merely examples, and anything that has substantially the same configuration as the technical idea described in the claims of the present invention and exhibits similar effects is included within the technical scope of the present invention.
1:探査システム
5:埋設物
10: 起振装置
20、21:受振装置
30:探査装置
41:反射点
G:地面
1: exploration system 5: buried object 10:
Claims (9)
前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の合計走時を算出する合計走時算出手段と、
前記合計走時に応じて、各前記1の受振点での前記波形データの振幅を重合する振幅重合手段と、
前記重合された振幅に応じて、前記埋設物での反射点を特定する反射点特定手段と、
を備え、
異なる複数の位置の前記起振装置に起振されて得られた前記波形データによって算出される複数の前記反射点から前記埋設物の形状を特定し、
前記合計走時算出手段が、前記埋設物の大きさより小さい間隔で設定された地下の散乱点を経由して、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の前記合計走時を算出し、
前記振幅重合手段が、前記合計走時における前記1の受振点での前記波形データの振幅を、前記散乱点の振幅として、前記散乱点毎に、各受振点からの前記散乱点の振幅を重合し、
前記反射点特定手段が、前記重合された散乱点の振幅に応じて、前記散乱点から前記反射点を特定することを特徴とする探査装置。 an acquisition means for acquiring waveform data from a plurality of points on the ground, the distance between the points being smaller than the size of the buried object, by an excitation device;
A total travel time calculation means for calculating a total travel time of an elastic wave from a vibration generating point of the vibration generating device to one of the vibration receiving points of the points;
an amplitude superposition means for superposing the amplitude of the waveform data at each of the one receiving points according to the total travel time;
a reflection point identifying means for identifying a reflection point at the buried object in accordance with the superimposed amplitude;
Equipped with
Identifying a shape of the buried object from a plurality of reflection points calculated from the waveform data obtained by exciting the vibration device at a plurality of different positions ;
the total travel time calculation means calculates the total travel time of the elastic waves from the excitation point of the excitation device to one of the receiving points among the points via underground scattering points set at intervals smaller than the size of the buried object;
The amplitude superimposition means superimposes the amplitude of the waveform data at the one receiving point during the total travel time as the amplitude of the scattering point, and superimposes the amplitude of the scattering point from each receiving point for each scattering point;
The exploration device , characterized in that the reflection point identifying means identifies the reflection point from the scattering points in accordance with the amplitude of the superimposed scattering points .
前記広域波形データから、前記埋設物の深さおよび地面上の位置を推定する埋設物位置推定手段と、
を更に備えたことを特徴とする請求項1から請求項3のいずれか1項に記載の探査装置。 a wide-area waveform data acquisition means for acquiring wide-area waveform data from a plurality of points spaced apart from each other by the vibration generated by the vibration generating device, the wide-area waveform data being acquired at intervals equal to or greater than the size of the buried object;
a buried object position estimating means for estimating the depth and the position on the ground of the buried object from the wide-area waveform data;
The exploration device according to any one of claims 1 to 3 , further comprising:
合計走時算出手段が、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の合計走時を算出する合計走時算出ステップと、
振幅重合手段が、前記合計走時に応じて、各前記1の受振点での前記波形データの振幅を重合する振幅重合ステップと、
反射点特定手段が、前記重合された振幅に応じて、前記埋設物での反射点を特定する反射点特定ステップと、
を含み、
異なる複数の位置の前記起振装置に起振されて得られた前記波形データによって算出される複数の前記反射点から前記埋設物の形状を特定し、
前記合計走時算出手段が、前記埋設物の大きさより小さい間隔で設定された地下の散乱点を経由して、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の前記合計走時を算出し、
前記振幅重合手段が、前記合計走時における前記1の受振点での前記波形データの振幅を、前記散乱点の振幅として、前記散乱点毎に、各受振点からの前記散乱点の振幅を重合し、
前記反射点特定手段が、前記重合された散乱点の振幅に応じて、前記散乱点から前記反射点を特定することを特徴とする探査方法。 an acquiring step in which the acquiring means acquires waveform data from a plurality of points on the ground, the intervals of which are smaller than the size of the buried object, by exciting the ground with a vibration device;
A total travel time calculation step in which a total travel time calculation means calculates a total travel time of an elastic wave from an excitation point of the excitation device to one of the receiving points of the points;
an amplitude stacking step in which an amplitude stacking means stacks the amplitude of the waveform data at each of the one receiving points according to the total travel time;
a reflection point identifying step in which a reflection point identifying means identifies a reflection point at the buried object in accordance with the superimposed amplitude;
Including,
Identifying a shape of the buried object from a plurality of reflection points calculated from the waveform data obtained by exciting the vibration device at a plurality of different positions ;
the total travel time calculation means calculates the total travel time of the elastic waves from the excitation point of the excitation device to one of the receiving points among the points via underground scattering points set at intervals smaller than the size of the buried object;
The amplitude superimposition means superimposes the amplitude of the waveform data at the one receiving point during the total travel time as the amplitude of the scattering point, and superimposes the amplitude of the scattering point from each receiving point for each scattering point;
The exploration method, wherein the reflection point identifying means identifies the reflection point from the scattering points in accordance with the amplitude of the superimposed scattering points .
前記探査装置は、
前記起振装置によって地面が起振され、埋設物の大きさより小さい間隔である前記地面における複数の地点からの波形データを、前記受振装置から取得する取得手段と、
前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の合計走時を算出する合計走時算出手段と、
前記合計走時に応じて、各前記1の受振点での前記波形データの振幅を重合する振幅重合手段と、
前記重合された振幅に応じて、前記埋設物での反射点を特定する反射点特定手段と、
を有し、
異なる複数の位置の前記起振装置に起振されて得られた前記波形データによって算出される複数の前記反射点から前記埋設物の形状を特定し、
前記合計走時算出手段が、前記埋設物の大きさより小さい間隔で設定された地下の散乱点を経由して、前記起振装置の起振点から前記地点うちの1の受振点までの弾性波の前記合計走時を算出し、
前記振幅重合手段が、前記合計走時における前記1の受振点での前記波形データの振幅を、前記散乱点の振幅として、前記散乱点毎に、各受振点からの前記散乱点の振幅を重合し、
前記反射点特定手段が、前記重合された散乱点の振幅に応じて、前記散乱点から前記反射点を特定することを特徴とする探査システム。 An exploration system including a vibration generating device that generates vibrations on the ground, a vibration receiving device that receives vibrations generated by the vibration generating device, and an exploration device,
The exploration device is
an acquisition means for acquiring waveform data from the vibration receiving device, the waveform data being obtained from a plurality of points on the ground that are spaced apart from each other by an interval smaller than the size of the buried object when the ground is excited by the vibration generating device;
A total travel time calculation means for calculating a total travel time of an elastic wave from a vibration generating point of the vibration generating device to one of the vibration receiving points of the points;
an amplitude superposition means for superposing the amplitude of the waveform data at each of the one receiving points according to the total travel time;
a reflection point identifying means for identifying a reflection point at the buried object in accordance with the superimposed amplitude;
having
Identifying a shape of the buried object from a plurality of reflection points calculated from the waveform data obtained by exciting the vibration device at a plurality of different positions ;
the total travel time calculation means calculates the total travel time of the elastic waves from the excitation point of the excitation device to one of the receiving points among the points via underground scattering points set at intervals smaller than the size of the buried object;
The amplitude superimposition means superimposes the amplitude of the waveform data at the one receiving point during the total travel time as the amplitude of the scattering point, and superimposes the amplitude of the scattering point from each receiving point for each scattering point;
An exploration system , characterized in that the reflection point identifying means identifies the reflection point from the scattering points according to the amplitude of the superimposed scattering points .
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021146336A JP7629376B2 (en) | 2021-09-08 | 2021-09-08 | Exploration device, exploration method, and exploration system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2021146336A JP7629376B2 (en) | 2021-09-08 | 2021-09-08 | Exploration device, exploration method, and exploration system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2023039260A JP2023039260A (en) | 2023-03-20 |
| JP7629376B2 true JP7629376B2 (en) | 2025-02-13 |
Family
ID=85600453
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2021146336A Active JP7629376B2 (en) | 2021-09-08 | 2021-09-08 | Exploration device, exploration method, and exploration system |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7629376B2 (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000178955A (en) | 1998-12-15 | 2000-06-27 | Shimizu Corp | Investigation method and investigation device for underground structures |
| JP2004279064A (en) | 2003-03-12 | 2004-10-07 | Yuzuru Ashida | 3D reflection data acquisition method using S wave |
| JP2009186449A (en) | 2008-02-04 | 2009-08-20 | Univ Of Tsukuba | Two-dimensional image data generation method and exploration device |
| JP2012042421A (en) | 2010-08-23 | 2012-03-01 | Toin Gakuen | Detection method employing sonic wave, non-contact acoustic detection system, program used for system, and recording medium with program recorded thereon |
-
2021
- 2021-09-08 JP JP2021146336A patent/JP7629376B2/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000178955A (en) | 1998-12-15 | 2000-06-27 | Shimizu Corp | Investigation method and investigation device for underground structures |
| JP2004279064A (en) | 2003-03-12 | 2004-10-07 | Yuzuru Ashida | 3D reflection data acquisition method using S wave |
| JP2009186449A (en) | 2008-02-04 | 2009-08-20 | Univ Of Tsukuba | Two-dimensional image data generation method and exploration device |
| JP2012042421A (en) | 2010-08-23 | 2012-03-01 | Toin Gakuen | Detection method employing sonic wave, non-contact acoustic detection system, program used for system, and recording medium with program recorded thereon |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2023039260A (en) | 2023-03-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20130094326A1 (en) | System and method for seismological sounding | |
| JP6902614B2 (en) | Internal structure grasping system | |
| US11841265B2 (en) | Heterogeneous subsurface imaging systems and methods | |
| WO2009032996A2 (en) | Seismic resonance imaging | |
| JP6531934B2 (en) | Hybrid surface wave search method and hybrid surface wave search system | |
| JP2014106128A (en) | Method of measuring natural ground elastic wave velocity | |
| JP7522504B2 (en) | Seismic exploration system and method | |
| US20200116881A1 (en) | Heterogeneous subsurface imaging systems and methods | |
| JP2022000666A (en) | Pit face front investigation method | |
| JP2013044612A (en) | Detection method and device of underground piping damage position | |
| Lancelle | Distributed acoustic sensing for imaging near-surface geology and monitoring traffic at Garner Valley, California | |
| JP7629376B2 (en) | Exploration device, exploration method, and exploration system | |
| JP2015092180A (en) | Geological exploration system | |
| JP2862171B2 (en) | Nondestructive method for measuring physical properties of formation using acoustic waves | |
| CN113126161B (en) | Method and system for predicting cave depth and size of karst cave based on shock waves | |
| JP2012117902A (en) | Ground survey device and ground survey method | |
| CN104020488A (en) | Wireless-distributed-type elastic wave reflector detection device, system and method | |
| CN108572390B (en) | Detection method for predicting shallow caves using surface wave spectral perturbations | |
| JPH09280906A (en) | Position locating system for falling rock source | |
| JP2006078243A (en) | Nondestructive shape diagnosis method and apparatus for buried object | |
| JP5679804B2 (en) | Geological exploration system | |
| JP3173790U (en) | Ground survey equipment | |
| JP2006275914A (en) | Method of imaging subsurface structure by reflection method seismic exploration | |
| Mihalache | Prediction of ground vibrations induced by impact driving of dolphin piles in Caland Canal | |
| CN113759422A (en) | Underground abnormal body detection method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20240405 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20241029 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20241031 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20241220 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20250121 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20250131 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7629376 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |