(57)【特許請求の範囲】
【請求項１】
将来の発破予定領域に設けられる複数の起爆孔での段発発破に伴って生じる、前記発破予定領域の周辺の所定地点における段発振動波を予測する振動予測方法であって、
少なくとも１つの発破実施領域において少なくとも１つの単発発破を行い、前記発破実施領域の周辺の少なくとも１つの測定地点における基本単発振動波を測定する測定工程と、
前記測定工程で測定した前記基本単発振動波の大きさと、前記起爆孔のそれぞれの装薬量および前記起爆孔から前記所定地点までのそれぞれの距離とを用いて、前記起爆孔での単発発破に伴って生じる前記所定地点における複数の予測単発振動波の大きさを算出し、算出した複数の大きさと、前記測定工程で測定した前記基本単発振動波の波形とに基づいて、複数の前記予測単発振動波を作成する作成工程と、
前記作成工程で作成した複数の前記予測単発振動波を、設定された起爆順序および設定された起爆秒時差に従って重ね合わせることで、前記段発振動波の大きさを算出する算出工程と、
を含み、
前記測定工程では、前記発破実施領域、前記単発発破の回数、および前記測定地点のうち少なくとも１つが複数であることにより複数の前記基本単発振動波を測定し、
前記作成工程では、前記測定工程で測定した複数の前記基本単発振動波の波形の全部または一部と、前記起爆孔から前記所定地点までのそれぞれの距離とを用いて、前記複数の予測単発振動波の波形を更に算出し、算出した複数の大きさと複数の波形とに基づいて、複数の前記予測単発振動波を作成することを特徴とする振動予測方法。
【請求項２】
前記測定工程では、前記発破実施領域、前記単発発破の回数、および前記測定地点のうち少なくとも１つが複数であることにより複数の前記基本単発振動波を測定し、
前記作成工程では、複数の前記基本単発振動波の大きさを与えるような、前記発破実施領域に設けられた起爆孔の装薬量と当該起爆孔から前記測定地点までの距離との近似関数の係数を求め、求めた係数を用いて、複数の前記予測単発振動波の大きさを算出することを特徴とする請求項１に記載の振動予測方法。
【請求項３】
将来の発破予定領域に設けられる複数の起爆孔での段発発破に伴って生じる、前記発破予定領域の周辺の所定地点における段発振動波を予測する振動予測方法であって、
少なくとも１つの発破実施領域において少なくとも１つの単発発破を行い、前記発破実施領域の周辺の少なくとも１つの測定地点における基本単発振動波を複数測定する測定工程と、
前記測定工程で測定した複数の前記基本単発振動波の波形の全部または一部と、前記起爆孔から前記所定地点までのそれぞれの距離とを用いて、前記起爆孔での単発発破に伴って生じる前記所定地点における複数の予測単発振動波の波形を算出し、算出した複数の波形と、前記測定工程で測定した前記基本単発振動波の大きさとに基づいて、複数の前記予測単発振動波を作成する作成工程と、
前記作成工程で作成した複数の前記予測単発振動波を、設定された起爆順序および設定された起爆秒時差に従って重ね合わせることで、前記段発振動波の大きさを算出する算出工程と、
を含み、
前記作成工程では、前記測定工程で測定した複数の前記基本単発振動波の大きさと、前記起爆孔のそれぞれの装薬量および前記起爆孔から前記所定地点までのそれぞれの距離とを用いて、前記起爆孔での単発発破に伴って生じる前記所定地点における複数の予測単発振動波の大きさを更に算出し、算出した複数の波形と複数の大きさとに基づいて、複数の前記予測単発振動波を作成することを特徴とする振動予測方法。
【請求項４】
前記作成工程では、前記測定工程で測定した複数の前記基本単発振動波の波形のうち、前記発破実施領域に設けられた起爆孔から前記測定地点までの距離が前記所定地点までの距離にもっとも近い２つの前記基本単発振動波の波形を用い、内挿または外挿により、複数の前記予測単発振動波の波形を算出することを特徴とする請求項３に記載の振動予測方法。
【請求項５】
前記作成工程では、前記測定工程で測定した複数の前記基本単発振動波の波形における、前記発破実施領域に設けられた起爆孔から前記測定地点までの距離と、前記所定地点までの距離との差に基づく前記基本単発振動波の波形の加重平均を求めることにより、複数の前記予測単発振動波の波形を算出することを特徴とする請求項３に記載の振動予測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、発破に伴って生じる振動を予測する振動予測方法に関する。
【背景技術】
【０００２】
トンネル掘削工事では、掘削対象の岩盤（切羽）が極めて堅硬な場合、発破工法が採用される。発破工法では、切羽に複数の装薬孔を穿ち、装薬孔に爆薬を装填し、同時または順次に爆薬を爆発させる。この爆発により、対象物を破砕する。近隣に住宅や構造物がある場合、発破工法によって生じる振動の対策が重要になる。
【０００３】
下記の特許文献１，２に記載されるように、発破に伴って生じる振動を予測する技術が知られている。特許文献１に記載される工法では、試験発破による地盤振動の実測結果に基づいて、本発破としての複数回の発破を行った場合の地盤振動を予測している。本発破において生じる地盤振動の予測では、まず各発破により生じる個々の地盤振動を予測し、各発破の時間差を考慮して、これらの地盤振動を合成する。試験発破の位置は、本発破の発破孔近傍に設けられ、試験発破における装薬量および装薬深度は、本発破で予定される複数回の発破のうちのいずれか１つに合わせられる。個々の地盤振動の予測は、試験発破による地盤振動の実測結果と、装薬量と地盤振動との関係（たとえば比例関係）とに基づいて行われる。また、応用例として、複数回の試験発破に基づく振動予測を行い、この予測と同じ条件で実際に試験発破を行って地盤振動を実測し、これらの結果から補正係数を求める方法が記載されている。
【０００４】
特許文献２に記載される工法では、切羽の発破対象面設けられた複数の装薬孔の側部にスロットを形成し、スロットにもっとも近い装薬孔とそのスロットとの間隔と、装薬孔間の間隔とをいずれも同じ一定の距離とし、試験的な起爆を実施して振動波を記録している。この振動波データに基づいて、遅延秒時間隔を設定した上で、シミュレーション用振動波データを生成する。そして、これらのデータを重ね合わせて、シミュレーション合成データを生成する。遅延秒時間隔を変えて複数のシミュレーション合成データを生成し、最大振幅がもっとも小さい合成データに対応する遅延秒時間隔を得て、この遅延秒時間隔を最適な起爆秒時間隔として採用する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平１１−１８１７５３号公報
【特許文献２】特開２００１−２１３００号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上記した従来の工法では、試験発破を行って振動を測定し、この測定結果に基づいて将来の発破における振動を予測している。しかしながら、将来の発破の振動を予測する際に、装薬量と振動との関係を用いているか、または試験発破と同じ振動をシミュレーションに用いているに過ぎないため、予測対象となり得る地点は、限られた地点となってしまう。
【０００７】
本発明は、将来の発破に伴って生じる、所定地点における振動を予測することができる振動予測方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明の振動予測方法は、将来の発破予定領域に設けられる複数の起爆孔での段発発破に伴って生じる、発破予定領域の周辺の所定地点における段発振動波を予測する振動予測方法であって、少なくとも１つの発破実施領域において少なくとも１つの単発発破を行い、発破実施領域の周辺の少なくとも１つの測定地点における基本単発振動波を測定する測定工程と、測定工程で測定した基本単発振動波の大きさと、起爆孔のそれぞれの装薬量および起爆孔から所定地点までのそれぞれの距離とを用いて、起爆孔での単発発破に伴って生じる所定地点における複数の予測単発振動波の大きさを算出し、算出した複数の大きさと、測定工程で測定した基本単発振動波の波形とに基づいて、複数の予測単発振動波を作成する作成工程と、作成工程で作成した複数の予測単発振動波を、設定された起爆順序および設定された起爆秒時差に従って重ね合わせることで、段発振動波の大きさを算出する算出工程と、を含み、測定工程では、発破実施領域、単発発破の回数、および測定地点のうち少なくとも１つが複数であることにより複数の基本単発振動波を測定し、作成工程では、測定工程で測定した複数の基本単発振動波の波形の全部または一部と、起爆孔から所定地点までのそれぞれの距離とを用いて、複数の予測単発振動波の波形を更に算出し、算出した複数の大きさと複数の波形とに基づいて、複数の予測単発振動波を作成することを特徴とする。
【０００９】
この振動予測方法では、発破実施領域において単発発破が行われ、測定地点において、基本単発振動波が測定される。次に、測定した基本単発振動波の大きさ（たとえば、振幅等）に基づいて、発破予定領域での単発発破を行った場合の所定地点における複数の予測単発振動波の大きさ（たとえば、振幅等）が算出される。この予測単発振動波の算出には、発破予定領域に設けられる起爆孔のそれぞれの装薬量と、起爆孔から所定地点までの距離とが用いられる。ある地層における発破を考えた場合、発破によって生じ得る振動波の大きさは、装薬量および距離と相関を持っている。この相関を用いることで、複数の予測単発振動波の大きさを精度良く算出することができる。さらに、予測単発振動波の大きさと、基本単発振動波の波形とに基づいて、複数の予測単発振動波が作成される。このように、基本単発振動波に関する２つの側面である、大きさと波形とを用いることで、発破予定領域で発破が行われた場合の予測単発振動波を精度良く推定できる。これらの予測単発振動波が起爆順序および起爆秒時差に従って重ね合わされ、段発振動波の大きさが算出されるので、将来の発破予定領域での段発発破に伴って生じる、所定地点における振動を予測することができる。起爆孔から所定地点までのそれぞれの距離が、複数の予測単発振動波の大きさおよび波形の算出に加味されることで、複数の予測単発振動波をより精度良く算出することができる。
【００１１】
測定工程では、発破実施領域、単発発破の回数、および測定地点のうち少なくとも１つが複数であることにより複数の基本単発振動波を測定し、作成工程では、複数の基本単発振動波の大きさを与えるような、発破実施領域に設けられた起爆孔の装薬量と当該起爆孔から測定地点までの距離との近似関数の係数を求め、求めた係数を用いて、複数の予測単発振動波の大きさを算出してもよい。複数の基本単発振動波に基づいて、装薬量および距離と基本単発振動波の大きさとの相関を示す近似関数の係数が求められるため、複数の予測単発振動波の大きさをより精度良く算出することができる。
【００１２】
本発明の振動予測方法は、将来の発破予定領域に設けられる複数の起爆孔での段発発破に伴って生じる、発破予定領域の周辺の所定地点における段発振動波を予測する振動予測方法であって、少なくとも１つの発破実施領域において少なくとも１つの単発発破を行い、発破実施領域の周辺の少なくとも１つの測定地点における基本単発振動波を複数測定する測定工程と、測定工程で測定した複数の基本単発振動波の波形の全部または一部と、起爆孔から所定地点までのそれぞれの距離とを用いて、起爆孔での単発発破に伴って生じる所定地点における複数の予測単発振動波の波形を算出し、算出した複数の波形と、測定工程で測定した基本単発振動波の大きさとに基づいて、複数の予測単発振動波を作成する作成工程と、作成工程で作成した複数の予測単発振動波を、設定された起爆順序および設定された起爆秒時差に従って重ね合わせることで、段発振動波の大きさを算出する算出工程と、を含み、作成工程では、測定工程で測定した複数の基本単発振動波の大きさと、起爆孔のそれぞれの装薬量および起爆孔から所定地点までのそれぞれの距離とを用いて、起爆孔での単発発破に伴って生じる所定地点における複数の予測単発振動波の大きさを更に算出し、算出した複数の波形と複数の大きさとに基づいて、複数の予測単発振動波を作成することを特徴とする。
【００１３】
この振動予測方法では、発破実施領域において単発発破が行われ、測定地点において、複数の基本単発振動波が測定される。次に、測定した複数の基本単発振動波の波形に基づいて、発破予定領域での単発発破を行った場合の所定地点における複数の予測単発振動波の波形が算出される。この予測単発振動波の算出には、複数の基本単発振動波の波形の全部または一部と、発破予定領域に設けられる起爆孔から所定地点までのそれぞれの距離とが用いられる。ある地層における発破を考えた場合、発破によって生じ得る振動波の波形は、起爆孔からの距離と相関を持っている。この相関を用いることで、複数の予測単発振動波の波形を精度良く算出することができる。さらに、予測単発振動波の波形と、基本単発振動波の大きさとに基づいて、複数の予測単発振動波が作成される。このように、基本単発振動波に関する２つの側面である、大きさと波形とを用いることで、発破予定領域で発破が行われた場合の予測単発振動波を精度良く推定できる。これらの予測単発振動波が起爆順序および起爆秒時差に従って重ね合わされ、段発振動波の大きさが算出されるので、将来の発破予定領域での段発発破に伴って生じる、所定地点における振動を予測することができる。起爆孔から所定地点までのそれぞれの距離が、複数の予測単発振動波の波形および大きさの算出に加味されることで、複数の予測単発振動波をより精度良く算出することができる。
【００１５】
作成工程では、測定工程で測定した複数の基本単発振動波の波形のうち、発破実施領域に設けられた起爆孔から測定地点までの距離が所定地点までの距離にもっとも近い２つの基本単発振動波の波形を用い、内挿または外挿により、複数の予測単発振動波の波形を算出してもよい。発破予定領域の起爆孔から所定地点までの距離にもっとも近い２つの基本単発振動波の波形が用いられることで、複数の予測単発振動波の波形をより精度良く算出することができる。
【００１６】
作成工程では、測定工程で測定した複数の基本単発振動波の波形における、発破実施領域に設けられた起爆孔から測定地点までの距離と、所定地点までの距離との差に基づく複数の基本単発振動波の波形の加重平均を求めることにより、複数の予測単発振動波の波形を算出してもよい。起爆孔からの距離の差に基づいて複数の基本単発振動波の波形の加重平均が求められることで、複数の予測単発振動波の波形をより精度良く算出することができる。
【発明の効果】
【００１８】
本発明によれば、将来の発破に伴って生じる、所定地点における振動を予測することができる。
【図面の簡単な説明】
【００１９】
【図１】本発明の一実施形態に係る振動予測方法が使用されるトンネル掘削現場を示す模式図である。
【図２】図１のトンネル掘削現場における発破パターンを示す模式図である。
【図３】単発発破と段発発破を組み合わせた発破により得られる振動波の一例である。
【図４】振動予測方法の手順を示すフローチャートである。
【図５】図５の手順の一部の詳細を示すフローチャートである。
【図６】１つの振動波の大きさ及び波形を説明するための図である。
【図７】（ａ）は基本単発振動波の波形とその内挿によって算出される予測単発振動波の波形とを示す図、（ｂ）は複数の基本単発振動波の波形とその加重平均によって算出される予測単発振動波の波形とを示す図である。
【図８】複数の基本単発振動波の近似曲線を示す図である。
【図９】予測単発振動波の間隔が起爆秒時差にまで縮められて、重ね合わされる手順を説明するための図である。
【図１０】（ａ）は、複数の起爆孔から測定地点への振動の伝搬経路を模式的に示す図、（ｂ）は（ａ）の各振動の測定地点への到達時間を模式的に示す図である。
【図１１】２つの所定地点において、段発振動派の振動速度の予測と実測を対比した図である。
【図１２】２つの所定地点において、段発振動派の周波数特性の予測と実測を対比した図である。
【発明を実施するための形態】
【００２０】
以下、本発明の実施形態について、図面を参照しながら説明する。なお、図面の説明において同一要素には同一符号を付し、重複する説明は省略する。
【００２１】
図１に示されるように、本実施形態の振動予測方法は、たとえば、堅硬な岩盤からなる地層１における発破工法に適用される。この発破工法は、たとえば地層１におけるトンネル掘削に適用される。トンネル掘削は、所定方向（図１において矢印で示す方向）に進行する。本実施形態の振動予測方法は、トンネル掘削を行うにあたり、将来の発破予定領域である切羽Ａａでの段発発破に伴って生じる、切羽Ａａの周辺の所定地点ＰＺにおける段発振動波（たとえば振幅や振動速度等）を予測するための方法である。振動予測の対象となる所定地点は特に限定されず、切羽Ａａの周辺の任意の地点とすることができる。所定地点は、たとえば、切羽Ａａの上方の地表部Ｇ上において住宅またはその他の構造物が存在する地点である。図１に示される測定地点Ｐ１〜Ｐ３（詳しくは後述）を振動予測の対象である所定地点とすることもできる。
【００２２】
本実施形態の振動予測方法では、トンネルＴＮＬの掘削において実施する発破に伴い実際に生じる複数の振動波（振動波データ）を用いて、トンネルＴＮＬの掘削の将来の段発発破に伴って生じる振動を予測する。「振動波」は、振幅や振動速度等によって表し得る大きさの情報と、波形の情報とを含む概念である。
【００２３】
図２に示されるように、発破を行う際、切羽Ａには、複数の起爆孔１０が設けられる。複数の起爆孔１０は、切羽Ａにおいて、段発発破を行うための複数の起爆孔群１０１〜１０８を成している。図２に示される起爆孔１０の配列パターンはあくまで一例であり、個数および配列パターンは適宜変更できる（以下、起爆孔１０の個数および配列パターンを含む概念を「発破パターン」という）。複数の起爆孔１０の段数および配列は、公知の方法に基づいて決定することができる。
【００２４】
切羽Ａには、更に、単発発破を行うための１または複数の起爆孔２０が設けられる。図２に示される例では、切羽Ａの中央部に設けられた起爆孔群１０１，１０２の左右の両側に、単発発破用の起爆孔２０が１つずつ配置されている。本実施形態の振動予測方法では、実際のトンネル掘削において振動データを収集し、収集した振動データは、将来の発破に伴って生じる振動の予測に用いられる。すなわち、１または複数の単発発破用の起爆孔２０は、振動データを収集するための起爆孔である。切羽Ａの発破では、まず起爆孔２０を用いた少なくとも１回の単発発破が行われ、それに続いて、予め定められた起爆秒時差をもって、起爆孔１０（起爆孔群１０１〜１０８）を用いた段発発破が行われる（図３参照）。単発発破と段発発破との間には、たとえば１００ｍｓ以上の間隔が設けられる。
【００２５】
各起爆孔１０，２０内には、たとえば複数の爆薬（いわゆる親ダイおよび増ダイ）と、奥側の爆薬（いわゆる親ダイ）に設けられた雷管と、爆薬の手前側に設けられた込め物とが設けられている。各起爆孔１０，２０における装薬量は、同一であってもよいし、起爆孔ごとに違っていてもよい。本実施形態に適用される発破工法では、１ｍｓ単位での起爆秒時設定が可能な電子雷管が用いられる。この電子雷管は、切羽Ａでの実施工に即し、起爆秒時設定を切羽Ａにて任意に変更・設定できる機能を有している。この電子雷管では、たとえば１ｍｓ〜３０ｍｓの範囲において、起爆秒時差を１ｍｓ単位で設定可能である。電子雷管を採用することにより、高精度の起爆秒時差で段発発破を行うことができる。その結果として、所定地点ＰＺで生じる振動を精度良く予測し、振動を低減することができる。
【００２６】
図１および図４を参照して、本実施形態の振動予測方法の手順について説明する。図１に示されるように、トンネルＴＮＬの掘削では、発破実施領域である切羽Ａ１，Ａ２，Ａ３…において順次段発発破が行われる。段発発破が行われ、トンネル掘削が進行する中で、単発発破による振動波を収集する。振動波は、切羽Ａ１，Ａ２，Ａ３…の周辺の少なくとも１つの測定地点において測定される。たとえば、図１に示される例では、測定地点Ｐ１〜Ｐ３の３つの測定地点が設定されている。複数の測定地点が設けられる場合、図１に示されるようにトンネル掘削の進行方向に離間する複数の地点を設定してもよいし、当該進行方向に交差する方向に離間する複数の地点を設定してもよい。測定地点は、任意に設定することができる。測定地点Ｐ１〜Ｐ３のそれぞれには振動計を設置しておき、各振動計により、切羽Ａ１，Ａ２，Ａ３…の発破に応じて発生する振動速度（変位速度）を測定する。なお、振動計によって、振幅を測定してもよいし、振動加速度（変位加速度）を測定してもよい。
【００２７】
図４に示されるように、まず、単発発破の実施によって、任意の地点における振動波を測定する（ステップＳ１、測定工程）。図２に示される発破パターンを採用する場合、単発発破が２回行われた後、段発発破が行われる。単発発破と段発発破の間隔、すなわち２回目の単発発破と最初の起爆孔群１０１の発破との間隔は、一定以上の時間とされる。これにより、単発振動波と段発振動波とは、互いに分離した状態で取得される（図３参照）。なお、単発振動波は、単発発破に伴って測定される振動波であり、段発振動波は、段発発破に伴って測定される振動波である。図２に示されるように、最初の２回の単発発破の振動波がそれぞれ独立して得られており、各振動波の振動速度および波形が得られている。振動速度は、ゼロ点を基準とした絶対値の最大値として求められる（一例として、図３中の０．３２８ｋｉｎｅと０．１４３ｋｉｎｅ）。
【００２８】
起爆孔２０は、切羽Ａに対して１つであってもよいし、３つ以上であってもよい。切羽Ａ１，Ａ２，Ａ３…と切羽Ａが変わるごとに、起爆孔２０の位置および個数を変更してもよい。切羽Ａの中で起爆孔２０が設けられる位置は、特に限定されない。１つの測定地点Ｐ１を基準に考えた場合、同一の切羽Ａ内において、各起爆孔２０から測定地点Ｐ１までの距離は異なっている。また、Ａ１，Ａ２，Ａ３…と掘削が進行することによっても、各起爆孔２０から測定地点Ｐ１までの距離は異なる。さらに、複数の測定地点Ｐ１〜Ｐ３が設けられる場合、起爆孔２０での１回の単発発破につき、測定地点Ｐ１〜Ｐ３の箇所数に応じた複数の単発振動波が得られる。複数の単発振動波は、起爆孔２０から測定地点Ｐ１〜Ｐ３までのそれぞれの距離を振動が伝搬した結果、発生した振動波であると言える。
【００２９】
このように、ステップＳ１では、少なくとも１つの切羽Ａにおいて少なくとも１つの単発発破を行い、少なくとも１つの測定地点Ｐ１〜Ｐ３における単発振動波を測定する。ステップＳ１で測定される単発振動波を基本単発振動波と呼ぶ。この基本単発振動波の収集が行われる切羽Ａ、単発発破の回数、および測定地点のいずれもが１つである場合は、収集される基本単発振動波は１つである。一方、切羽Ａ、単発発破の回数、および測定地点のうち少なくとも１つが複数である場合、収集される基本単発振動波は複数である。基本単発振動波は１つであってもよいが、複数である方が好ましい。
【００３０】
続いて、将来の単発発破による所定地点ＰＺにおける複数の振動波を作成する（ステップＳ２、作成工程）。ステップＳ２で作成される単発振動波を予測単発振動波と呼ぶ。予測単発振動波は、基本単発振動波の大きさ及び波形と、基本単発振動波に対応する起爆孔の装薬量と、起爆孔から測定地点までの距離とに基づいて作成される。続いて、トンネルＴＮＬ掘削の将来の段発発破による、所定地点ＰＺにおける振動波を予測および評価する（ステップＳ３）。ステップＳ３で予測・評価される振動波を段発振動波と呼ぶ。以下、ステップＳ２およびステップＳ３の各手順について詳述する。
【００３１】
ステップＳ２では、図５に示されるステップＳ２１〜Ｓ２６が行われる。まず、ステップＳ１で測定した基本単発振動波の中から、２つ以上の基本単発振動波を選定する（ステップＳ２１）。ここで、図６に示されるように、振動波は、時間に対する大きさの変化として測定されるため、時間の関数で表すことができる。１つの振動波は、振幅Ｘと波形Ｙを有している。本実施形態の振動予測方法では、振動波のこれら２つの側面に着目して、予測単発振動波を算出している。なお、振幅Ｘは、振動速度や振動加速度等、他の「大きさ」の指標に置き換えて以下の手順を行うこともできる。
【００３２】
具体的には、ステップＳ２１では、ステップＳ１で測定した複数の基本単発振動波のうち、単発発破が実施された切羽Ａの起爆孔２０から測定地点Ｐ１〜Ｐ３までの距離が、切羽Ａａの各起爆孔１０から所定地点ＰＺまでの距離にもっとも近い２つの基本単発振動波を選定する。たとえば、図７（ａ）に示されるように、所定地点Ｃの距離ｄｃにもっとも近い距離を有するのが測定地点Ａ（距離ｄａを有する）と、測定地点Ｂ（距離ｄｂを有する）である場合、これらの２地点Ａ，Ｃで測定された基本単発振動波ａ（ｔ），ｂ（ｔ）が選定される。切羽Ａａに設けられる予定のすべての起爆孔１０について、基本単発振動波を選定する。
【００３３】
次に、ステップＳ２１で選定した基本単発振動波に基づいて、所定地点ＰＺにおける波形を算出する（ステップＳ２２）。ここでは、たとえば、上記のように選定した２つの基本単発振動波ａ（ｔ），ｂ（ｔ）の波形を用い、内挿または外挿により、予測単発振動波の波形を算出する。図７（ａ）に示される例では、距離ｄａと距離ｄｂの間に距離ｄｃが存在する。この場合、下記式（１）に基づいて、内挿によって予測単発振動波の波形ｃ（ｔ）が算出される。
【数１】
0006484089.tif 000002
【００３４】
なお、ステップＳ１で周波数スペクトルａ（ω），ｂ（ω）が得られる場合、これらの周波数スペクトルを用いて、下記式（２）に基づいて、内挿によって予測単発振動波の周波数スペクトルｃ（ω）を算出してもよい。
【数２】
0006484089.tif 000003
【００３５】
また、距離ｄａと距離ｄｂの間に距離ｄｃが存在する場合は、所定地点Ｃの距離ｄｃよりも大きい距離ｄｂを有する測定地点Ｂに関して、距離の差Ｌ２は負の値となる。距離ｄａと距離ｄｂの間に距離ｄｃが存在しない場合、たとえば測定地点Ｂの距離ｄｂが所定地点Ｃの距離ｄｃよりも小さい場合は、距離の差Ｌ２は正の値となり、上記式（１）は外挿を表す式となる。周波数スペクトルを用いる場合も同様である。
【００３６】
以上のステップＳ２１およびＳ２２により、２以上の基本単発振動波の波形に基づいて、予測単発振動波の波形が算出される。切羽Ａａに設けられる予定のすべての起爆孔１０について、予測単発振動波の波形が算出される。
【００３７】
なお、上記のステップＳ２１，Ｓ２２において、単発発破が実施された切羽Ａの起爆孔２０から測定地点Ｐ１〜Ｐ３までの距離と、切羽Ａａの各起爆孔１０から所定地点ＰＺまでの距離との差に基づく複数の基本単発振動波の波形の加重平均を求めることにより、複数の予測単発振動波の波形を算出してもよい。この場合、図７（ｂ）に示されるように、たとえば所定地点Ｃの距離ｄｃに距離が近い地点Ａ，Ｃ，Ｄで測定された基本単発振動波ａ（ｔ），ｂ（ｔ），ｄ（ｔ）が選定される。そして、下記式（３）に基づいて、距離の差に基づく加重平均によって、予測単発振動波の波形ｃ（ｔ）が算出される。
【数３】
0006484089.tif 000004
【００３８】
なお、ステップＳ１で周波数スペクトルａ（ω），ｂ（ω），ｄ（ω）が得られる場合、これらの周波数スペクトルを用いて、下記式（４）に基づいて、荷重平均によって予測単発振動波の周波数スペクトルｃ（ω）を算出してもよい。
【数４】
0006484089.tif 000005
【００３９】
続いて、ステップＳ１で測定した複数の基本単発振動波から、距離と振動速度の近似曲線を算出し、平均振動速度と最大振動速度に対応する係数を算出する（ステップＳ２３）。このステップＳ２３では、複数の基本単発振動波の大きさを与えるような、単発発破が実施された切羽Ａの起爆孔２０の装薬量と、当該起爆孔２０から測定地点Ｐ１〜Ｐ３までの距離との近似関数の係数を求める。近似曲線は、図８に示される近似式で近似され得る。ここで、Ｖは振動速度であり、Ｒは距離であり、Ｗは起爆孔２０の装薬量である。αおよびβは係数である。横軸ＳＤはＳｃａｌｅ  Ｄｉｓｔａｎｃｅを意味し、装薬量の０．５乗で正規化された距離である。
【００４０】
ここで、振動速度は、ある正規化距離ＳＤに対して正規分布を取ると仮定する。平均振動速度に対応する近似曲線は、正規分布における平均の近似曲線である（５０％ライン。図８において実線Ｒで示されるライン）。そこで、この平均の近似曲線の係数αおよびβを算出する。また、最大振動速度に対応する近似曲線は、正規分布におけるたとえば９５％の近似曲線である（９５％上限ライン）。そこで、この確率９５％の近似曲線の係数αおよびβを算出する。なお、図８に示される近似式に限られず、Ｖ＝Ｋ・Ｗｍ／Ｒｎなる近似式を利用してもよい（Ｋ，ｍ，ｎは係数）。このように、複数のシードウェーブ（Ｍｕｌｔｉ  Ｓｅｅｄ  Ｗａｖｅ）を用いて、振動波の大きさの分布を出す。
【００４１】
次に、算出した係数を用いて、所定地点ＰＺにおける平均振動速度と最大振動速度を算出する（ステップＳ２４）。このステップＳ２４では、ステップＳ２３で算出した平均の近似曲線の係数αおよびβと、切羽Ａａの各起爆孔１０の装薬量と、起爆孔１０から所定地点ＰＺまでの距離とを用いて、図８に示される近似式により、起爆孔１０での単発発破に伴って生じる所定地点ＰＺにおける予測単発振動波の大きさ（たとえば振動速度）を算出する。なお、切羽Ａの起爆孔１０では、実際には段発発破が行われるが、ここでの計算上、単発発破が行われるものと仮定する。また、ステップＳ２３で算出した確率９５％の近似曲線の係数αおよびβと、切羽Ａａの各起爆孔１０の装薬量と、起爆孔１０から所定地点ＰＺまでの距離とを用いて、図８に示される近似式により、起爆孔１０での単発発破に伴って生じる所定地点ＰＺにおける予測単発振動波の大きさ（たとえば振動速度）を算出する。切羽Ａａに設けられる予定のすべての起爆孔１０について、予測単発振動波の平均振動速度と最大振動速度とを算出する。
【００４２】
次に、各起爆孔１０の装薬量を適宜設定して、ステップＳ２４で算出した平均振動速度および最大振動速度となるような分布で、複数の予測単発振動波の振動速度を算出する（ステップＳ２５）。これにより、図８に示される基本単発振動波と同じ分布（正規分布）を有する複数の予測単発振動波の振動速度が算出される。以上のステップＳ２３〜Ｓ２５により、複数の基本単発振動波の大きさに基づいて、複数の予測単発振動波の大きさが算出される。切羽Ａａに設けられる予定のすべての起爆孔１０について、予測単発振動波の大きさが算出される。
【００４３】
次に、ステップＳ２２で算出した予測単発振動波の波形と、ステップＳ２５で算出した予測単発振動波の振動速度とに基づいて、予測単発振動波を作成する（ステップＳ２６、作成工程）。このステップＳ２６では、ステップＳ２４で算出した最大振動速度となるように、予測単発振動波の波形を拡幅または縮小する。
【００４４】
所定地点ＰＺにおける振動波を予測および評価するステップＳ３では、図５に示されるステップＳ３１〜Ｓ３４が行われる。まず、切羽Ａａにおける起爆孔１０の起爆順序を適宜設定し、上記のステップＳ２（ステップＳ２１〜Ｓ２６）で作成した複数の予測単発振動波を、設定した起爆順序に並べる（ステップＳ３１）。たとえば、図９の上段に示されるように、単発振動波Ｗ１，Ｗ２，Ｗ３，Ｗ４，Ｗ５が時系列で並べられる。
【００４５】
次に、起爆秒時差を適宜設定し、設定した起爆秒時差に従って、複数の単発振動波を重ね合わせる（ステップＳ３２）。ここでは、公知の手法（たとえばモンテカルロ法等）で振動波を重ね合わせることができる。たとえば、図９の中段、下段に示されるように、単発振動波Ｗ１，Ｗ２，Ｗ３，Ｗ４，Ｗ５が重ね合わされて、段発振動波Ｗｔが得られる。
【００４６】
そして、ステップＳ３２で算出された段発振動波Ｗｔの最大振動速度を算出する（ステップＳ３３、算出工程）。ここでは、図６に示されるように、ゼロ点を基準とした絶対値の最大値が、最大振動速度として算出される。算出された最大振動速度は、予測最大振動速度である。
【００４７】
さらに、ステップＳ３３で算出した予測最大振動速度は許容値以下であるか否かを判断する（ステップＳ３４）。予測最大振動速度は許容値以下であると判断される場合（ステップＳ３４；Ｙｅｓ）、振動予測は終了する。予測最大振動速度は許容値を超えていると判断される場合（ステップＳ３４；Ｎｏ）、そのような段発発破を切羽Ａａで行うことはできない。そこでステップＳ２４に戻り、切羽Ａａの発破パターン、各起爆孔１０の装薬量、起爆孔１０の発破順序、起爆秒時差を再設定することにより、予測単発振動波の予測・評価を繰り返す。
【００４８】
なお、切羽Ａａの発破パターン、各起爆孔１０の装薬量、起爆孔１０の発破順序、起爆秒時差を変更して、ステップＳ２１〜ステップＳ３３の予測単発振動派の算出手順を繰り返し、複数回シミュレーションしてもよい。本実施形態の振動予測方法では、複数回のシミュレーションにより、将来の発破で生じ得る、もっとも確からしいシードウェーブを作り出すことが望ましい。
【００４９】
また、図１０（ａ）に示されるように、切羽Ａと測定地点との間に２種類の地質Ｉ，ＩＩがある場合は、各起爆孔１０の位置と測定地点とを結ぶ線（伝搬経路）がいずれの地質を主に伝搬するかを考慮し、地質によって異なる伝搬速度を考慮してもよい。たとえば、Ａの起爆孔１０は、主として地質Ｉ内を伝搬する。Ｂの起爆孔１０は、地質Ｉと地質ＩＩをほぼ均等な長さだけ伝搬する。Ｃの起爆孔１０は、主として地質ＩＩ内を伝搬する。このように、複数の予測単発振動波を重ね合わせる際、地質の違いによる伝搬速度の違いを考慮してもよい。図１０（ｂ）に示されるように、振動が測定地点に到達する時間は、伝搬経路と地質との関係で異なっている。このように、地質が異なる場合は伝搬速度も異なると考えられるため、上記の方法によれば、複数の予測単発振動波をより精度良く算出することができる。なお、伝搬速度の推定には、たとえばトンネルトモグラフィを利用することができる。
【００５０】
以上説明した本実施形態の振動予測方法によれば、図１に示されるように、発破実施領域である切羽Ａ１，Ａ２，Ａ３…において単発発破が行われ、測定地点Ｐ１〜Ｐ３において、複数の基本単発振動波が測定される。次に、測定した複数の基本単発振動波の波形および大きさに基づいて、発破予定領域である切羽Ａａでの単発発破を行った場合の所定地点ＰＺにおける複数の予測単発振動波の波形および大きさが算出される。この予測単発振動波の算出には、複数の基本単発振動波の波形の全部または一部と、切羽Ａａに設けられる起爆孔１０から所定地点ＰＺまでのそれぞれの距離とが用いられる。ある地層１における発破を考えた場合、発破によって生じ得る振動波の波形および大きさは、起爆孔１０からの距離と相関を持っている。この相関を用いることで、複数の予測単発振動波の波形および大きさを精度良く算出することができる。さらに、予測単発振動波の波形および大きさに基づいて、複数の予測単発振動波が作成される。このように、基本単発振動波に関する２つの側面である、波形と大きさとを用いることで、切羽Ａａで発破が行われた場合の予測単発振動波を精度良く推定できる。これらの予測単発振動波が起爆順序および起爆秒時差に従って重ね合わされ、段発振動波の大きさが算出されるので、切羽Ａａ（将来の発破予定領域）での段発発破に伴って生じる、所定地点ＰＺにおける振動を予測することができる。
【００５１】
また、複数の基本単発振動波に基づいて、装薬量および距離と基本単発振動波の大きさとの相関を示す近似関数の係数α，βが求められるため、複数の予測単発振動波の大きさをより精度良く算出することができる。
【００５２】
もっとも近い２つの基本単発振動波の波形の内挿または外挿により、複数の予測単発振動波の波形を算出するので、切羽Ａａの起爆孔１０から所定地点ＰＺまでの距離にもっとも近い２つの基本単発振動波の波形が用いられる。よって、複数の予測単発振動波の波形をより精度良く算出することができる。
【００５３】
なお、距離の差に基づく複数の基本単発振動波の波形の加重平均を求める場合も、複数の予測単発振動波の波形をより精度良く算出することができる。
【００５４】
以上、本発明の一実施形態について説明したが、本発明は上記実施形態に限られない。たとえば、上記実施形態では、複数の基本単発振動波の大きさを用いて予測単発振動波の大きさを推定したが、予測単発振動波の大きさを推定する際に、１つの基本単発振動波の大きさを用いてもよい。予測単発振動波の大きさまたは波形のいずれかを算出せずに、基本単発振動波の大きさまたは波形のいずれかをそのまま用いてもよい。
【００５５】
本発明は、トンネル掘削に限られず、他の発破工事にも適用可能である。たとえば、本発明を明かり工事に適用してもよい。発破実施領域と発破予定領域とは、同じ地層に設けられていることが好ましい。
【００５６】
次に、図１１および図１２を参照して、本発明の実施例による振動予測の試験結果について説明する。図１１は、２つの所定地点において、段発振動派の振動速度の予測と実測を対比した図である。図１２は、２つの所定地点において、段発振動派の周波数特性の予測と実測を対比した図である。この試験では、実際のトンネル掘削における波形予測の運用を考慮し，１〜５回までの発破で収録した基本単発振動波を基に、６回目の段発発破（段発振動波）の推定を行った．図１１および図１２に示されるように、波形、周波数特性ともに、実測と予測は概ね良く一致していることが確認された。
【符号の説明】
【００５７】
１…地層、１０…起爆孔、２０…起爆孔、Ａ、Ａ１，Ａ２，Ａ３…切羽（発破実施領域）、Ａａ…切羽（発破予定領域）、Ｐ１〜Ｐ３…測定地点、ＰＺ…所定地点、ＴＮＬ…トンネル、Ｗ１〜Ｗ５…単発振動波、Ｗｔ…段発振動波、Ｘ…振幅（大きさ）、Ｙ…波形。
【図１】
0006484089.tif 000006
【図２】
0006484089.tif 000007
【図３】
0006484089.tif 000008
【図４】
0006484089.tif 000009
【図５】
0006484089.tif 000010
【図６】
0006484089.tif 000011
【図７】
0006484089.tif 000012
【図８】
0006484089.tif 000013
【図９】
0006484089.tif 000014
【図１０】
0006484089.tif 000015
【図１１】
0006484089.tif 000016
【図１２】
0006484089.tif 000017