JP2020165925A - Face forward velocity estimation method in face forward survey system - Google Patents
Face forward velocity estimation method in face forward survey system Download PDFInfo
- Publication number
- JP2020165925A JP2020165925A JP2019069126A JP2019069126A JP2020165925A JP 2020165925 A JP2020165925 A JP 2020165925A JP 2019069126 A JP2019069126 A JP 2019069126A JP 2019069126 A JP2019069126 A JP 2019069126A JP 2020165925 A JP2020165925 A JP 2020165925A
- Authority
- JP
- Japan
- Prior art keywords
- face
- wave
- elastic
- elastic wave
- waveform
- 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.)
- Pending
Links
Images
Abstract
Description
本発明は、トンネル等の掘削を行う場合の切羽前方探査における、切羽前方速度推定方法に関するものである。 The present invention relates to a method for estimating the forward velocity of the face in the forward exploration of the face when excavating a tunnel or the like.
トンネル工事においては、トンネルの掘削や支保を安全かつ効率良く施工するために、地山の弾性波速度分布を推定することが行われ、また、掘削面である切羽前方の地質変化を把握し、その対策を講じることが安全かつ迅速な掘削作業を行う上で必要となるため、切羽前方の探査が行われる。 In tunnel construction, in order to construct tunnel excavation and support safely and efficiently, the elastic wave velocity distribution of the ground is estimated, and the geological change in front of the face, which is the excavation surface, is grasped. Since it is necessary to take such measures for safe and quick excavation work, exploration in front of the face is carried out.
切羽前方の探査を行う手法については、従来から種々の手法が提案されている。切羽掘削に支障を与えることなく、切羽前方の地質構造の３次元的評価を可能にする切羽前方探査装置として、特開２００１−９９９４５号公報（特許文献１）に開示された装置がある。特許文献１では、トンネル坑内において発生させた弾性波が切羽前方の弾性波反射面において鏡面反射した弾性波（反射波）を、トンネル坑内に設置した複数の地震計により計測し、３次元空間に反射面を仮定し、測定波形を重合処理することにより、反射面を推定する方法が提案されている。
Various methods have been conventionally proposed for exploration in front of the face. As a face forward exploration device that enables three-dimensional evaluation of the geological structure in front of the face without hindering face excavation, there is a device disclosed in Japanese Patent Application Laid-Open No. 2001-99945 (Patent Document 1). In
また、弾性波の３成分受振によるトンネル切羽前方の高精度イメージングに関する技術として、非特許文献１に示されたものがあり、この非特許文献１では、トンネル坑内において発生させた弾性波が切羽前方の弾性波反射面において鏡面反射した弾性波（反射波）を、トンネル坑内に設置した３成分地震計により計測し、測定された波形を重合処理した後、３成分データにより重み付けを行うことで、３次元的に反射面を推定する方法が提案されている。
Further, as a technique related to high-precision imaging of the front of the tunnel face by receiving three components of elastic waves, there is one shown in Non-Patent
さらに、トンネル坑内切羽前方の反射面に凹凸があっても、反射面を３次元的に推定可能にする切羽前方速度推定方法として、特開２０１５−１５８４３７号公報（特許文献２）に開示された技術がある。特許文献２では、発破により切羽から発生する弾性波を３成分地震計で計測して、多成分記録装置に記録することを切羽掘削毎に繰り返し、切羽掘削毎に記録された弾性波から直達波を抽出して、この直達波に基づいて切羽のトンネル坑内の弾性波速度分布を推定し、切羽掘削毎に記録された弾性波から反射波を抽出して、この反射波に基づいて反射・回折点及び反射面の３次元位置空間を推定する方法が提案されている。
Further, as a method for estimating the forward velocity of the face, which makes it possible to estimate the reflective surface three-dimensionally even if the reflective surface in front of the face in the tunnel is uneven, it is disclosed in Japanese Patent Application Laid-Open No. 2015-158437 (Patent Document 2). There is technology. In
しかしながら、上述した特許文献１及び非特許文献１に開示された技術を含むこれまでの切羽前方探査技術においては、トンネルの施工管理の面で十分に成果が挙がっているとはいえなかった。
However, it cannot be said that the conventional face forward exploration techniques including the techniques disclosed in
一般に、トンネルの施工管理の上では下記のような要望に応える必要がある。
Ａ．トンネルに損傷を与えないこと、
Ｂ．工事を出来るだけ中断しないこと、
Ｃ．切羽前方探査の探査結果を迅速に入手できること、
Ｄ．追加のコスト、及び手間をかけないこと、
Ｅ．施工現場に常駐する一般技術者操作できること、
等である。
In general, it is necessary to meet the following requests in tunnel construction management.
A. Do not damage the tunnel,
B. Do not interrupt the construction as much as possible,
C. To be able to quickly obtain the exploration results of the face forward exploration,
D. No additional cost and hassle,
E. Being able to operate a general engineer resident at the construction site,
And so on.
本発明は、上記従来からの要望に応えるためになされたものであって、トンネルの施工管理がより一層効率よく行える切羽前方探査における切羽前方速度推定方法を提供することを目的としている。 The present invention has been made in order to meet the above-mentioned conventional demands, and an object of the present invention is to provide a method for estimating the forward speed of a face in forward exploration of a face, which enables more efficient construction management of a tunnel.
上記課題を解決するため、本発明の切羽前方探査における前方速度推定方法は、切羽から離れた後方所定の位置に、発破信号検知手段及び弾性波記録手段を設置すると共に、１または複数の地震計をトンネル坑内壁面に設置して、トンネル坑内の切羽に発破孔を設けて爆薬を装填し、発破器から発生し爆薬を起爆させる発破信号を計測し、爆薬を起爆させた際に前記切羽から発生する弾性波を前記地震計でそれぞれ受振、計測し、前記弾性波記録手段に記録することを切羽掘削毎に繰り返し、切羽掘削毎に前記弾性波記録手段に記録された弾性波から反射波を抽出して、当該反射波に基づいて反射面の３次元位置を推定する切羽前方速度推定方法であって、切羽掘削毎に前記弾性波記録手段に記録された弾性波から直達波と反射波とを分離抽出する弾性波分離抽出ステップと、分離抽出された弾性波についてエンベロープ計算を行うエンベロープ演算ステップと、エンベロープ計算された弾性波についてピーク抽出を行うピーク抽出ステップと、ピーク抽出された直達波のピーク値で反射波振幅を基準化する弾性波基準化ステップと、反射波のピーク付近の波形を切り出す波形切出しステップと、反射波の反射係数をエンベロープの比で算出する反射係数演算ステップと、相関係数を計算する相関係数演算ステップと、相関係数の正負により、反射波の位相を判定する位相を判定ステップと、反射係数から速度分布を算出し、当該速度分布から弾性波の反射面の３次元位置を推定する反射面推定ステップと、を有することを特徴とするものである。 In order to solve the above problems, in the forward speed estimation method in the forward exploration of the face of the present invention, a burst signal detecting means and an elastic wave recording means are installed at a predetermined position behind the face, and one or more seismometers are installed. Is installed on the wall surface of the tunnel mine, a rupture hole is provided in the face inside the tunnel mine, and the explosive is loaded. The elastic waves to be generated are received and measured by the seismometer, and recorded in the elastic wave recording means is repeated for each face excavation, and the reflected waves are extracted from the elastic waves recorded in the elastic wave recording means for each face excavation. Then, it is a face forward velocity estimation method that estimates the three-dimensional position of the reflecting surface based on the reflected wave, and the direct wave and the reflected wave are separated from the elastic wave recorded in the elastic wave recording means for each face excavation. Elastic wave separation and extraction step for separation and extraction, envelope calculation step for performing envelope calculation for separated and extracted elastic wave, peak extraction step for peak extraction for envelope-calculated elastic wave, and peak of peak-extracted direct wave The phase relationship between the elastic wave standardization step that standardizes the reflected wave amplitude by the value, the waveform cutting step that cuts out the waveform near the peak of the reflected wave, and the reflection coefficient calculation step that calculates the reflection coefficient of the reflected wave by the ratio of the envelope. The correlation coefficient calculation step for calculating the number, the phase for determining the phase of the reflected wave based on the positive and negative of the correlation coefficient, and the velocity distribution for calculating the velocity distribution from the reflection coefficient, and the reflective surface of the elastic wave from the velocity distribution. It is characterized by having a reflection surface estimation step for estimating a three-dimensional position.
本発明の１つの態様は、前記弾性波分離抽出を行った後に、受振波形を積分して変位波形に変換し、幾何減衰補正するステップと、初動時刻をそろえるようにシフト（０ウェイ変換）して直達波の振幅の変位波形を得るステップと、得られた変位波形データをスタックするステップと、を含み、その次段階でエンベロープ演算ステップを実行する、ことを特徴とするものである。 In one aspect of the present invention, after performing the elastic wave separation extraction, the vibration receiving waveform is integrated and converted into a displacement waveform, and the step of correcting the geometrical attenuation is shifted (0-way conversion) so that the initial motion time is aligned. It is characterized by including a step of obtaining a displacement waveform of the amplitude of the direct wave and a step of stacking the obtained displacement waveform data, and executing an envelope calculation step in the next step.
本発明の別の態様は、前記弾性波分離抽出を行った後に、受振波形を積分して変位波形に変換し、幾何減衰補正するステップと、変位波形に対して２ウェイ変換を行い、この変換により反射波の振幅の変位波形を得るステップと、この変位波形データにフィルター処理を施すとともにスタックするステップと、を含み、その次段階でエンベロープ演算ステップを実行する、ことを特徴とするものである。 In another aspect of the present invention, after performing the elastic wave separation extraction, a step of integrating the vibration receiving waveform and converting it into a displacement waveform to correct geometrical attenuation, and performing a 2-way conversion on the displacement waveform are performed, and this conversion is performed. It is characterized in that it includes a step of obtaining a displacement waveform of the amplitude of the reflected wave and a step of filtering and stacking the displacement waveform data, and an envelope calculation step is executed in the next step. ..
本発明の第１の特徴の切羽前方速度推定方法及びこれに用いる切羽前方速度推定システムによれば、簡易なアルゴリズムにより解析を容易に行い、切羽前方速度を推定することができる。また、これらの結果として、本発明によれば、装置コストを抑制しつつ、より簡単なアルゴリズムによるデータ解析でトンネル坑内の切羽前方の反射面を３次元的に精度良く推定し得る切羽前方速度推定方法及びこれに用いる切羽前方速度推定システムを実現することができる。 According to the face forward velocity estimation method of the first feature of the present invention and the face forward velocity estimation system used therefor, analysis can be easily performed by a simple algorithm and the face forward velocity can be estimated. Further, as a result of these, according to the present invention, the front velocity of the face can be estimated three-dimensionally and accurately by data analysis by a simpler algorithm while suppressing the cost of the device. A method and a face forward velocity estimation system used for this can be realized.
本発明者らは、切羽での掘削発破を起振源とした切羽前方速度推定システム（トンネルフェイステスター：ＴＦＴ探査）を開発し、トンネル掘削現場への適用を重ねてきた。ＴＦＴ探査の内容は、国土交通省新技術情報提供システム（登録番号：ＴＨ−１７０００３−Ａ）で参照することができる。 The present inventors have developed a face forward velocity estimation system (tunnel face tester: TFT exploration) using the blasting of excavation at the face as a trigger source, and have repeatedly applied it to a tunnel excavation site. The contents of the TFT exploration can be referred to by the Ministry of Land, Infrastructure, Transport and Tourism's new technology information provision system (registration number: TH-170003-A).
以下の説明では、切羽前方に探知された速度不連続面の向こう側（奥側）の速度を推定する方法として考案した処理を組み合わせたアルゴリズムを提示し、その実用性ないしは有用性についてシミュレーション波形を用いて検討する。また、上記アルゴリズムをＴＦＴ探査システムに組み込んで実データへ適用した切羽前方速度推定方法を提案する。 In the following explanation, we present an algorithm that combines processes devised as a method for estimating the velocity on the other side (back side) of the velocity discontinuity detected in front of the face, and simulate the practicality or usefulness of the algorithm. Consider using. We also propose a face forward velocity estimation method in which the above algorithm is incorporated into a TFT exploration system and applied to actual data.
以下、本発明の切羽前方探査システムにおける前方速度推定方法、及びこれに用いる切羽前方探査における前方速度推定システムの実施の形態について、図面を参照して説明する。 Hereinafter, the forward speed estimation method in the face forward exploration system of the present invention and the embodiment of the forward speed estimation system in the face forward exploration used thereto will be described with reference to the drawings.
図１は、本発明の一実施の形態に係る切羽前方探査システムにおけるデータ収集部のハードウェア構成を表すブロック図である。図２は、上記実施の形態のデータ収集部で収集されたデータを受け取ってデータの処理及び解析を行うデータ解析部のハードウェア構成を表すブロック図である。図３は上記実施の形態において実行される切羽前方探査の処理手順を説明するフローチャートである。 FIG. 1 is a block diagram showing a hardware configuration of a data collection unit in the face forward exploration system according to the embodiment of the present invention. FIG. 2 is a block diagram showing a hardware configuration of a data analysis unit that receives data collected by the data collection unit of the above embodiment, processes and analyzes the data. FIG. 3 is a flowchart illustrating a processing procedure of the face forward exploration executed in the above embodiment.
図１及び図２において、切羽前方速度推定システムは、発破によって発生した弾性波信号を記録するデータ収集部Ａ（図１に示される）と、弾性波解析ユニットとしてのＰＣ（パーソナルコンピュータ）１０（図２に示される）とを備えて構成される。 In FIGS. 1 and 2, the face forward velocity estimation system includes a data collection unit A (shown in FIG. 1) that records elastic wave signals generated by bursting, and a PC (personal computer) 10 (personal computer) 10 as an elastic wave analysis unit. (Shown in FIG. 2).
データ収集部Ａは、トンネル坑内に設置され、発破器から発生する発破信号を計測する発破信号検知器１と、発破信号を出力する装置１−１とを内蔵するトリガーユニットＢと、発破を振源とする弾性波を受振し、受振し弾性波信号を記録する受振ユニットＣとから成る。受振ユニットＣは、発破を振源とする弾性波を受振、計測し弾性波信号を出力する地震計２と、発破信号及び地震計２の弾性波信号を記録する弾性波デジタル記録装置３と、弾性波デジタル記録装置３に記録された弾性波データを外部記録媒体へ出力し記録するためのデータ出力部４と、弾性波デジタル記録装置３に記録された弾性波データを上記ＰＣ１０、或いは他のユニット、例えばデータを格納するためのサーバへ送信するための通信部５とを備えて成る。
The data collection unit A is installed in the tunnel mine, and has a trigger unit B incorporating a
トリガーユニットＢは、発破を起爆する発破器の爆薬に点火電流を供給するための発破母線に設置され、発破点火時に発破器から出力される電流を検出し、電流検出信号を出力する非接触式の電流センサと、発破信号を出力する装置とを備えており、電流センサーは一般に市販される非接触式の直流電流センサが採用される。 The trigger unit B is installed on a blasting bus for supplying an ignition current to the explosive of the blasting device that detonates the blasting, and is a non-contact type that detects the current output from the blasting device at the time of blasting ignition and outputs a current detection signal. The current sensor is equipped with a device for outputting a blasting signal, and a non-contact type DC current sensor generally available on the market is adopted as the current sensor.
受振ユニットＣは、トンネル坑内の支保工ロックボルトの頭部に設置される。ロックボルトは通常、トンネルの壁面を３〜４ｍ程度削孔した孔に挿入されて、この孔にグラウト材を注入してトンネル坑壁の地山に一体化されている。ロックボルトの頭部には受振ユニットＣがクランプ、ねじなどにより圧接して構成されており、該受振ユニットＣ内に内蔵される形で地震計２が設置されている。なお、地震計２単体では指向性を持たないが、設置のし方により、例えば、少なくとも３個以上の地震計２を所定位置に配置することにより、反射面の走向・傾斜を推定可能にするものである。したがって、図１において、地震計２は１個が示されているが、複数個設置されていてもよい。また、本発明において、地震計２の型式は単成分測定式であっても良いし、３次元測定式であっても良い。
The vibration receiving unit C is installed at the head of the support lock bolt in the tunnel mine. The lock bolt is usually inserted into a hole in which the wall surface of the tunnel is drilled by about 3 to 4 m, and a grout material is injected into the hole to be integrated with the ground of the tunnel wall. A vibration receiving unit C is pressed against the head of the lock bolt by a clamp, a screw, or the like, and a
弾性波デジタル記録装置３は、受振ユニットＣに内蔵されており、地震計２と接続されており、発破信号検知器１から出力される発破信号と地震計２からそれぞれ出力される弾性波信号とを記録可能である。記録装置として、例えば、一般に市販される８チャンネルのマルチチャンネルレコーダが採用される。また、ハードディスクのような設置型記録装置や光ディスク、ＣＤ或いはＤＶＤのような外部取り出し可能な記録媒体を用いた記録装置が使用されてもよい。
The elastic wave
データ出力部４は、弾性波デジタル記録装置３に接続され、記録されたデータを外部記録媒体へ出力する。通信部５は、弾性波デジタル記録装置３に接続され、記録されたデータを他のユニット、例えば弾性波解析ユニットやデータを格納するためのサーバへ送信する。この場合、データ出力部４としては、ＳＤメモリカード他各種のメモリカード用又はＰＣカード用のカードスロット、ＵＳＢメモリ用の端子、各種のカードメモリリーダー／ライター、各種のディスクドライバーなどの各種のインタフェース機器が用いられる。また、通信部５は、弾性波デジタル記録装置３に記録されたデータを通信により送出したり外部からのデータ要求を受信したりするために、インターネットなどの公衆通信ネットワークに接続される。
The
弾性波解析ユニットとしてのＰＣ１０は、図２に示すように、弾性波解析のための各種指示が入力される操作入力部１１と、弾性波解析用の各種アプリケーション（ソフトウェア）及びデータが格納される記憶部１２と、弾性波解析処理のための情報或いは処理結果が表示される表示部１３と、外部記録媒体に格納された弾性波データ等を読み出すデータ読取部１４と、データ収集部Ａから送付されたデータを受信する通信部１５と、上記各機能部の動作をコントロールしまた各種演算処理を行う制御部１６とを備える。
As shown in FIG. 2, the
操作入力部１１は、キーボード、タッチパネル、音声入力マイクロホンなどのデータ入力機器により構成され、この操作入力部１１により、制御部１６における各種処理動作に必要なコマンドおよびデータが入力される。
The
記憶部１２は、少なくとも読み出し専用メモリ（ＲＯＭ）と、ランダムアクセスメモリ（ＲＡＭ）と、計測データメモリとを有しており、それぞれのメモリが必要に応じて使用される。ここで、読み出し専用メモリ（ＲＯＭ）は、弾性波デジタル記録装置３に切羽掘削毎に記録される弾性波から直達波を抽出し、当該直達波に基づいて切羽のトンネル坑内の弾性波速度を測定する手段としての直達波処理プログラム、及び弾性波デジタル記録装置３に切羽掘削毎に記録される弾性波から反射波を抽出し、当該反射波に基づいて弾性波の反射・回折点及び反射面の位置を推定する手段としての反射波・回折波処理プログラムを含む各種のプログラムが格納されている。またランダムアクセスメモリ（ＲＡＭ）に対しては、処理動作に際してデータの書き込み、読み出しが実行される。そして計測データメモリには、発破信号を計測して得られた発破計測データ及び弾性波を計測して得られた弾性波計測データが格納される。
The
表示部１３は液晶その他のディスプレイ機器からなり、弾性波の解析により得られた坑内弾性波速度及び切羽前方の弾性波反射面などこの弾性波探査システムの動作中における種々の状態情報や処理情報などが表示される。データ読取部１４としては各種のメモリカードスロット、各種のカードメモリリーダー／ライター、各種のディスクドライバーなど各種のインタフェース機器が用いられる。通信部１５はデータ収集部Ａから通信により送られてきた記録データを受信するためにインターネットなどの公衆通信ネットワークに接続される。制御部１６は、例えば、ＣＰＵなどからなり、３次元弾性波などのデータ解析処理を実行したり、解析結果をサーバへ送付したりする。
The
以上の構成を有する切羽前方探査システムにおいて、本実施の形態の前提となる弾性波探査方法、及び切羽前方速度推定方法を実現する動作について説明する。なお、弾性波探査方法（弾性波のデータを集める操作）に関する以下の説明は、複数の地震計２を用いて一般的に行われると想定される手法を説明するものであり、本実施形態に特有の手法ではないことに留意すべきである。ここで、図３は、本発明で用いる切羽前方速度推定方法の動作を説明するフローチャートである。図４は探査時の探査状況を模式的に例示する説明図である。
In the face forward exploration system having the above configuration, the operation of realizing the elastic wave exploration method and the face forward velocity estimation method which are the premise of the present embodiment will be described. The following description of the elastic wave exploration method (operation of collecting elastic wave data) describes a method generally assumed to be performed using a plurality of
上記動作においては、先ず、図３に示すように、トンネル弾性波探査方法の処理概要は、まずステップＳ１で、トンネル坑内の切羽から離れた後方所定の位置に発破信号検知器１及び多成分デジタル記録装置３を設置すると共に、各実施形態の態様に応じて複数個の単成分地震計２をトンネル坑内壁面に設置する。次にステップＳ２で、切羽に発破孔を設けて爆薬を装填する。そしてステップＳ３では、発破器から発生し爆薬を起爆させる発破信号を計測するとともに、爆薬を起爆させた際に切羽を振源として発生する弾性波を各単成分地震計２で受振、計測し、多成分デジタル記録装置３に記録することを切羽掘削毎に繰り返す。これにより、多地点起振の弾性波の波形が得られる。またステップＳ４では、切羽掘削毎に多成分デジタル記録装置３に記録された弾性波から直達波を抽出して、当該直達波に基づいて切羽のトンネル坑内の弾性波速度分布を推定する。さらにステップＳ５では、切羽掘削毎に多成分デジタル記録装置３に記録された弾性波から反射波を抽出して、当該反射波に基づいて反射・回折点及び反射面の３次元空間位置を推定する。
In the above operation, first, as shown in FIG. 3, the processing outline of the tunnel elastic wave exploration method is first described in step S1 by the
上記各ステップ動作のうちステップＳ３までの動作が切羽前方速度推定動作の中の「探査」動作に相当する。 Of the above step operations, the operation up to step S3 corresponds to the "exploration" operation in the face forward velocity estimation operation.
１．探査
まず、図４に示すように、掘削完了区間の切羽後方所定の位置（例えば、切羽から３０ｍ（メートル）後方）の坑壁壁面に、発破信号検知器１と発破信号出力部１−１とを内蔵したトリガーユニットＢ、及び弾性波デジタル記録装置３と地震計２とを内蔵した受振ユニットＣを設置する。なお、同図では地震計２が２台設置されているが、台数は限定されない。また、以下の説明では、Ｘ軸、Ｙ軸及びＺ軸により３次元空間を表すが、トンネル掘削方向（トンネル軸に平行な方向）をＸ軸方向とし、切羽に向かって左右方向で該Ｘ軸と直交する方向をＹ軸方向とし、また、切羽に向かって上下方向でＸ軸及びＹ軸とそれぞれ直交する方向をＺ軸方向とする。
1. 1. Exploration First, as shown in FIG. 4, the
また、ここで、切羽と受振ユニットＣとの間の距離は、測量機器などで測定され、距離データがＰＣ１０に入力される。発破は瞬発電気雷管及びＤＳ電気雷管を用いて起爆させることとし、切羽に発破孔を削孔し、電気雷管を装着した爆薬を装填する。そして、発破母線に発破信号検知器１を装着する。
Further, here, the distance between the face and the vibration receiving unit C is measured by a surveying instrument or the like, and the distance data is input to the
このようにして、切羽の掘削に当たって発破を行い、発破信号及び、地中の断層７等の脆弱部から反射されて届いた弾性波（反射波）信号を弾性波デジタル記録装置３に記録する。この場合、発破スイッチのＯＮ操作のみ行えばよく、この操作により、発破器から出力された電流が発破信号検知器１によって検出され、切羽に装填された爆薬は起爆し切羽が爆破される。この爆発によって切羽から発生した弾性波は地山を伝播して切羽後方の地震計２に到達し、弾性波が地震計２により受振、計測される。そして、発破信号検知器１から出力された発破信号と地震計２から出力された弾性波信号が弾性波デジタル記録装置３に自動的に記録される。かかるデータの収録後、弾性波デジタル記録装置３に記録されたデータは、ＳＤカードなどを介してＰＣ１０に入力され、データ解析処理が行われる。
In this way, blasting is performed when excavating the face, and the blasting signal and the elastic wave (reflected wave) signal reflected and received from a fragile portion such as a
このような切羽の掘削と共に行う発破信号の計測及び記録、弾性波の受振、計測及び記録を、弾性波の受振、計測点を変えないで、切羽掘削サイクル毎に繰り返す。これにより、多地点起振の波形が得られたことになる。そして、弾性波デジタル記録装置３に記録された各切羽の弾性波データは、通信手段、ＳＤカードなどを介して順次ＰＣ１０に入力され、ＰＣ１０内に格納された直達波処理プログラム及び反射波処理プログラムにより解析されて、その解析結果により、トンネル坑内の弾性波速度を推定し、切羽前方の反射・回折点を推定する。
The measurement and recording of the blasting signal and the reception, measurement and recording of the elastic wave performed together with the excavation of the face are repeated for each face excavation cycle without changing the reception of the elastic wave and the measurement point. As a result, a waveform of multipoint oscillation is obtained. Then, the elastic wave data of each face recorded in the elastic wave
図３に示された各ステップ動作のうちステップＳ４、５の動作が切羽前方速度推定動作の中の「解析」動作もしくは切羽前方速度推定動作に相当する。 Of the step movements shown in FIG. 3, the movements of steps S4 and S5 correspond to the "analysis" movement or the face forward speed estimation movement in the face forward speed estimation movement.
２．解析
この解析動作においては、上述のようにして多地点起振の波形が得られたならば、直達波と反射波に波動場を分離して反射波の発生地点を速度不連続面の位置とする。切羽から不連続面までの弾性波の進行速度は初動走時の傾きから得られる速度と同一とみなすことが可能であるが、速度不連続面の向こう側の速度が問題となる。
2. Analysis In this analysis operation, if the waveform of multipoint oscillation is obtained as described above, the wave field is separated into the direct wave and the reflected wave, and the generation point of the reflected wave is set as the position of the velocity discontinuity surface. To do. The traveling speed of elastic waves from the face to the discontinuous surface can be regarded as the same as the speed obtained from the inclination at the time of initial running, but the speed on the other side of the speed discontinuous surface becomes a problem.
ここでは、反射波の振幅情報から切羽前方の速度を推定することを考える。一般に、平面境界への垂直入射では反射波の振幅ＡRは以下の式で表される。
ここで、Ｚ1 は入射側のインピーダンスを示し、Ｚ2 は透過側のインピーダンスを示す。式（１）の右辺は、入射振幅Ａと反射係数との積である。式（１）からわかる通り、入射側のインピーダンスと透過側のインピーダンスとの間で、
Ｚ2 ＜Ｚ1
の場合（不連続面の向こう側が低インピーダンスである）、反射波形は正負が逆転する。
Here, it is considered to estimate the velocity in front of the face from the amplitude information of the reflected wave. Generally, at normal incidence to the plane boundary amplitude A R of the reflected wave is expressed by the following equation.
Here, Z 1 indicates the impedance on the incident side, and Z 2 indicates the impedance on the transmission side. The right side of equation (1) is the product of the incident amplitude A and the reflection coefficient. As can be seen from equation (1), between the impedance on the incident side and the impedance on the transmissive side,
Z 2 <Z 1
In the case of (the other side of the discontinuity surface has low impedance), the positive and negative of the reflected waveform are reversed.
図５は速度構造と反射係数、および想定される波形の関係を表す図である。ｉ層中の入射波振幅をＡT(i) 、反射波振幅をＡR(I) とすれば、Ｚ(I+1) は次の式で計算することができる。
すなわち、地山の密度は変化しないと仮定した場合、Ｖ１は初動走時の傾きで与えられるので、反射波と直達波の振幅の比が得られればＶ２は計算することができる。図６は速度推定の基本的動作の流れを表すフローチャートである。図６において、速度推定の動作は、まずステップＳ１１で前処理動作を行う。この前処理動作では、速度フィルターやバンドパスフィルターを用いて各種フィルター処理を行ったり、振幅復元処理を行う。次にステップＳ１２で、直達波と反射波を分離抽出処理する。この処理ステップでは、ＦＫフィルターや、τ-Pフィルターなどによるフィルター処理、および各種重合処理が含まれる。そしてステップＳ１３では、直達波振幅で反射波振幅を基準化処理する。この基準化処理により、起振力のバランスが補正される。ステップＳ１４では、反射係数列を計算処理し、またステップＳ１５では上記計算の結果得られた値を速度分布に変換処理する。上記反射係数列の計算処理及び速度分布への変換処理に関しては、コンボリューションモデル等のフォワードモデルを設定してインバーションによって求めることが多い。 That is, assuming that the density of the ground does not change, V1 is given by the slope at the time of the initial run, so V2 can be calculated if the ratio of the amplitudes of the reflected wave and the direct wave is obtained. FIG. 6 is a flowchart showing the flow of the basic operation of speed estimation. In FIG. 6, as for the speed estimation operation, the preprocessing operation is first performed in step S11. In this preprocessing operation, various filter processings are performed using a velocity filter and a bandpass filter, and amplitude restoration processing is performed. Next, in step S12, the direct wave and the reflected wave are separated and extracted. This treatment step includes a filter treatment with an FK filter, an τ-P filter, and various polymerization treatments. Then, in step S13, the reflected wave amplitude is standardized by the direct wave amplitude. By this standardization process, the balance of the vibrating force is corrected. In step S14, the reflection coefficient sequence is calculated, and in step S15, the value obtained as a result of the calculation is converted into a velocity distribution. Regarding the calculation process of the reflection coefficient sequence and the conversion process to the velocity distribution, a forward model such as a convolution model is often set and obtained by inversion.
単純化した速度推定方法
本発明が対象とする切羽前方速度推定システムを利用するユーザーには物理探査や地球物理学の知識を要求しない。したがって、解析パラメータを変えて試行錯誤的に解析を行うことは想定しない。また、現場に設置できる一般的なパーソナルコンピュータで１時間以内に結果を出すことが要件となっている。図７は切羽前方速度推定システムに組み込んだプログラムにより図６に示された動作を具体的に実行する切羽前方速度推定の処理動作を表すフローチャートである。
Simplified Velocity Estimating Method No knowledge of geophysics or geophysics is required for users who use the face forward velocity estimation system of the present invention. Therefore, it is not assumed that the analysis is performed by trial and error by changing the analysis parameters. In addition, it is a requirement that a general personal computer that can be installed in the field produces results within one hour. FIG. 7 is a flowchart showing a processing operation of face forward speed estimation that specifically executes the operation shown in FIG. 6 by a program incorporated in the face forward speed estimation system.
図７において、速度推定の動作は、まずステップＳ２１で初動読取りを行う（ステップＳ２１）。初動読取りにより得られたデータ（受振波形）は１／Ｒ補正・積分がさ変位波形に変換され、データの補正が行われる（ステップＳ２２）。このデータ補正処理以降の処理は、補正されたデータを０ウェイ変換するか、それとも２ウェイ変換するかにより処理の流れが異なる。 In FIG. 7, the speed estimation operation first performs initial reading in step S21 (step S21). The data (received waveform) obtained by the initial reading is converted into a 1 / R correction / integral displacement waveform, and the data is corrected (step S22). The flow of processing after this data correction processing differs depending on whether the corrected data is converted into 0-way or 2-way conversion.
ここでは、まず、０ウェイ変換する場合の処理の流れを説明する。「０ウェイ変換」処理では、上述のように受振波形を積分して変位波形に変換し、幾何減衰補正（ステップＳ２２）をした後に続く処理として、初動時刻をそろえるようにシフト（０ウェイ変換）して直達波の振幅の変位波形（データ）を得る（ステップＳ２３）。この変位波形データはスタックされ（ステップＳ２４）、次いでエンベロープ計算が行われる（ステップＳ２５）。さらに、エンベロープを基に初動付近の変位波形のピークが抽出され（ステップＳ２６）、０ウェイ変換された変位波形のピーク値で基準化が行われる（ステップＳ２７）。その後、ピーク付近のスタック波形の切り出しが行われる（ステップＳ２８）。次に変位波形の相関係数を計算し、相関係数の正、または負により位相を判断する（ステップＳ２９）。そして、弾性波の伝播媒体（地山）の密度が一定と仮定して反射係数から速度分布を計算して（ステップＳ３０）、一連の切羽前方速度推定の処理動作を終了する。 Here, first, the flow of processing in the case of 0-way conversion will be described. In the "0-way conversion" process, the vibration-receiving waveform is integrated and converted into a displacement waveform as described above, and as a process following the geometric attenuation correction (step S22), the shift (0-way conversion) is performed so that the initial motion times are aligned. Then, the displacement waveform (data) of the amplitude of the direct wave is obtained (step S23). The displacement waveform data are stacked (step S24), and then the envelope calculation is performed (step S25). Further, the peak of the displacement waveform near the initial motion is extracted based on the envelope (step S26), and the peak value of the displacement waveform converted to 0 way is used for standardization (step S27). After that, the stack waveform near the peak is cut out (step S28). Next, the correlation coefficient of the displacement waveform is calculated, and the phase is determined based on whether the correlation coefficient is positive or negative (step S29). Then, assuming that the density of the elastic wave propagation medium (ground) is constant, the velocity distribution is calculated from the reflection coefficient (step S30), and a series of processing operations for face forward velocity estimation is completed.
他方、上記ステップＳ２２に続く処理として、２ウェイ変換する場合の処理の流れを説明する。上述のように受振波形を積分して変位波形に変換し、幾何減衰補正（ステップＳ２２）をした後に続く処理として、２ウェイ変換が行われ、この変換により反射波の振幅の変位波形（データ）を得る（ステップＳ３１）。この変位波形データはフィルター処理が施されるとともにスタックされ（ステップＳ３２）、次いでエンベロープ計算が行われる（ステップＳ３３）。さらに、エンベロープを基に０ウェイ変換された変位波形のピーク値で基準化が行われ（ステップＳ３４）、次いで主要イベントの抽出が行われる（ステップＳ３５）。その後、ピーク付近のスタック波形の切り出しが行われる（ステップＳ３６）。その後は、上述の０ウェイ変換した場合の処理の流れと同様に、変位波形の相関係数を計算し、相関係数の正、または負により位相を判断する（ステップＳ２９）。そして、弾性波の伝播媒体（地山）の密度が一定と仮定して反射係数から速度分布を計算して（ステップＳ３０）、一連の切羽前方速度推定の処理動作を終了する。 On the other hand, as the process following step S22, the flow of the process in the case of 2-way conversion will be described. As described above, the vibration receiving waveform is integrated and converted into a displacement waveform, and as a subsequent process after the geometric attenuation correction (step S22), 2-way conversion is performed, and this conversion performs the displacement waveform (data) of the amplitude of the reflected wave. (Step S31). The displacement waveform data is filtered and stacked (step S32), and then envelope calculation is performed (step S33). Further, the peak value of the displacement waveform converted to 0 way based on the envelope is standardized (step S34), and then the main event is extracted (step S35). After that, the stack waveform near the peak is cut out (step S36). After that, the correlation coefficient of the displacement waveform is calculated and the phase is determined by the positive or negative of the correlation coefficient in the same manner as the processing flow in the case of the 0-way conversion described above (step S29). Then, assuming that the density of the elastic wave propagation medium (ground) is constant, the velocity distribution is calculated from the reflection coefficient (step S30), and a series of processing operations for face forward velocity estimation is completed.
上記のような速度推定の処理動作について、アルゴリズムには以下のような工夫を加えた。
（１）主要なイベントを抽出するためにエンベロープを使用した。
（２）反射位相の判定を容易にするため、変位波形を使用した。
（３）反射係数はエンベロープの比で算出した。
（４）反射位相の判定は波形の相関係数を利用した。
Regarding the processing operation of speed estimation as described above, the following device was added to the algorithm.
(1) An envelope was used to extract the major events.
(2) A displacement waveform was used to facilitate the determination of the reflection phase.
(3) The reflectance coefficient was calculated by the ratio of the envelopes.
(4) The correlation coefficient of the waveform was used to determine the reflection phase.
シミュレーション波形での推定結果
３次元粘性弾性差分法によって波形を計算し本手法を適用したモデルを図８及び図９に示す。モデルＡは単純な２層構造、モデルＢは低速度帯狭在モデルで、かまぼこ型断面のトンネル空洞をモデル化している。トンネル軸と速度境界面とは直交しており、受振点は、Ｘ＝１２０（ｍ）で、トンネル切羽に向かって右壁中央、坑壁杭から１（ｍ）奥とし、起振点はトンネル切羽中心の３（ｍ）奥で、切羽は、Ｘ＝１５０（ｍ）からＸ＝１９６（ｍ）まで、２（ｍ）毎、合計２４ケースでモデル化し、波形を計算した。震源関数は中心周波数２００Ｈｚのリッカーウェーブレットとした。
Estimated Results with Simulation Waveforms Figures 8 and 9 show models to which the waveform was calculated by the three-dimensional viscoelastic difference method and this method was applied. Model A is a simple two-layer structure, and model B is a low-speed band narrowing model, which models a tunnel cavity with a semi-cylindrical cross section. The tunnel axis and the velocity boundary surface are orthogonal, the vibration receiving point is X = 120 (m), the center of the right wall toward the tunnel face, 1 (m) behind the tunnel wall pile, and the vibration receiving point is the tunnel. At the depth of 3 (m) in the center of the face, the face was modeled from X = 150 (m) to X = 196 (m) every 2 (m) in a total of 24 cases, and the waveform was calculated. The epicenter function was a licker wavelet with a center frequency of 200 Hz.
上記モデルＡ、モデルＢに関し、図８は、３次元粘性弾性差分法によって波形を計算し、本発明の切羽前方速度推定手法を適用した場合の、反射波の速度が４ｋｍ／ｓである層と、２ｋｍ／ｓである層の２層構造のモデル（モデルＡ）を示す図である。また、図９は、３次元粘性弾性差分法によって波形を計算し、本発明の切羽前方速度推定手法を適用した場合の、反射波の速度が４ｋｍ／ｓである２つの層の間に２ｋｍ／ｓである層が狭在する構造のモデル（モデルＢ）を示す図である。 Regarding the above models A and B, FIG. 8 shows a layer in which the velocity of the reflected wave is 4 km / s when the waveform is calculated by the three-dimensional viscoelastic difference method and the face forward velocity estimation method of the present invention is applied. It is a figure which shows the model (model A) of the two-layer structure of the layer which is 2km / s. Further, FIG. 9 shows a waveform calculated by the three-dimensional viscoelastic finite difference method, and 2 km / s between two layers having a reflected wave velocity of 4 km / s when the face forward velocity estimation method of the present invention is applied. It is a figure which shows the model (model B) of the structure in which the layer which is s is narrowly squeezed.
図１０（ａ）にモデルＡで計算された速度振幅波形を示す。これを積分して変位波形に変換し幾何減衰補正をする。その後、初動時刻をそろえるようにシフトしたものを図１０（ｂ）に示し、２ウェイ変換したものを図１０（ｃ）に示す。なお図１０（ａ）、図１０（ｂ）、図１０（ｃ）において、横軸は時間を表し、その単位は秒である。 FIG. 10A shows the velocity amplitude waveform calculated by the model A. This is integrated and converted into a displacement waveform to correct the geometric attenuation. After that, the one shifted so that the initial motion times are aligned is shown in FIG. 10 (b), and the two-way conversion is shown in FIG. 10 (c). In FIGS. 10 (a), 10 (b), and 10 (c), the horizontal axis represents time, and the unit thereof is seconds.
図１１に図１０（ｂ）をスタックした波形と、そのエンベロープを示し、図１２に図１０（ｃ）をスタックした波形と、そのエンベロープを示す。図１１を直達波とみなし、図１２を反射波とみなす。 FIG. 11 shows a waveform in which FIG. 10 (b) is stacked and its envelope, and FIG. 12 shows a waveform in which FIG. 10 (c) is stacked and its envelope. FIG. 11 is regarded as a direct wave, and FIG. 12 is regarded as a reflected wave.
図１２の星印及び△印は、検出された反射波ピークの位置を示している。これらの印のうち、△印は同位相を表す一方、星印は逆位相を表す。図１２中の０．０２秒付近および０．０７秒付近の同位相のピークは直達波の除去が不完全なための誤認識である。 The asterisk and Δ in FIG. 12 indicate the positions of the detected reflected wave peaks. Of these marks, the Δ mark represents the same phase, while the star mark represents the opposite phase. The peaks of the same phase around 0.02 seconds and around 0.07 seconds in FIG. 12 are erroneous recognitions due to incomplete removal of the direct wave.
図１３に、計算された反射係数列と推定速度モデルを示す。誤認識があるものの、０．０４秒付近の速度境界は明瞭に捕らえられている。 FIG. 13 shows the calculated reflectance coefficient sequence and the estimated velocity model. Although there is a misunderstanding, the velocity boundary around 0.04 seconds is clearly captured.
図１４に、モデルＢの波形と反射波抽出の処理結果を示す。０．０５秒付近に、図１２では存在しなかった同位相の反射波ピークが検出されている。 FIG. 14 shows the waveform of model B and the processing result of the reflected wave extraction. Around 0.05 seconds, a reflected wave peak of the same phase, which did not exist in FIG. 12, was detected.
図１５に、モデルＢで反射係数列と推定速度モデルを示す。誤認識があるものの、０．０４秒付近の低速度帯は明瞭に捕らえられている。 FIG. 15 shows a reflectance coefficient sequence and an estimated velocity model in model B. Although there is a misunderstanding, the low speed band around 0.04 seconds is clearly captured.
本発明において考案されたアルゴリズムをＴＦＴ探査システムに組み込んで実測データに適用した。このときの波形例を図１６に示す。また上記実測データに適用した場合の直達波抽出および反射波抽出の処理結果を図１７に示す。上記アルゴリズムを適用したトンネルは上記実測データを取得したときは施工中であり、工事の進行に伴って反射波抽出の処理結果の妥当性の評価が検討される。 The algorithm devised in the present invention was incorporated into the TFT exploration system and applied to the measured data. An example of the waveform at this time is shown in FIG. Further, FIG. 17 shows the processing results of the direct wave extraction and the reflected wave extraction when applied to the measured data. The tunnel to which the above algorithm is applied is under construction when the above-mentioned actual measurement data is acquired, and the evaluation of the validity of the reflected wave extraction processing result is examined as the construction progresses.
切羽での掘削発破を起振源とするトンネルの切羽前方探査において、切羽前方の速度推定アルゴリズムを考案し、シミュレーション波形にて検討を行った結果、適用可能と判断できる結果を得た。また、切羽前方探査システムに組み込んで実測データへの適用を試みることが可能である。本手法は、施工管理上有用な情報を必要なタイミングで提供可能な方法である。 We devised a velocity estimation algorithm in front of the face in the forward exploration of the tunnel using excavation and blasting at the face as the trigger, and examined it using simulation waveforms. As a result, we obtained results that could be judged to be applicable. In addition, it is possible to incorporate it into the face forward exploration system and try to apply it to actual measurement data. This method is a method that can provide useful information for construction management at the required timing.
Ａ データ収集部
Ｂ トリガーユニット
Ｃ 受振ユニット
１０ ＰＣ（弾性波解析ユニット）
１ 発破信号検知器
１−１ 発破信号出力部
２ 地震計
３ 弾性波デジタル記録装置
４ データ出力部
５ 通信部
１１ 操作入力部
１２ 記憶部
１３ 表示部
１４ データ読取部
１５ 通信部
１６ 制御部
A Data collection unit B Trigger unit C
1 Break signal detector 1-1 Break
Claims (3)
切羽掘削毎に前記弾性波記録手段に記録された弾性波から直達波と反射波とを分離抽出する弾性波分離抽出ステップと、
分離抽出された弾性波についてエンベロープ計算を行うエンベロープ演算ステップと、
エンベロープ計算された弾性波についてピーク抽出を行うピーク抽出ステップと、
ピーク抽出された直達波のピーク値で反射波振幅を基準化する弾性波基準化ステップと、
反射波のピーク付近の波形を切り出す波形切出しステップと、
反射波の反射係数をエンベロープの比で算出する反射係数演算ステップと、
相関係数を計算する相関係数演算ステップと、
相関係数の正負により、反射波の位相を判定する位相を判定ステップと、
反射係数から速度分布を算出し、当該速度分布から弾性波の反射面の３次元位置を推定する反射面推定ステップと、
を有することを特徴とする切羽前方速度推定方法。 Blasting signal detecting means and elastic wave recording means are installed at a predetermined position behind the face, and one or more seismographs are installed on the inner wall surface of the tunnel to provide a blasting hole in the face inside the tunnel to explode. Is loaded, the blasting signal generated from the blasting device to detonate the explosive is measured, and the elastic waves generated from the face when the explosive is detonated are received and measured by the seismograph, respectively, and used as the elastic wave recording means. Recording is repeated for each face blasting, the reflected wave is extracted from the elastic wave recorded in the elastic wave recording means for each face blasting, and the three-dimensional position of the reflecting surface is estimated based on the reflected wave. It is a speed estimation method
An elastic wave separation / extraction step for separating and extracting direct waves and reflected waves from elastic waves recorded in the elastic wave recording means for each face excavation.
Envelope calculation step to calculate envelope for separated and extracted elastic waves,
A peak extraction step for peak extraction of envelope-calculated elastic waves,
An elastic wave standardization step that standardizes the reflected wave amplitude with the peak value of the peak-extracted direct wave,
A waveform cutting step that cuts out the waveform near the peak of the reflected wave, and
The reflectance coefficient calculation step that calculates the reflectance coefficient of the reflected wave by the ratio of the envelope,
Correlation coefficient calculation step to calculate the correlation coefficient and
The phase for determining the phase of the reflected wave is determined by the positive / negative of the correlation coefficient.
A reflection surface estimation step that calculates the velocity distribution from the reflection coefficient and estimates the three-dimensional position of the reflection surface of the elastic wave from the velocity distribution.
A method for estimating the forward velocity of the face, which comprises.
受振波形を積分して変位波形に変換し、幾何減衰補正するステップと、
初動時刻をそろえるようにシフト（０ウェイ変換）して直達波の振幅の変位波形を得るステップと、
得られた変位波形データをスタックするステップと、を含み、
その次段階でエンベロープ演算ステップを実行する、
ことを特徴とする請求項１に記載の切羽前方速度推定方法。 After performing the elastic wave separation extraction,
The step of integrating the received vibration waveform, converting it to a displacement waveform, and correcting the geometric damping,
A step to obtain a displacement waveform of the amplitude of the direct wave by shifting (0-way conversion) so that the initial motion time is aligned, and
Including a step of stacking the obtained displacement waveform data,
Perform the envelope calculation step in the next step,
The method for estimating the forward speed of the face according to claim 1.
受振波形を積分して変位波形に変換し、幾何減衰補正するステップと、
変位波形に対して２ウェイ変換を行い、この変換により反射波の振幅の変位波形を得るステップと、
この変位波形データにフィルター処理が施すとともにスタックするステップと、を含み、
その次段階でエンベロープ演算ステップを実行する、
ことを特徴とする請求項１に記載の切羽前方速度推定方法。 After performing the elastic wave separation extraction,
The step of integrating the received vibration waveform, converting it to a displacement waveform, and correcting the geometric damping,
The step of performing 2-way conversion on the displacement waveform and obtaining the displacement waveform of the amplitude of the reflected wave by this conversion,
This displacement waveform data includes steps to filter and stack.
Perform the envelope calculation step in the next step,
The method for estimating the forward speed of the face according to claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019069126A JP2020165925A (en) | 2019-03-29 | 2019-03-29 | Face forward velocity estimation method in face forward survey system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019069126A JP2020165925A (en) | 2019-03-29 | 2019-03-29 | Face forward velocity estimation method in face forward survey system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2020165925A true JP2020165925A (en) | 2020-10-08 |
Family
ID=72715189
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019069126A Pending JP2020165925A (en) | 2019-03-29 | 2019-03-29 | Face forward velocity estimation method in face forward survey system |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2020165925A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7470614B2 (en) | 2020-10-22 | 2024-04-18 | 株式会社安藤・間 | A method for estimating three-dimensional velocity structure ahead of the face in a face-ahead detection system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015158437A (en) * | 2014-02-25 | 2015-09-03 | 株式会社安藤・間 | Tunnel elastic wave survey method and tunnel elastic wave survey system used therefor |
| CN106932839A (en) * | 2015-12-31 | 2017-07-07 | 中国石油天然气股份有限公司 | It is a kind of to rock density and the method and device of rate curve adjustment |
-
2019
- 2019-03-29 JP JP2019069126A patent/JP2020165925A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2015158437A (en) * | 2014-02-25 | 2015-09-03 | 株式会社安藤・間 | Tunnel elastic wave survey method and tunnel elastic wave survey system used therefor |
| CN106932839A (en) * | 2015-12-31 | 2017-07-07 | 中国石油天然气股份有限公司 | It is a kind of to rock density and the method and device of rate curve adjustment |
Non-Patent Citations (1)
| Title |
|---|
| 山本浩之、他: "トンネル前方探査における速度推定手法の試み", トンネル工学報告集 第２８巻, vol. Ｉ−３０, JPN6021011425, 16 January 2019 (2019-01-16), pages 1 - 6, ISSN: 0004601091 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7470614B2 (en) | 2020-10-22 | 2024-04-18 | 株式会社安藤・間 | A method for estimating three-dimensional velocity structure ahead of the face in a face-ahead detection system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2750200C (en) | Processing time series data embedded in high noise | |
| RU2570221C2 (en) | Determining position of geologic horizon relative to manifestation of seismic pulse in seismic data | |
| JP6280768B2 (en) | Tunnel elastic wave exploration method and tunnel elastic wave exploration system used therefor | |
| US11112513B2 (en) | Method and device for estimating sonic slowness in a subterranean formation | |
| KR101547508B1 (en) | Apparatus for predicting front geological features and the method thereof | |
| US20150300161A1 (en) | Down Hole Subsurface Wave System with Drill String Wave Discrimination and Method of Using Same | |
| JP7190015B2 (en) | Face forward exploration method | |
| CN100456035C (en) | Blast source based dam nondestructive detecting system | |
| EP2960683A1 (en) | Calibration of pressure gradient recordings | |
| CN105308479A (en) | Multi-parameter inversion through offset dependent elastic fwi | |
| EA025309B1 (en) | Geological mapping method and system and machine-readable medium | |
| US7889596B2 (en) | Method and a system for determining the position of a drill bit | |
| EP0299862A2 (en) | Method for evaluating parameters related to the elastic properties of subsurface earth formations | |
| WO2021033503A1 (en) | Seismic observation device, seismic observation method, and recording medium in which seismic observation program is recorded | |
| US20150268365A1 (en) | Method to characterize geological formations using secondary source seismic data | |
| JP5757474B2 (en) | Fault survey apparatus and fault survey method | |
| JP2020165925A (en) | Face forward velocity estimation method in face forward survey system | |
| JP6764725B2 (en) | Tunnel elastic wave exploration method and tunnel elastic wave exploration system used for this | |
| CN110687607A (en) | Stoneley wave detection method and system | |
| WO2019224354A1 (en) | Method for establishing position and timing of seismic recording devices that are deployed on a seafloor | |
| RU2736810C1 (en) | Formation well evaluation system and method | |
| US20170336523A1 (en) | Seismic signal processing method, apparatus and system | |
| JP6111107B2 (en) | Elastic wave exploration method | |
| JP7470614B2 (en) | A method for estimating three-dimensional velocity structure ahead of the face in a face-ahead detection system | |
| JP4069177B2 (en) | Fault fracture zone determination method and determination support device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A80 | Written request to apply exceptions to lack of novelty of invention |
Free format text: JAPANESE INTERMEDIATE CODE: A80 Effective date: 20190422 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200708 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210330 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210531 |
|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20210610 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20210611 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210811 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20210531 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20210928 |