JP2019113514A - Tunnel cut-out probe method - Google Patents
Tunnel cut-out probe method Download PDFInfo
- Publication number
- JP2019113514A JP2019113514A JP2018009056A JP2018009056A JP2019113514A JP 2019113514 A JP2019113514 A JP 2019113514A JP 2018009056 A JP2018009056 A JP 2018009056A JP 2018009056 A JP2018009056 A JP 2018009056A JP 2019113514 A JP2019113514 A JP 2019113514A
- Authority
- JP
- Japan
- Prior art keywords
- vibration data
- time
- geophone
- drilling
- 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.)
- Granted
Links
Images
Landscapes
- Excavating Of Shafts Or Tunnels (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
本発明は、トンネル切羽前方探査方法に関し、特に、打撃振動により生じた弾性波を用いてトンネルの切羽前方地山の状態を探査するトンネル切羽前方探査方法に関するものである。 The present invention relates to a tunnel face forward search method, and more particularly, to a tunnel face forward search method for searching for a state of a front face of a tunnel face using elastic waves generated by impact vibration.
トンネル切羽前方地山の状態を的確に予測することは、安全で経済的に施工するために必要不可欠である。そして、切羽前方地山の状態を探査する技術として、切羽から水平ボーリングを行い、試料（コア）を採取して物性を観察する技術や、ボーリング孔内にカメラを挿入して、直接地山の状態を観察する技術などが用いられていた。しかし、これらの技術は、費用が高額であること、作業期間中は工事ができないことなどの問題がある。そこで、コアを採取しないノンコアボーリングによる探査が行われることが多くなった。 Accurately predicting the condition of the front of the tunnel face is essential for safe and economical construction. Then, as a technique for detecting the state of the face in front of the face, horizontal boring is carried out from the face, and a sample (core) is collected to observe physical properties, or a camera is inserted into the borehole, Techniques for observing the state have been used. However, these techniques have problems such as high cost and inability to perform construction during the work period. Therefore, exploration by non-core boring which does not collect the core is often performed.
ノンコアボーリングによる探査である削孔検層法では、インパルス振動が用いられる場合が多い。具体的には、台車の上に多数の掘削機が取り付けられた大型の油圧式掘削機で時間をおいて打撃を繰り返すなど、本来の掘削作業ではない計測用の振動を発生させることで探査を行うものである。そして、油圧式掘削機の油圧ドリフタに設置された加速度計などを用いて打撃時間を特定して、切羽表面に設置した受振器に到達した時間の遅れと波形から地山の状態を特定している。 In the drilling logging method which is the exploration by non-core boring, impulse vibration is often used. Specifically, a large hydraulic excavator with many excavators mounted on a bogie repeatedly strikes at time, etc. It is something to do. Then, the impact time is specified using an accelerometer or the like installed in the hydraulic drifter of the hydraulic excavator, and the state of the earth is specified from the delay in time and waveform that reached the geophone installed on the face surface. There is.
このノンコアボーリングによる探査は、山岳トンネルの施工に用いられる大型の油圧式掘削機を利用できること、削孔速度や掘削機の油圧データを用いて地山の硬軟の程度を定量的に評価できるなどの利点がある。 In this non-core boring survey, it is possible to use a large hydraulic excavator used for construction of mountain tunnels, and to quantitatively evaluate the degree of hardness of the ground using drilling speed and hydraulic data of the drilling machine. There is an advantage.
ここで、ノンコアボーリングによるトンネル切羽前方探査については、特許文献１（特許第２７１０７４０号明細書）および特許文献２（特許第４１５７６３５号明細書）に記載された技術が知られている。 Here, as for the tunnel face forward exploration by non-core boring, the technology described in patent documents 1 (patent 2710740 specification) and patent documents 2 (patent 4157635 specification) is known.
特許文献１には、切羽表面に複数の受振器を設置し、油圧ドリフタによって切羽前方を削孔する際に生じる弾性波を計測する技術が記載されている。当該技術では、油圧ドリフタに取り付けられた加速度計の打撃データと各受振器の計測データとを比較して、スペクトル解析とジオモトグラフィ解析を行い、切羽前方地山の地質断面や物性分布を比較している。
また、特許文献２には、トンネルボーリングマシン（ＴＢＭ）に複数の加速度計を設置し、ＴＢＭによって掘削する際に生じる弾性波がカッタ前方の弾性波反射面で反射して帰ってきた波のデータに相互相関関数を用いて地山の硬軟を解析する技術が記載されている。また、先行ボーリング掘削を行う際に、油圧ドリフタと切羽表面に加速度計を設け、削孔ビットが打撃した際に発生した弾性波から地質状況を判別する技術が記載されている。
Further, in
なお、特許文献３（特開２０１６−０１７９００号公報）には、油圧ドリフタと削孔ビットとをつなぐシャフト上に複数の受振器を設置して、その計測値の相互相関関数を求める技術が記載されている。 Patent Document 3 (Japanese Patent Laid-Open No. 2016-017900) describes a technique in which a plurality of geophones are installed on a shaft connecting a hydraulic drifter and a drilling bit, and the cross-correlation function of the measured values is obtained. It is done.
ここで、ノンコアボーリングの探査においては、打撃の開始時間は油圧ドリフタの振動を基準としているが、得られた計測データをそのまま他の計測データ（切羽表面の受振器の計測データなど）と比較することは好ましくないことが分かった。 Here, in the non-core boring survey, the impact start time is based on the vibration of the hydraulic drifter, but the obtained measurement data is compared with other measurement data (measurement data of a surface striker on the face surface, etc.) as it is Was found to be undesirable.
すなわち、油圧ドリフタは、その構造上、大きなノイズが発生しやすく、真の打撃時間が特定しにくいからである。具体的には、掘削機先端の削孔ビットが岩盤を打撃した振動よりも油圧ドリフタから発生するノイズの方が大きくなるために、加速度計による計測波の最初のピークはノイズであって、その前の小さい波が打撃振動であり、当該打撃振動を捉えにくいからである。 That is, the hydraulic drifter is likely to generate a large noise due to its structure, and it is difficult to identify the true impact time. Specifically, the first peak of the measurement wave by the accelerometer is noise because the noise generated from the hydraulic drift is larger than the vibration of the drilling bit at the tip of the drilling machine hitting the bedrock. It is because it is hard to catch the said striking vibration since the small wave in front is striking vibration.
しかしながら、前述した特許文献１および特許文献２に記載の技術では、何れも油圧ドリフタに取り付けられた加速度計の計測値について、打撃時間を特定するための補正等は行っていない。
However, in the techniques described in
また、特許文献３に記載の技術では、シャフトに受振器を設置することは容易ではないことから、実工事でシャフト上の受振器から計測データを得ることは困難である。また、孔内に破砕したズリが残っている場合、削孔ビットと岩の間に入るなどしてノイズが発生することがあり、正確な解析を行うためには切羽表面の加速度計によるデータも補正する必要がある。
Moreover, in the technique described in
本発明は、上述の技術的背景からなされたものであって、トンネルの切羽前方地山の探査精度を向上させることができるトンネル切羽前方探査方法を提供することを目的とする。 The present invention has been made from the above technical background, and it is an object of the present invention to provide a tunnel face front exploration method capable of improving the exploration accuracy of the tunnel front face mine.
上記課題を解決するため、請求項１に記載の本発明のトンネル切羽前方探査方法は、トンネル切羽前方の地山に対し、掘削機のロッドの先端に設けられた削孔ビットにより所定の削孔長おきに所定の時間間隔で複数回打撃して弾性波を発生させ、前記掘削機に設置された第１の受振器による各打撃で分割された前記弾性波の複数の振動データ、およびトンネル切羽に設置された複数の第２の受振器による各打撃ごとに分割された前記弾性波の複数の振動データを削孔長ごとに取得し、前記第２の受振器で取得された複数の前記振動データに対して、当該振動データの相互相関関数が最大となる前記振動データを基準にして初動時刻をシフトさせた後にスタッキング処理を行い、スタッキングされた振動データを削孔長ごとに求め、前記第１の受振器で取得された前記振動データに対して、前記第２の受振器の前記振動データに対応してシフトしてスタッキング処理を行い、スタッキングされた振動データを削孔長ごとに求め、予め求められた前記掘削機が発振する発振時間と前記第１の受振器におけるスタッキングされた前記振動データの初動時刻を用いて、前記削孔ビットが地山を打撃した時刻でのそれぞれの前記第２の受振器におけるスタッキングされた前記振動データから初動時刻を読み取って走時曲線を求めてトモグラフィ解析を行い弾性波速度分布を求める、ことを特徴とする。
In order to solve the above-mentioned problems, according to the tunnel face front search method of the present invention according to
請求項２に記載の発明のトンネル切羽前方探査方法は、上記請求項１に記載の発明において、何れか一つの前記第２の受振器を基準受振器とし、当該基準受振器で取得された振動データに対して、相互相関関数が最大となる前記振動データを基準にして初動時刻をシフトさせた後にスタッキング処理を行ってスタッキングされた振動データを求め、前記基準受振器以外の前記第２の受振器で取得された前記振動データに対しては、前記基準受振器の前記振動データに対応してシフトしてスタッキング処理を行い、スタッキングされた振動データを求める、ことを特徴とする。
In the tunnel face forward search method of the invention according to
請求項３に記載の発明のトンネル切羽前方探査方法は、上記請求項２に記載の発明において、前記基準受振器は、前記掘削機に最も近い位置に設置された前記第２の受振器である、ことを特徴とする。
In the tunnel face forward exploration method according to
請求項４に記載の発明のトンネル切羽前方探査方法は、上記請求項１〜３の何れか一項に記載の発明において、前記スタッキング処理では、前記相互相関関数が所定値を下回る前記振動データは除外する、ことを特徴とする。
In the tunnel face forward search method of the invention according to
請求項５に記載の発明のトンネル切羽前方探査方法は、上記請求項１〜４の何れか一項に記載の発明において、前記走時曲線を求める際には、前記振動データの初動時刻を局所定常ＡＲ（ａｕｔｏｒｅｇｒｅｓｓｉｖｅ：自己回帰）モデルを用いて自動的に読み取る、ことを特徴とする。 In the tunnel face forward search method according to the fifth aspect of the present invention, in the invention according to any one of the first to fourth aspects, when obtaining the travel time curve, the initial motion time of the vibration data is locally It is characterized in that it is automatically read using a stationary AR (autoregressive) model.
本発明によれば、予め掘削機の発振時間を求めておき、ロッドの先端に設けられた削孔ビットが地山を打撃した際の弾性波がそれぞれの第２の受振器に届くまでの時間を第１の受振器と第２の受振器との初動時刻と発振時間とから計測し、削孔ビットと第２の受振器との距離から弾性波の速度を求めているので、トンネルの切羽前方地山の探査精度を向上させることが可能になる。 According to the present invention, the oscillation time of the drilling machine is obtained in advance, and the time until the elastic wave when the drilling bit provided at the tip of the rod strikes the ground reaches the respective second geophones Since the velocity of the elastic wave is determined from the distance between the drilling bit and the second geophone, by measuring the initial motion time and oscillation time of the first geophone and the second geophone, the face of the tunnel It will be possible to improve the search accuracy of the foreland.
以下、本発明の一例としての実施の形態について、図面に基づいて詳細に説明する。なお、実施の形態を説明するための図面において、同一の構成要素には原則として同一の符号を付し、その繰り返しの説明は省略する。 Hereinafter, an embodiment as an example of the present invention will be described in detail based on the drawings. Note that, in the drawings for describing the embodiments, the same components are denoted by the same reference symbols in principle, and the repetitive description thereof will be omitted.
図１は本発明の一実施の形態であるトンネル切羽前方探査を行う装置のレイアウトを示す概念図、図２は本発明の一実施の形態であるトンネル切羽前方探査に用いられる油圧ドリフタによる地山打撃を示す概念図である。 FIG. 1 is a conceptual view showing a layout of an apparatus for performing tunnel face forward search according to an embodiment of the present invention, and FIG. 2 is a ground with a hydraulic drifter used for tunnel face forward search according to an embodiment of the present invention. It is a conceptual diagram which shows a hit.
図１に示すように、本実施の形態であるトンネル切羽前方探査には、地山Ｇを掘削するための掘削機１０と、掘削機１０に取り付けられたパイロットセンサである加速度計（第１の受振器）Ｐと、切羽Ｓに設置された複数（本実施の形態では７台）の速度計（第２の受振器）Ｒと、加速度計Ｐおよび速度計Ｒで計測された振動を記録する計測装置Ｍとが用いられる。
As shown in FIG. 1, in the tunnel face forward search according to the present embodiment, an
また、図１および図２に示すように、掘削機１０は、油圧ドリフタ１１と、スリーブ１２を介して油圧ドリフタ１１のシャンクロッド１３に取り付けられるとともに先端に削孔ビット１４が固定されたロッド１５とを備えている。さらに、油圧ドリフタ１１は、シャンクロッド１３を打撃するピストン１６と、シャンクロッド１３を回転させるロータ１７とを有している。
Further, as shown in FIGS. 1 and 2, the
このような掘削機１０により、油圧ドリフタ１１内のピストン１６が作動油によりシリンダ１８内を移動し、シャンクロッド１３を打撃する。シャンクロッド１３は、打撃による打撃力のほかロータ１７からの回転力と油圧式掘削機の推力（フィード圧）Ｆをスリーブ１２を介してロッド１５に伝達し、ロッド１５は打撃力と回転力を削孔ビット１４に伝達する。先端の削孔ビット１４が打撃力を直接地山Ｇに加えることにより、地山Ｇが破砕される。
With such an
本実施の形態のトンネル切羽前方探査では、ロッド１５の先端に設けられた削孔ビット１４が地山Ｇを打撃した際の弾性波が各速度計Ｒに届くまでの時間を計測し、削孔ビット１４と速度計Ｒとの距離から弾性波の速度を求めている。その際、削孔ビット１４は削孔に伴って前進し、これにより削孔ビット１４と速度計Ｒとを結ぶ線が辺がする。これにより、得られた弾性波の速度からトモグラフィ解析を行って途中の岩盤の物性を予測している。
In the tunnel face forward search of the present embodiment, the time taken for the elastic wave to reach each speedometer R when the
ここで、図３に本実施の形態によるトンネル切羽前方探査の探査原理を示す。 Here, FIG. 3 shows the search principle of the tunnel face forward search according to the present embodiment.
図３において、油圧ドリフタ１１内のピストン１６がシャンクロッド１３を打撃することにより発生した振動（弾性波）は、打撃位置Ａから油圧ドリフタ１１を介して時間ＴＡＢ後に加速度計Ｐの位置Ｂに到達し、同じ時間ＴＡＣでロッド１５上の位置Ｃまで到達する。その後、当該振動が時間（掘削機１０が発振する時間：発振時間）ＴＣＤ後に削孔ビット１４の先端Ｄに到達すると岩石を破砕して地山Ｇに振動を発生させ、その振動が時間（伝播時間）ＴＤＥ後に位置Ｅに到達して速度計Ｒで受信される。
In FIG. 3, the vibration (elastic wave) generated by the
したがって、速度計Ｒの初動時刻ＴＥ、加速度計Ｐの初動時刻ＴＢ、前述した発振時間ＴＣＤが分かれば、伝播時間ＴＤＥは、 Therefore, if the initial motion time T E of the speedometer R, the initial motion time T B of the accelerometer P, and the oscillation time T CD described above are known, the propagation time T DE is
ＴＤＥ＝（ＴＥ−ＴＢ）−ＴＣＤ
で求めることができる。
T DE = (T E −T B ) −T CD
You can ask for
以上のことから、本実施の形態のトンネル切羽前方探査は、発振時間の把握、削孔振動の計測、削孔振動の波形処理、トモグラフィ解析、弾性波速度分布の把握の流れで行われる。 From the above, the tunnel face forward search of this embodiment is performed by the flow of grasping oscillation time, measurement of drilling vibration, waveform processing of drilling vibration, tomographic analysis, grasping elastic wave velocity distribution.
すなわち、発振時間の把握として、削孔ビット１４が地山Ｇを打撃するまでの時間（発振時間）を求める。次に、削孔振動の計測として、地山Ｇを削孔する際に油圧ドリフタ１１内のピストン１６がシャンクロッド１３を打撃することにより発生する振動を加速度計Ｐにより計測すると同時に、シャンクロッド１３およびロッド１５を介して伝達されたピストン１６の打撃力で削孔ビット１４が地山Ｇを打撃する振動（弾性波）を複数の速度計Ｒにより計測する。次に、削孔振動の波形処理として、計測された振動データの初動位置を調整した後にスタッキング処理する。そして、発振時間と弾性波が速度計Ｒに到達するまでの時間（走時）と削孔ビット１４の位置データとを用い、スタッキングされた振動データから走時曲線を求めてトモグラフィ解析を行い、切羽前方地山の弾性波速度分布を把握する。
That is, as a grasp of the oscillation time, the time (oscillation time) until the
さて、前述のように、地山Ｇの弾性波速度を求めるために、削孔ビット１４が地山Ｇを打撃するまでの時間である発振時間を事前に把握する。ここで、回転するロッド１５や地山Ｇを打撃する削孔ビット１４に直接センサを取り付けて発振時間を計測することは困難である。そこで、本実施の形態では、ピストン１６がシャンクロッド１３を打撃した振動は、ロッド１５を介して削孔ビット１４に伝播すると同様にピストン１６やシャンクロッド１３を覆う油圧ドリフタ１１のフレームにも間接的に伝播することに着目し、油圧ドリフタ１１のフレームの振動から削孔ビット１４が地山Ｇを打撃する発振時間を求めている。
Now, as described above, in order to obtain the elastic wave velocity of the ground G, the oscillation time which is the time until the
つまり、油圧ドリフタ１１内のピストン１６がシャンクロッド１３を打撃することにより発生した振動（弾性波）は、打撃位置から油圧ドリフタ１１を介して所定時間後に加速度計Ｐに伝播し、同時刻にロッド１５上のある位置（図３の位置Ｃ）まで伝播する。その後、当該振動が削孔ビット１４に到達すると岩石を破砕して地山Ｇに振動を発生させ、その振動が速度計Ｒで受信される。
That is, the vibration (elastic wave) generated by the
そこで、図４に示すように、油圧ドリフタ１１と削孔ビット１４に受振器として加速度計Ｐ１、Ｐ２をそれぞれ取り付けてピストン１６の打撃による振動を発生させ、加速度計Ｐ１に弾性波が到達する時間と加速度計Ｐ２に弾性波が到達する時間との時間差を把握する。この時間差は前述した発振時間になる。
Therefore, as shown in FIG. 4, accelerometers P1 and P2 are attached to the
例えば、油圧ドリフタ１１内の打撃位置から油圧ドリフタ１１の後部に設けられた加速度計１までの距離を０．６ｍ、当該打撃位置から削孔ビット１４に設けられた加速度計２までの距離を４．１ｍ（シャンクロッド１３の長さ０．４４ｍ＋ロッド１５の長さ３．６６ｍ）としたとき、削孔ビット１４を地山Ｇに押し当てて打撃したときの加速度計Ｐ１と加速度計Ｐ２とで計測された振動波形の一例を図５（ａ）および図５（ｂ）に示す。
For example, the distance from the striking position in the
図示するように、長さ３．６６ｍのロッド１５を使用した場合は、０．７０９ｍｓの時間差（つまり発振時間）になっているので、加速度計Ｐ１の振動波形の初動時刻が分かれば削孔ビット１４の打撃時刻を求めることができる。
As shown in the figure, when using a
なお、２本目以降のロッド１５を継ぎ足す場合には、同様にして、ロッド１５の一本あたりの弾性波の伝搬時間を求める。したがって、ロッド１５を継ぎ足した場合には、継ぎ足したロッドの本数に対応した伝播時間だけ、発振時間が長くなる。例えば継ぎ足すロッド１５の長さを３．０５ｍ、伝搬時間が０．６１１ｍｓとした場合、ロッド１５を２本継ぎ足すと、発振時間は、０．７０９ｍｓ＋（０．６１１ｍｓ×２）＝１．９３１ｍｓとなる。
When the second and
以上のようにして発振時間を求めたならば，次に、削孔振動の計測を行う。すなわち、油圧ドリフタ１１内のピストン１６によりシャンクロッド１３を所定の時間間隔で複数回打撃し、トンネルの切羽前方の地山Ｇに対して削孔ビット１４により所定の削孔長おきに弾性波を発生させる。そして、掘削機１０に設置された加速度計Ｐで各打撃ごとに分割された弾性波の振動データを取得するとともに、トンネル切羽Ｓに設置された複数の速度計Ｒで各打撃ごとに分割された弾性波の振動データを削孔長おきに取得する。
After the oscillation time is determined as described above, next, the measurement of drilling vibration is performed. That is, the
なお、削孔ビット１４による打撃は、例えば１秒間に５０〜６０回に設定される。但し、岩盤の硬軟などの要因によって、打撃回数や打撃間隔は各打撃で同一とはならない。よって、所定の時間間隔とは、一定の時間間隔という意味を有するものではない。
The impact by the
加速度計Ｐで取得された１０回の打撃による振動データの振動波形の一例を図６に、当該振動データの初動付近の振動波形の一例を図７に示す。また、速度計Ｒで取得された１０回の打撃による振動データの振動波形の一例を図８に、当該振動データの初動付近の振動波形の一例を図９に示す。なお、本実施の形態では、各速度計Ｒにおける各削孔長で１０回の打撃による１０種類の振動データを取得しているが、取得する振動データの数は複数であればよく、取得数は限定されない。 FIG. 6 shows an example of the vibration waveform of the vibration data of the ten impacts obtained by the accelerometer P, and FIG. 7 shows an example of the vibration waveform near the initial movement of the vibration data. Further, FIG. 8 shows an example of a vibration waveform of vibration data of ten impacts obtained by the speedometer R, and FIG. 9 shows an example of a vibration waveform near the initial movement of the vibration data. In the present embodiment, although 10 types of vibration data are obtained by 10 impacts with each drilling length in each speedometer R, the number of acquired vibration data may be plural, and the number of acquired Is not limited.
図７に示すように、ピストン１６がシャンクロッド１３を打撃することにより発生した振動が掘削機１０のフレームを伝播したものであることから、波形がパルス状にならないため、波の立ち上がりを読み取る初動に最大０．３２ｍｓのばらつきが生じている。また、図８に示すように、１回目と３回目の打撃では、削孔ビット１４が地山Ｇを打撃する際に破砕したずりが削孔ビット１４と地山Ｇとの間に挟まったことが原因と思われるノイズ（異常な部分）が発生している。さらに、図９に示すように、前述した加速度計Ｐの初動誤差により、受振波形の初動位置に最大０．３ｍｓのばらつきが生じている。
As shown in FIG. 7, since the vibration generated by the
さて、これら複数回（本実施の形態では１０回）の打撃による振動データの振動波形をそのままでスタッキング処理（波形データを時系列で足し合わせる処理）すると、Ｓ／Ｎ比が低下して初動時刻の読み取りに誤差を生じる。そこで、本実施の形態では、各振動波形間の相互相関関数を求め、相互相関関数が最大となる振動データを基準にして初動時刻を調整した後にスタッキング処理を行い、スタッキングされた振動データを求めている。スタッキング処理によりランダムノイズが打ち消し合って小さくなるために、振動データの振動波形のＳ／Ｎ比が向上する。その際、相互相関関数が所定値（本実施の形態では０．９以下）を下回る振動データは除外している。 Now, if the vibration waveform of the vibration data due to the plurality of impacts (10 times in this embodiment) is stacked as it is (processing to add waveform data in time series), the S / N ratio decreases and the initial motion time Error in the reading of Therefore, in the present embodiment, the cross correlation function between the vibration waveforms is determined, and the initial movement time is adjusted based on the vibration data for which the cross correlation function is maximum, the stacking process is performed, and the stacked vibration data is determined. ing. Since the random noise is canceled out and reduced by the stacking process, the S / N ratio of the vibration waveform of the vibration data is improved. At that time, vibration data whose cross correlation function is less than a predetermined value (in the present embodiment, 0.9 or less) is excluded.
相互相関関数が最大となる振動データを求める際に使用する速度計である基準速度計（基準受振器）Ｒは任意の一つでよい。この場合、掘削機１０に最も近い位置に設置された速度計を基準速度計Ｒとすることが望ましい。これは、削孔ビット１４からの振動は距離の２乗に反比例して小さくなるので、遠い速度計Ｒでは伝達した弾性波とノイズとの区別がしにくくなるからである。基準速度計Ｒによる以上の処理は、削孔長ごとの振動データについて行う。
The reference speedometer (reference geophone) R, which is a speedometer used when obtaining vibration data for maximizing the cross correlation function, may be any one. In this case, it is desirable that the speedometer installed at a position closest to the
なお、必ずしも相互相関関数が所定値を下回る振動データを除外する必要はなく、除外する場合でも、その閾値は自由に設定することができる。 In addition, it is not necessary to necessarily exclude the vibration data whose cross correlation function falls below a predetermined value, and even in the case of exclusion, the threshold can be set freely.
ここで、相互相関関数は、数１〜数３から計算することができる。
なお、ｙ（ｉ）はｉ打撃による振動データの振幅値、Ｎはデータ数、μ（ｉ）はｉ打撃による振動データの平均値、ｋはラグ（ｉ打撃とｊ打撃の振動データ間のシフト時間）、Ｃ（ｉ，ｊ）はｉ打撃とｊ打撃による振動データの相互共分散関数、Ｒ（ｉ，ｊ）はｉ打撃とｊ打撃による振動データの相互相関関数である。 Note that y (i) is the amplitude value of vibration data due to i impact, N is the number of data, μ (i) is the average value of vibration data due to i impact, k is lag (shift between vibration data of i impact and j impact Time), C (i, j) is a mutual covariance function of vibration data of i and j hits, and R (i, j) is a cross correlation function of vibration data of i and j hits.
打撃１０回分の振動データの振動波形に対して求めた相互相関関数およびその最大値の一例を図１０に、ｉ打撃が１から１０の相関関数のうち平均値が最も大きくなるラグ（ここではｉ＝８）を基準にして初動時刻をシフトさせた振動データの波形およびこれらの波形をスタッキング処理して求められた加速データの波形の一例を図１１に示す。 An example of the cross-correlation function and its maximum value obtained for the vibration waveform of the vibration data of 10 hits is shown in FIG. 10, and i is the lag in which the average value is the largest among the correlation functions 1 to 10 (here i FIG. 11 shows an example of the waveform of vibration data whose initial movement time is shifted with reference to 8) and the waveform of acceleration data obtained by stacking these waveforms.
なお、時刻をシフトさせるとは、上記のような方法で選定した振動データを基準にしてその他の振動データの時間軸を移動させ、初動時刻を合わせることをいう。 Note that shifting the time means moving the time axis of the other vibration data with reference to the vibration data selected by the above method, and adjusting the initial movement time.
図１０において、１回目と３回目の打撃による波形の相関値が０．７８、０．８８で０．９より小さくなっている。これは、前述の図８において、１回目の打撃による波形には１０ｍｓ付近、３回目の打撃による波形には７ｍｓ付近において他の打撃の波形に見られない振動を観測していることから、この影響を受けて相関値が低下していると推察される。したがって、０．９より小さくなる１回目と３回目の打撃による波形を除外してスタッキング処理をしている。 In FIG. 10, the correlation values of the waveforms of the first and third impacts are smaller than 0.9 at 0.78 and 0.88. This is because, in the above-mentioned FIG. 8, vibrations not observed in the waveforms of other strikes are observed in the vicinity of 10 ms in the waveform of the first impact and in the vicinity of 7 ms in the waveform of the third impact. It is inferred that the correlation value is lowered due to the influence. Therefore, the stacking process is performed excluding the waveforms of the first and third hits smaller than 0.9.
このようにして、基準速度計Ｒについてスタッキングされた振動データを求めたならば、基準速度計Ｒ以外の速度計Ｒに受信された各打撃に対応する振動データに対して、基準速度計Ｒの振動データに対応してシフトしてスタッキング処理を行い、スタッキングされた振動データを削孔長ごとに求める。 In this way, when the stacked vibration data is determined for the reference speedometer R, the vibration data corresponding to each impact received by the speedometer R other than the reference speedometer R Shifting and stacking are performed according to the vibration data, and the stacked vibration data is obtained for each drilling length.
続いて、加速度計Ｐで取得された１０回の打撃による振動データについてもスタッキング処理を行い、スタッキングされた加速データを削孔長ごとに求める。このとき、前述した基準速度計Ｒの振動データに対応してシフトしてスタッキング処理を行う。 Subsequently, stacking processing is also performed on vibration data of ten impacts acquired by the accelerometer P, and acceleration data thus stacked is determined for each drilling length. At this time, the stacking process is performed by shifting in accordance with the vibration data of the reference speedometer R described above.
このように、加速度計Ｐや基準速度計Ｒ以外の速度計Ｒで取得された振動データを基準速度計Ｒの振動データに対応してシフトさせるのは、次のような理由による。 The reason for shifting the vibration data acquired by the speedometer R other than the accelerometer P and the reference speedometer R according to the vibration data of the reference speedometer R in this way is as follows.
すなわち、加速度計Ｐで取得された各打撃の振動データと速度計Ｒ（基準速度計Ｒを含んだ全ての速度計Ｒ）で取得された各打撃の振動データとは、各回の打撃を検出するための切換タイミングの同期がとられている。具体的には、図７や図９などにおいて、打撃番号１（１回目の打撃）の振動データを検出しているときに、どのタイミングで打撃番号２（２回目の打撃）の振動データに切り換えるかという切換タイミングの同期がとられており、例えば、１５ｍｓ間隔で打撃番号が順次切り替えられるようになっている。したがって、加速度計Ｐや基準速度計Ｒ以外の速度計Ｒで取得された各打撃の振動データにおける初動の読み取り誤差と基準速度計Ｒで取得された各打撃の振動データにおける初動の読み取り誤差は同じになる。 That is, the vibration data of each impact acquired by the accelerometer P and the vibration data of each impact acquired by the speedometer R (all speedometers R including the reference speedometer R) detect each impact. The timing of switching for this is synchronized. Specifically, when vibration data of impact number 1 (first impact) is detected in FIGS. 7 and 9, etc., the vibration data of impact number 2 (second impact) is switched at any timing. The switching timing of the key is synchronized, and for example, the impact numbers are sequentially switched at intervals of 15 ms. Therefore, the reading error of the initial movement in the vibration data of each impact acquired by the speedometer R other than the accelerometer P and the reference speedometer R is the same as the reading error of the initial movement in the vibration data of each impact acquired by the reference speedometer R become.
そこで、基準速度計Ｒの各振動データと同じ時間だけ加速度計Ｐや基準速度計Ｒ以外の速度計Ｒで取得された各振動データをシフトさせれば、同一条件で初動時刻をシフトさせたことになるからである。 Therefore, if each vibration data acquired by the speedometer R other than the accelerometer P and the reference speedometer R is shifted by the same time as each vibration data of the reference speedometer R, the initial movement time is shifted under the same condition It is because
なお、振動データの切換タイミングについては、前述のように予め設定された時間間隔に基づいた方法に限定されるのではなく、これ以外の方法であってもよい。例えば、加速度計Ｐの振動をトリガとして切り換えるようにしてもよい。具体的には、加速度計Ｐの最大振幅の例えば２０％を閾値に設定しておき、これを超える時刻の例えば２ｍｓ前を切換タイミングとし、そこから所定時間の振動データを計測するようにしてもよい。 In addition, about the switching timing of vibration data, it is not limited to the method based on the time interval preset as mentioned above, You may be methods other than this. For example, the vibration of the accelerometer P may be switched as a trigger. Specifically, for example, 20% of the maximum amplitude of the accelerometer P is set as a threshold value, and for example, 2 ms before the time exceeding this is used as the switching timing, from which vibration data of a predetermined time is measured. Good.
加速度計Ｐで取得された振動データを基準速度計Ｒで取得された振動データに対応してシフトさせた振動データの波形およびこれらの波形をスタッキング処理して求められた加速データの波形の一例を図１２に示す。 Waveforms of vibration data obtained by shifting vibration data acquired by the accelerometer P in accordance with vibration data acquired by the reference speedometer R and an example of a waveform of acceleration data obtained by stacking these waveforms It is shown in FIG.
さて、このようにして各速度計Ｒでの削孔長おきのスタッキングされた振動データ、および加速度計Ｐでのスタッキングされた振動データが求められたならば、前述のようにして予め求められた発振時間（本実施の形態では、０．７０９ｍｓ＋継ぎ足したロッド１５の一本あたり０．６１１ｍｓの時間）と、加速度計Ｐにおけるスタッキングされた振動データの初動時刻を用いて、削孔ビット１４が地山Ｇを打撃した時刻での速度計Ｒにおけるスタッキングされた振動データ（各速度計Ｒでの所定の削孔長ごとのスタッキングされた振動データ）から初動時刻を読み取って走時曲線を求め、トモグラフィ解析を行い弾性波速度分布を求める。
Now, if the stacked vibration data for every drilling length in each speedometer R and the stacked vibration data for the accelerometer P are thus obtained, they are obtained in advance as described above. Using the oscillation time (in this embodiment, 0.709 ms + time of 0.611 ms for each added rod 15) and the initial movement time of the stacked vibration data in the accelerometer P, the
ここで、本実施の形態では、走時曲線を求める際、スタッキングされた振動データの波形に対し、ＡＩＣ（Ａｋａｉｋｅ’ｓ Ｉｎｆｏｒｍａｔｉｏｎ Ｃｒｉｔｅｒｉｏｎ：赤池情報量規準）が最も小さくなる時刻を初動時刻として、当該初動時刻を局所定常ＡＲ（ａｕｔｏｒｅｇｒｅｓｓｉｖｅ：自己回帰）モデルを用いて自動的に読み取る。 Here, in the present embodiment, when obtaining the travel time curve, the time at which AIC (Akaike's Information Criterion: Akaike's information criterion) becomes the smallest with respect to the waveform of the stacked vibration data is taken as the initial action time. The initial motion time is automatically read using a local stationary AR (autoregressive) model.
スタッキングされた加速度計Ｐの振動データの波形および当該波形から求めた初動位置の一例を図１３（ａ）（ｂ）に示す。また、スタッキングされた速度計Ｒの振動データの波形および当該波形から求めた初動位置の一例を図１４（ａ）（ｂ）に示す。 FIGS. 13A and 13B show an example of the waveform of vibration data of the stacked accelerometers P and an initial movement position obtained from the waveform. Moreover, an example of the waveform of vibration data of the stacked speedometer R and an initial movement position obtained from the waveform is shown in FIGS. 14 (a) and 14 (b).
これらの図面から、加速度計Ｐの初動時刻は２．０ｍｓ、速度計Ｒの初動時刻は５．２ｍｓであることから、両者の初動の時間差は３．２ｍｓと求められる。また、掘削機１０における弾性波の伝播時間（ロッド１５を継ぎ足さない場合）は０．７０９ｍｓであるから、削孔ビット１４から速度計Ｒまでの弾性波の伝播時間は２．４９１ｍｓとなる。
From these figures, since the initial motion time of the accelerometer P is 2.0 ms and the initial motion time of the speedometer R is 5.2 ms, the time difference between the two initial motions is determined to be 3.2 ms. Further, since the propagation time of the elastic wave in the excavator 10 (in the case where the
なお、初動時刻はＡＩＣ以外の手法によってもよく、また初動時刻は局所定常ＡＲモデルを用いた自動読み取り以外で求めてもよい。 The initial movement time may be determined by a method other than AIC, and the initial movement time may be determined by means other than automatic reading using a local steady state AR model.
加速度計Ｐおよび速度計Ｒで取得された弾性波の振動データをスタッキング処理し、削孔ビット１４が地山Ｇを打撃した時刻を基準にして０．２ｍの削孔間隔で波形を整理したグラフの一例を図１５および図１６に示す。図１５は掘削機１０に最も近い位置に設置された速度計Ｒ、図１６は掘削機１０から最も遠い位置に設置された速度計Ｒでのグラフである。
A graph in which the vibration data of elastic waves acquired by the accelerometer P and the speedometer R is subjected to stacking processing, and the waveforms are arranged at a drilling interval of 0.2 m based on the time when the
これらの図面において、波形ごとに初動の位置（点）をピッキングし、この点を通る線を引くことにより求まる走時曲線（削孔ビット１４の打撃により発生した弾性波が速度計Ｒに到達した時間と削孔長との関係を表す曲線）を示している。削孔長０．２ｍにおける走時は、図１５のグラフで０．６３０ｍｓ、図１６のグラフで２．３６４ｍｓであり、削孔位置から離れた速度計Ｒで得られるデータほど走時が大きくなっている。
In these drawings, the position (point) of the initial movement is picked for each waveform, and a travel time curve (elastic wave generated by the impact of the
走時曲線が求められたならば、トモグラフィ解析を行って切羽前方地山Ｇの弾性波速度分布を求める。トモグラフィ解析で求められたトンネル切羽前方の弾性波速度分布の一例を図１７に示す。トモグラフィ解析では、発振点（削孔位置）から受振点（速度計Ｒの設置位置）に至る波線の通過するセルに対して弾性波速度が反復修正される。そのため、発振点と受振点で囲まれる台形領域の底辺を三角形とする切羽前方の領域には波線が通過しないため解析対象外となる。図１７は、解析対象外部分を除いた波線が通過するセルの弾性波速度分布の結果を示している。図１７により、切羽付近では弾性波速度が３．５ｋｍ／ｓ〜４．５ｋｍ／ｓの地山Ｇが分布しているが、切羽前方の４ｍ付近から５ｋｍ／ｓ程度の速度の速い硬質な地山Ｇが左から右に徐々に拡大して現れること、その地山Ｇが切羽前方に約５ｍの厚さで断面方向から２０°の傾きで層状に存在することが予想できる。 After the travel time curve is determined, tomographic analysis is performed to determine the elastic wave velocity distribution of the face G in front of the face. An example of the elastic wave velocity distribution in front of the tunnel face obtained by tomographic analysis is shown in FIG. In tomographic analysis, the elastic wave velocity is iteratively corrected with respect to a cell through which a broken line from the oscillation point (the drilling position) to the receiving point (the installation position of the speedometer R) passes. Therefore, the broken line does not pass through the area in front of the face where the base of the trapezoidal area surrounded by the oscillation point and the receiving point is a triangle, so that it is not analyzed. FIG. 17 shows the result of the elastic wave velocity distribution of the cell through which the dashed line excluding the non-analysis part passes. According to Fig. 17, a ground G with an elastic wave velocity of 3.5km / s to 4.5km / s is distributed near the face, but a hard ground with a high speed of about 5km / s from around 4m in front of the face It can be expected that the mountain G appears gradually enlarged from the left to the right, and that the ground mountain G exists in a layer with a thickness of about 5 m in front of the face and a slope of 20 ° from the cross sectional direction.
このように、本実施の形態によれば、予め掘削機１０の発振時間を求めておき、ロッド１５の先端に設けられた削孔ビット１４が地山Ｇを打撃した際の弾性波が各速度計Ｒに届くまでの時間を加速度計Ｐと速度計Ｒとの初動時刻と発振時間とから計測し、削孔ビット１４と速度計Ｒとの距離から弾性波の速度を求めているので、トンネルの切羽前方地山の探査精度を向上させることが可能になる。
As described above, according to the present embodiment, the oscillation time of the
以上本発明者によってなされた発明を実施の形態に基づき具体的に説明したが、本明細書で開示された実施の形態はすべての点で例示であって、開示された技術に限定されるものではない。すなわち、本発明の技術的な範囲は、前記の実施の形態における説明に基づいて制限的に解釈されるものでなく、あくまでも特許請求の範囲の記載に従って解釈されるべきであり、特許請求の範囲の記載技術と均等な技術および特許請求の範囲の要旨を逸脱しない限りにおけるすべての変更が含まれる。 Although the invention made by the inventor has been specifically described based on the embodiments, the embodiments disclosed in the present specification are illustrative in all points and limited to the disclosed technology. is not. That is, the technical scope of the present invention is not to be interpreted restrictively based on the description in the above embodiment, and should be interpreted in accordance with the description of the claims only, and the claims It is intended to include all modifications without departing from the scope of the claims and the technical equivalents of the recited technology.
たとえば、本実施の形態では、何れか一つの速度計Ｒを基準速度計Ｒとしておき、この基準速度計Ｒで取得された振動データに対して、相互相関関数が最大となる振動データでスタッキング処理を行ってスタッキングされた振動データを求め、それ以外の速度計Ｒで取得された振動データに対しては、基準速度計Ｒの振動データに対応してシフトしてスタッキング処理を行い、スタッキングされた振動データを求めている。しかしながら、基準速度計Ｒを設定することなく、全ての速度計Ｒで取得された振動データに対して、相互相関関数が最大となる振動データでスタッキング処理を行ってスタッキングされた振動データを求めるようにしてもよい。 For example, in the present embodiment, any one speedometer R is used as a reference speedometer R, and the vibration processing data acquired by the reference speedometer R is subjected to stacking processing using vibration data that maximizes the cross correlation function. To determine the stacked vibration data, and for the other vibration data acquired by the speedometer R, the stacking process is performed by shifting in accordance with the vibration data of the reference speedometer R and stacking Vibration data is being sought. However, without setting the reference speedometer R, the vibration data acquired by all the speedometers R is subjected to the stacking process with the vibration data with the largest cross correlation function to obtain the stacked vibration data You may
また、本実施の形態では、掘削機１０には第１の受振器として加速度計Ｐを、切羽Ｓには第２の受振器として速度計Ｒを設置しているが、削孔ビット１４の打撃による振動の初動時刻を検出することができる限り、第１の受振器および第２の受振器としては、様々な検出手段を用いることができる。したがって、例えば、第１の受振器として速度計を用い、第２の受振器として加速度計を用いてもよい。
Further, in the present embodiment, although the accelerometer P is installed as a first geophone in the
さらに、弾性波速度分布から岩盤の性状を検討する際に、削孔時に破砕したずりや打撃時の反発力などを参考にしてもよい。 Furthermore, when examining the properties of the rock from the elastic wave velocity distribution, it is possible to refer to the shattering force at the time of drilling, the repulsive force at the time of impact, and the like.
さらに、本実施の形態では、削孔時の連続した打撃データを使用したが、例えばあるピッチだけ削孔するごとに一旦停止し、所定回数だけ連続して打撃したデータを用いるようにしてもよい。 Furthermore, in the present embodiment, continuous impact data at the time of drilling is used, but for example, it may be stopped once every drilling by a certain pitch, and data may be used continuously struck a predetermined number of times. .
そして、本実施の形態における説明で用いた様々の数値や速度計Ｒの設置数、グラフや表などは何れも一例に過ぎず、本発明がこれらによって限定されるものではない。 The various numerical values used in the description of the present embodiment, the number of installed speedometers R, graphs, tables, and the like are merely examples, and the present invention is not limited by these.
以上のように、本発明に係るトンネル切羽前方探査方法は、地山に対してトンネルを掘削する際に、打撃振動により生じた弾性波を用いてトンネル切羽前方の地山を探査するトンネル切羽前方探査方法に適用して有効である。 As described above, according to the tunnel face forward search method according to the present invention, when excavating the tunnel to the ground, the tunnel face forward to search the ground in front of the tunnel face using elastic waves generated by the impact vibration. It is effective when applied to the exploration method.
１０ 掘削機
１１ 油圧ドリフタ
１２ スリーブ
１３ シャンクロッド
１４ 削孔ビット
１５ ロッド
１６ ピストン
１７ ロータ
１８ シリンダ
Ｇ 地山
Ｍ 計測装置
Ｐ 加速度計（第１の受振器）
Ｐ１、Ｐ２ 加速度計（受振器）
Ｒ 速度計（第２の受振器）、基準速度計（基準受振器）
Ｓ 切羽
10
P1, P2 accelerometer (receiver)
R Speedometer (second geophone), Reference speedometer (reference geophone)
S face
Claims (5)
前記掘削機に設置された第１の受振器による各打撃で分割された前記弾性波の複数の振動データ、およびトンネル切羽に設置された複数の第２の受振器による各打撃ごとに分割された前記弾性波の複数の振動データを削孔長ごとに取得し、
前記第２の受振器で取得された複数の前記振動データに対して、当該振動データの相互相関関数が最大となる前記振動データを基準にして初動時刻をシフトさせた後にスタッキング処理を行い、スタッキングされた振動データを削孔長ごとに求め、
前記第１の受振器で取得された前記振動データに対して、前記第２の受振器の前記振動データに対応してシフトしてスタッキング処理を行い、スタッキングされた振動データを削孔長ごとに求め、
予め求められた前記掘削機が発振する発振時間と前記第１の受振器におけるスタッキングされた前記振動データの初動時刻を用いて、前記削孔ビットが地山を打撃した時刻でのそれぞれの前記第２の受振器におけるスタッキングされた前記振動データから初動時刻を読み取って走時曲線を求めてトモグラフィ解析を行い弾性波速度分布を求める、
ことを特徴とするトンネル切羽前方探査方法。 The rock surface in front of the tunnel face is impacted with a drilling bit provided at the tip of the rod of the drilling machine a plurality of times at predetermined time intervals at predetermined drilling lengths to generate elastic waves.
A plurality of vibration data of the elastic wave divided by each impact by the first geophone installed in the excavator, and a plurality of vibration data divided by each of the second geophones installed in the tunnel face Acquiring a plurality of vibration data of the elastic wave for each drilling length;
Stacking processing is performed on a plurality of the vibration data acquired by the second geophone after shifting the initial movement time with reference to the vibration data that maximizes the cross correlation function of the vibration data, and stacking Calculated vibration data for each drilling length,
The vibration data acquired by the first geophone is shifted in accordance with the vibration data of the second geophone to perform stacking processing, and the stacked vibration data is generated for each drilling length. Ask for
Each of the first at the time when the drilling bit strikes the ground using the oscillation time at which the drilling machine oscillates and the initial movement time of the stacked vibration data in the first geophone, which are determined in advance. The initial motion time is read from the stacked vibration data in the two geophones, the travel time curve is determined, the tomographic analysis is performed, and the elastic wave velocity distribution is determined.
A tunnel face forward search method characterized by
前記基準受振器以外の前記第２の受振器で取得された前記振動データに対しては、前記基準受振器の前記振動データに対応してシフトしてスタッキング処理を行い、スタッキングされた振動データを求める、
ことを特徴とする請求項１記載のトンネル切羽前方探査方法。 Any one of the second geophones is used as a reference geophone, and with respect to the vibration data acquired by the reference geophone, the initial movement time is shifted based on the vibration data for which the cross correlation function is maximum The stacking process is performed later to find the stacked vibration data,
The vibration data acquired by the second geophone other than the reference geophone is shifted corresponding to the vibration data of the reference geophone to perform stacking processing, and the stacked vibration data is Ask,
The tunnel face front exploration method according to claim 1 characterized by things.
ことを特徴とする請求項２記載のトンネル切羽前方探査方法。 The reference geophone is the second geophone installed at a position closest to the excavator.
The tunnel face front exploration method according to claim 2 characterized by things.
ことを特徴とする請求項１〜３の何れか一項に記載のトンネル切羽前方探査方法。 The stacking process excludes the vibration data in which the cross correlation function falls below a predetermined value.
The tunnel face front exploration method according to any one of claims 1 to 3, characterized in that:
ことを特徴とする請求項１〜４の何れか一項に記載のトンネル切羽前方探査方法。 When obtaining the travel time curve, the initial motion time of the vibration data is automatically read using a local stationary AR (autoregressive) model.
The tunnel face front exploration method according to any one of claims 1 to 4, characterized in that:
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017246137 | 2017-12-22 | ||
| JP2017246137 | 2017-12-22 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2019113514A true JP2019113514A (en) | 2019-07-11 |
| JP6945917B2 JP6945917B2 (en) | 2021-10-06 |
Family
ID=67221493
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018009056A Active JP6945917B2 (en) | 2017-12-22 | 2018-01-23 | Tunnel face forward exploration method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6945917B2 (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000170478A (en) * | 1998-12-07 | 2000-06-20 | Taisei Corp | Method for probing front of tunnel facing |
| JP2001099945A (en) * | 1999-09-30 | 2001-04-13 | Kajima Corp | Survey instrument for surveying front side of working face, and recording medium |
| JP2004286527A (en) * | 2003-03-20 | 2004-10-14 | National Institute Of Advanced Industrial & Technology | Initial motion detection method of multicomponent ae waveform |
| JP2010038790A (en) * | 2008-08-06 | 2010-02-18 | Kajima Corp | Elastic wave probe system |
| JP2014106128A (en) * | 2012-11-28 | 2014-06-09 | Hazama Ando Corp | Method of measuring natural ground elastic wave velocity |
| JP2014535044A (en) * | 2011-09-28 | 2014-12-25 | サウジ アラビアン オイル カンパニー | Reservoir characteristics prediction using least square support vector machine |
| JP2016017900A (en) * | 2014-07-10 | 2016-02-01 | 大成建設株式会社 | Method of exploring ahead of tunnel face |
| US20170199290A1 (en) * | 2014-07-21 | 2017-07-13 | Cgg Services Sas | Systems and methods for attenuating noise using interferometric estimation |
| JP2017156106A (en) * | 2016-02-29 | 2017-09-07 | 株式会社奥村組 | Tunnel face front survey method |
-
2018
- 2018-01-23 JP JP2018009056A patent/JP6945917B2/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000170478A (en) * | 1998-12-07 | 2000-06-20 | Taisei Corp | Method for probing front of tunnel facing |
| JP2001099945A (en) * | 1999-09-30 | 2001-04-13 | Kajima Corp | Survey instrument for surveying front side of working face, and recording medium |
| JP2004286527A (en) * | 2003-03-20 | 2004-10-14 | National Institute Of Advanced Industrial & Technology | Initial motion detection method of multicomponent ae waveform |
| JP2010038790A (en) * | 2008-08-06 | 2010-02-18 | Kajima Corp | Elastic wave probe system |
| JP2014535044A (en) * | 2011-09-28 | 2014-12-25 | サウジ アラビアン オイル カンパニー | Reservoir characteristics prediction using least square support vector machine |
| JP2014106128A (en) * | 2012-11-28 | 2014-06-09 | Hazama Ando Corp | Method of measuring natural ground elastic wave velocity |
| JP2016017900A (en) * | 2014-07-10 | 2016-02-01 | 大成建設株式会社 | Method of exploring ahead of tunnel face |
| US20170199290A1 (en) * | 2014-07-21 | 2017-07-13 | Cgg Services Sas | Systems and methods for attenuating noise using interferometric estimation |
| JP2017156106A (en) * | 2016-02-29 | 2017-09-07 | 株式会社奥村組 | Tunnel face front survey method |
Also Published As
| Publication number | Publication date |
|---|---|
| JP6945917B2 (en) | 2021-10-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101547508B1 (en) | Apparatus for predicting front geological features and the method thereof | |
| JP6131027B2 (en) | Measurement method of natural ground elastic wave velocity | |
| JP5839271B2 (en) | Tunnel face forward exploration method | |
| JP3380795B2 (en) | Rock bed exploration method | |
| JP6522918B2 (en) | Elastic wave velocity measurement method | |
| JP4157635B2 (en) | Tunnel face forward exploration method | |
| US7178626B2 (en) | Method of seismic evaluation of subterranean strata | |
| CN115390129A (en) | In-situ acoustic penetration device with built-in longitudinal and transverse wave transmitting and receiving transducers | |
| JP2017049198A (en) | Working face front survey device and working face front survey method | |
| JP2710740B2 (en) | Front face exploration method | |
| JP6602675B2 (en) | Geological exploration method in front of the face | |
| JP5997521B2 (en) | Face investigation method using shield machine | |
| JP4071988B2 (en) | Ground survey method using S-wave amplitude associated with impact penetration | |
| JP5258734B2 (en) | Tunnel front face exploration method and exploration system | |
| JP5940303B2 (en) | Tunnel face forward exploration method | |
| JPH09242459A (en) | Ground survey device and method | |
| JP4260329B2 (en) | Geological exploration method in front of tunnel face | |
| JPH08226975A (en) | Method for surveying geology in front of face of tunnel | |
| JP6945917B2 (en) | Tunnel face forward exploration method | |
| JP3632023B2 (en) | Ground survey method based on drilling data of rock drill | |
| JP2017141598A (en) | Teacher data making method and ground evaluation method | |
| JP5186538B2 (en) | Natural mountain exploration method | |
| JP6653072B2 (en) | Method for evaluating the number of impacts of hydraulic hammer, method for exploring front ground using the same, and system for exploring front ground | |
| Tsukamoto et al. | Predicting seismic velocity distribution in the ground ahead of tunnel face using drilling vibration of hydraulic rock drill | |
| JP2850572B2 (en) | Geological survey equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200624 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20210526 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210604 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210719 |
|
| 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: 20210913 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210913 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6945917 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |