JPH0814480B2 - 発破方法 - Google Patents
発破方法Info
- Publication number
- JPH0814480B2 JPH0814480B2 JP1121420A JP12142089A JPH0814480B2 JP H0814480 B2 JPH0814480 B2 JP H0814480B2 JP 1121420 A JP1121420 A JP 1121420A JP 12142089 A JP12142089 A JP 12142089A JP H0814480 B2 JPH0814480 B2 JP H0814480B2
- Authority
- JP
- Japan
- Prior art keywords
- blasting
- blast
- time series
- vibration
- sequence
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
本発明は、多数の爆薬を起爆させるに当たり、周辺へ
の駆音及び振動等の弊害少なく起爆させることができる
発破方法に関する。
の駆音及び振動等の弊害少なく起爆させることができる
発破方法に関する。
従来のいわゆる多段式発破は、多数の爆薬を一定の時
間間隔または全くランダムな時間間隔で起爆させてい
た。第２図はトンネル空間をつくる場合の従来の等間隔
発破パターンを示す。この図で丸印は爆薬のセット位置
を示し、それを接続する線は同時に起爆されるグループ
を示す。従来では、、……の時系列で全て等時
間間隔に起爆させていた。 いま第４図に示すようにＮ段の発破の発破時系列を等
間隔τとした場合、音源（音圧）もしくは振動源（圧
力）の時間波形Ｐ（ｔ）は次式で表される。 Ｐ（ｔ）＝as（ｔ）＋（ｔ−τ）＋as（ｔ−２τ）＋…
…＋as｛（ｔ−（Ｎ−１）τ｝ ……（１） ここでａは発破の規模に比例した定数である。 （１）式をフーリエ変換すると、 ここでωは角周波数、δｎは発破時系列、 である。 すなわち、Ｋ（ω）は発破時系列δｎによって決まる
スペクトルである。例えばωが零もしくは２π／τとな
れば、 Ｋ（ω）＝Ｎ ……（４） となり、一発の破砕のＮ倍の圧力となる。特にω＝２π
／τのときのスペクトルピークは甚大となる。 第６図は等時間間隔で起爆させた場合のモデル図で、
同図（ａ）に示す音源（音圧）もしくは振動源（圧力）
の時間波形Ｐ（ｔ）のように、例えば８個の爆薬を一定
の時間間隔τで起爆させたとき、各爆破によって伝搬す
る音圧もしくは振動をそれぞれ（ｂ），（ｃ）……
（ｉ）とすると、観測される時間波形Pm（ｔ）は（ｊ）
のようになる。
間間隔または全くランダムな時間間隔で起爆させてい
た。第２図はトンネル空間をつくる場合の従来の等間隔
発破パターンを示す。この図で丸印は爆薬のセット位置
を示し、それを接続する線は同時に起爆されるグループ
を示す。従来では、、……の時系列で全て等時
間間隔に起爆させていた。 いま第４図に示すようにＮ段の発破の発破時系列を等
間隔τとした場合、音源（音圧）もしくは振動源（圧
力）の時間波形Ｐ（ｔ）は次式で表される。 Ｐ（ｔ）＝as（ｔ）＋（ｔ−τ）＋as（ｔ−２τ）＋…
…＋as｛（ｔ−（Ｎ−１）τ｝ ……（１） ここでａは発破の規模に比例した定数である。 （１）式をフーリエ変換すると、 ここでωは角周波数、δｎは発破時系列、 である。 すなわち、Ｋ（ω）は発破時系列δｎによって決まる
スペクトルである。例えばωが零もしくは２π／τとな
れば、 Ｋ（ω）＝Ｎ ……（４） となり、一発の破砕のＮ倍の圧力となる。特にω＝２π
／τのときのスペクトルピークは甚大となる。 第６図は等時間間隔で起爆させた場合のモデル図で、
同図（ａ）に示す音源（音圧）もしくは振動源（圧力）
の時間波形Ｐ（ｔ）のように、例えば８個の爆薬を一定
の時間間隔τで起爆させたとき、各爆破によって伝搬す
る音圧もしくは振動をそれぞれ（ｂ），（ｃ）……
（ｉ）とすると、観測される時間波形Pm（ｔ）は（ｊ）
のようになる。
図６（ｋ）は（ｊ）をフーリエ変換して周波数特性を
みたものであるが、多段発破の時間間隔τに相当する周
波数1/τ（Hz）のピークが大きい周波数特性となってお
り、このような特性を持つ音波が周辺民家に伝搬するこ
とで建具をがたつかせるなどの被害をもたらしていた。 本発明は、多段発破においてもその基本時間間隔τに
相当する周波数1/τ（Hz）のピークを極力小さくできる
斬新な発破方法を提供することを目的とする。
みたものであるが、多段発破の時間間隔τに相当する周
波数1/τ（Hz）のピークが大きい周波数特性となってお
り、このような特性を持つ音波が周辺民家に伝搬するこ
とで建具をがたつかせるなどの被害をもたらしていた。 本発明は、多段発破においてもその基本時間間隔τに
相当する周波数1/τ（Hz）のピークを極力小さくできる
斬新な発破方法を提供することを目的とする。
本発明の発破方法は、基本周期をτとするＭ系列（Ma
ximum unll period sequenceもしくはMaximum−length
linear shift register sequenceの略称）に従って多数
の爆薬を起爆することを特徴とする。 Ｍ系列は情報処理分野等で利用されている人為的な時
系列の一つであり、ここでいう時系列とは、確立的な現
像の時間的変化を連続的または一定間隔をおいて不連続
に観測して得た系列をいう。Ｍ系列は一般にはPn−１の
周期をもつ系列が考えられるが、コンピュータ利用の利
便性から2n−１の周期をもつ２値（０と１）を利用す
る。 一般に周期Ｎ＝2n−１のＭ系列は、ｎ次の既約多項式
でしかも素なものを使って、次のようにすればよい。疎
な既約多項式をｆ（ｘ）として、 の関係を、和を“２を法とする和”に置き換え、冪を添
字とみなし、 とする。これをXkについて解いた式 Xk＝C1Xk-1＋C2Kk-2＋……＋CnXk-n ……（７） で与えられ、初期値として、Xk（Ｋ＝0,……ｎ−１）に
はすべてが０以外なら何を与えてもよい。表１は本発明
で用いるＭ系列を発生するための既約多項式の一例であ
る。 例えば、ｎ＝５の場合、Ｍ系列を（８）式を用いて０
または１の値をとるXkの値を次々に求める。ここで＋の
信号は２を法とする和を表すものとする（１を“真”、
０を“偽”に対応させればこれは排他的論理和となり計
算により簡単に求めることができる）。 Xk＝Xk-2＋Xk-5 ……（８） 初期値として、X0＝０、X1＝０、X2＝０、X3＝０、X4
＝１とおいてみると、 なる系列が得られ、これは周期25−１＝31をもってい
る。並んだ５つの組合わせは00000以外のすべてが一度
ずつ表れている。この系列を装置で発生させるには、５
個の遅延要素（shift register）と排他的論理和の回路
で第12図のように構成すればよい。Ｍ系列は、このよう
にｎ個の記憶要素と線形な演算回路とで構成できる種々
の系列の中で最大の周期をもっており、最大周期列とを
呼ばれる。
ximum unll period sequenceもしくはMaximum−length
linear shift register sequenceの略称）に従って多数
の爆薬を起爆することを特徴とする。 Ｍ系列は情報処理分野等で利用されている人為的な時
系列の一つであり、ここでいう時系列とは、確立的な現
像の時間的変化を連続的または一定間隔をおいて不連続
に観測して得た系列をいう。Ｍ系列は一般にはPn−１の
周期をもつ系列が考えられるが、コンピュータ利用の利
便性から2n−１の周期をもつ２値（０と１）を利用す
る。 一般に周期Ｎ＝2n−１のＭ系列は、ｎ次の既約多項式
でしかも素なものを使って、次のようにすればよい。疎
な既約多項式をｆ（ｘ）として、 の関係を、和を“２を法とする和”に置き換え、冪を添
字とみなし、 とする。これをXkについて解いた式 Xk＝C1Xk-1＋C2Kk-2＋……＋CnXk-n ……（７） で与えられ、初期値として、Xk（Ｋ＝0,……ｎ−１）に
はすべてが０以外なら何を与えてもよい。表１は本発明
で用いるＭ系列を発生するための既約多項式の一例であ
る。 例えば、ｎ＝５の場合、Ｍ系列を（８）式を用いて０
または１の値をとるXkの値を次々に求める。ここで＋の
信号は２を法とする和を表すものとする（１を“真”、
０を“偽”に対応させればこれは排他的論理和となり計
算により簡単に求めることができる）。 Xk＝Xk-2＋Xk-5 ……（８） 初期値として、X0＝０、X1＝０、X2＝０、X3＝０、X4
＝１とおいてみると、 なる系列が得られ、これは周期25−１＝31をもってい
る。並んだ５つの組合わせは00000以外のすべてが一度
ずつ表れている。この系列を装置で発生させるには、５
個の遅延要素（shift register）と排他的論理和の回路
で第12図のように構成すればよい。Ｍ系列は、このよう
にｎ個の記憶要素と線形な演算回路とで構成できる種々
の系列の中で最大の周期をもっており、最大周期列とを
呼ばれる。
すなわち、本発明は上記（２）式におけるδｎとして
Ｍ系列またはその一部を使用する。この発破方法の原理
について説明すると、例えば起爆段階が８段の爆破で
は、Ｍ系列として、N/2＋１＝（2n−１）/2＋１≧８と
なるｎ次の既約多項式を用いる。例えば、ｎ＝４とし
て、表１よりＭ系列を発生する既約多項式として（９）
式を用いると、 ｆ（ｘ）＝X4＋X3＋１ ……（９） これから、 Xk＝Xk-3＋Xk-4 ……（10） となり、初期値として、X0＝１、X1＝１、X2＝１、X3＝
１とすると、 となる。 ここで１のときのみ起爆するものとすると第５図のよ
うな発破時系列となる。 基本周期τ＝100msとすれば、爆破時系列は、 0,100,200,300,700,1000,1100,1300……（ms） となる。 第３図はこの方法で第２図に示したときと同様のトン
ネルを造る例の発破パターンを示し、，，，，
，，，，，，の時系列で起爆する。 第７図（ａ）はこの場合のモデル図で、各段の爆発に
よって伝搬する音圧もしくは振動をそれぞれ（ｂ）、
（ｃ）……（ｉ）とすると、観測される時間波形Pm
（ｔ）は（ｊ）のようになる。第７図（ｋ）は同（ｊ）
をフーリエ変換して周波数特性をみたものであるが、第
６図に示す通常の発破を比べて、多段発破の基本周期τ
に相当する周波数1/τ（Hz）のピークが小さくなる。 なお、０と１の組合わせで求めたＭ系列では、１のと
きばかりでなく０のときに起爆させてもよい。例えば前
述の例では、 であるから、基本周期τ＝100msとすると、 0,100,200,400,500,800,1000,1500……（ms） となる。 このように多段の発破の爆破時系列をＭ系列とするこ
とで、爆破の基本周期τに依存するスペクトルの卓越を
軽減し、あらゆる角種は数ωにおいて一様なスペクトル
にすることが可能である。すなわち、 となり、上記（４）式とのレベル差をとれば、その効果
Ｅは、 を得る。例えば、発破の段数がＮ＝８の場合、Ｅ＝−1
7.6dbとなる。Ｎが大きいと、 となる。
Ｍ系列またはその一部を使用する。この発破方法の原理
について説明すると、例えば起爆段階が８段の爆破で
は、Ｍ系列として、N/2＋１＝（2n−１）/2＋１≧８と
なるｎ次の既約多項式を用いる。例えば、ｎ＝４とし
て、表１よりＭ系列を発生する既約多項式として（９）
式を用いると、 ｆ（ｘ）＝X4＋X3＋１ ……（９） これから、 Xk＝Xk-3＋Xk-4 ……（10） となり、初期値として、X0＝１、X1＝１、X2＝１、X3＝
１とすると、 となる。 ここで１のときのみ起爆するものとすると第５図のよ
うな発破時系列となる。 基本周期τ＝100msとすれば、爆破時系列は、 0,100,200,300,700,1000,1100,1300……（ms） となる。 第３図はこの方法で第２図に示したときと同様のトン
ネルを造る例の発破パターンを示し、，，，，
，，，，，，の時系列で起爆する。 第７図（ａ）はこの場合のモデル図で、各段の爆発に
よって伝搬する音圧もしくは振動をそれぞれ（ｂ）、
（ｃ）……（ｉ）とすると、観測される時間波形Pm
（ｔ）は（ｊ）のようになる。第７図（ｋ）は同（ｊ）
をフーリエ変換して周波数特性をみたものであるが、第
６図に示す通常の発破を比べて、多段発破の基本周期τ
に相当する周波数1/τ（Hz）のピークが小さくなる。 なお、０と１の組合わせで求めたＭ系列では、１のと
きばかりでなく０のときに起爆させてもよい。例えば前
述の例では、 であるから、基本周期τ＝100msとすると、 0,100,200,400,500,800,1000,1500……（ms） となる。 このように多段の発破の爆破時系列をＭ系列とするこ
とで、爆破の基本周期τに依存するスペクトルの卓越を
軽減し、あらゆる角種は数ωにおいて一様なスペクトル
にすることが可能である。すなわち、 となり、上記（４）式とのレベル差をとれば、その効果
Ｅは、 を得る。例えば、発破の段数がＮ＝８の場合、Ｅ＝−1
7.6dbとなる。Ｎが大きいと、 となる。
次に、本発明の実施例について図面を参照に詳述す
る。 上述した通り、発破時系列δｎによって生ずるスペク
トルＫ（ω）は上記（３）式で表される。本発明は、基
本的にはδｎとしてＭ系列信号またはその一部分を採用
することによって、多段発破の重畳スペクトルピークを
低減しようとするものであるが、あらゆる周波数におい
て衝撃波を最小限に留めるためには、次式を満足するよ
うなＫ（ω）を見出す必要がある。 Ｓ（ω）Ｋ（ω）Ｈ（ω）Ｑ（ω）＝const（min）……
（14） すなわち、 となる。ここで、 Ｓ（ω）：予め観測して求める単一発破の音・振動スペ
クトル； Ｈ（ω）：発破位置から観測点までの音・振動の伝達関
数（予め観測できる）； Ｑ（ω）：人または建物等の音・振動感度を表す伝達関
数（既存のデータから推定する）； である。この場合、問題となる周波数範囲に限って|T
（ω）｜≠０なる条件を付加しておく。 従って、あらゆる周波数において衝撃波を最小限に留
めるためには、上記（15）のＫ（ω）が近似的に満たさ
れるような時系列δｎを選ぶことに帰着する。 第１図はこのような思想に基づき、コンピュータで使
用して発破時系列を制御する本発明の一例のシステム構
成図である。先ず、Ｋ（ω）を除く上記ファクタＳ
（ω）,H（ω）,Q（ω）を予め決定するため、現場で単
一の発破をかけ、任意の観測点に設置した音・振動検出
部（マイクロホンまたは振動ピックアップ手段）１でそ
の伝播してきた音・振動を電気信号として検出し、これ
をA/Dコンバータ２によってデジタル信号に変換して取
り込む。そのデジタル量から、音・振動スペクトルＳ
（ω）測定部３で単一発破の音・振動スペクトルＳ
（ω）を測定するとともに、伝達関数Ｈ（ω）測定部４
で発破位置から観測点までの伝達関数Ｈ（ω）を測定す
る。なお、音・振動検出部１の伝達関数は便宜上理想的
で一定と仮定する。 測定したＳ（ω）とＨ（ω）に基づいて、人または建
物等の音・振動感度を表す伝達関数Ｑ（ω）を推定す
る。すなわち、既存のデータから数段階の伝達関数Q
1（ω）,Q2（ω）……QX（ω）を想定してこれをメモリ
５に予め記憶しておき、そのうちからＳ（ω）及びＨ
（ω）の測定値に適合する伝達関数Ｑ（ω）を伝達関数
Ｑ（ω）選択部６において選択する。 一方、キーボード７から入力した発破段数Ｎから、前
述のようにＭ系列発生部８においてＭ系列信号を発生す
る。そして、基本周期τを例えば段階的に代えたτ1,τ
２……τｘごとに、Ｍ系列による発破時系列データ（時
間間隔はそれぞれの基本周期の整数倍）を発破時系列生
成部９で生成するとともに、その各発破時系列によって
生ずるスペクトルデータK1（ω）,K2（ω）……K
x（ω）をスペクトル生成部10において上記（３）式に
基づいて生成し、これらのデータをメモリ５に予め記憶
しておく。 上記のように測定したＳ（ω）とＨ（ω）及び選択し
たＱ（ω）から単一発破の総合スペクトルＴ（ω） Ｔ（ω）＝Ｓ（ω）Ｈ（ω）Ｑ（ω） 及びその逆関数T-1（ω）をＴ（ω）及びT-1（ω）演算
部11において演算する。この後、上記（８）式を満足す
るＫ（ω）をＫ（ω）演算部12において上記（15）式に
従って求め、その求めたＫ（ω）と、メモリ５に記憶さ
れているスペクトルデータK1（ω）,K2（ω）……K
x（ω）とを照合し、これらスペクトルデータのうちか
らＫ（ω）に最も近似するスペクトルデータを探す。そ
して、そのスペクトルデータに対応する発破時系列デー
タを発破時系列選択部13によってメモリ５から選択す
る。最後に、その選択した発破時系列データをキーボー
ド７からの指示に従い出力し、D/Aコンバータ14によっ
てアナログ電気信号に変換して起爆装置15へ入力し、選
択した発破時系列δｎに従って起爆させる。 第９図及び第11図は、上記第３図に示した発破パター
ンにより基本周期を250msとしてＭ系列に従って起爆
（各発破点の火薬量0.8kg）させ、異なる観測点Ａ（ト
ンネル抗口外約5km）と観測点Ｂ（同約150km）で音圧を
観測し、その波形をフーリエ変換したパワースペクトル
の波形図である。なお、この場合はＳ（ω）,H（ω）,Q
（ω）については何等考慮していない。 第８図及び第10図は上記第２図の発破パターンにより
等時間間隔（250ms）で従来通りに起爆させ、同様に観
測したパワースペクトルの波形図である。 これらの図を比較すれば分かるように、いずれの場合
もＡ点、Ｂ点ともに約1.5Hzを基本周期として3.0,4.5,
6.0,7.5及び9.0（Hz）でピークが現れている。最大値は
従来法では3.0Hz前後、本発明の場合は4.5Hz前後にみら
れる。パワースペクトルの形は本発明の場合が鋭いが、
その数が従来法よりも増えている。また、一般に10Hz以
下の低周波領域で本発明の場合が従来法よりもレベルが
低く、Ａ点では3.0Hz付近で本発明の方がレベルが平均
約6dB低い。また、Ｂ点においては、本発明の方が1.5Hz
付近で5dB程度、3.0Hz付近で10dB程度低くなっている。 この実測データは上述した理論と一致するものであ
り、従って発破時間間隔をＭ系列に従って基本周期の整
数倍とし、さらに第１図に示した実施例のように上記Ｓ
（ω）,H（ω）,Q（ω）のファクタを加味した発破時系
列制御を行うことによって、特に周辺への被害が大きい
低周波領域で衝撃波を一層効果的に低減できるものであ
る。
る。 上述した通り、発破時系列δｎによって生ずるスペク
トルＫ（ω）は上記（３）式で表される。本発明は、基
本的にはδｎとしてＭ系列信号またはその一部分を採用
することによって、多段発破の重畳スペクトルピークを
低減しようとするものであるが、あらゆる周波数におい
て衝撃波を最小限に留めるためには、次式を満足するよ
うなＫ（ω）を見出す必要がある。 Ｓ（ω）Ｋ（ω）Ｈ（ω）Ｑ（ω）＝const（min）……
（14） すなわち、 となる。ここで、 Ｓ（ω）：予め観測して求める単一発破の音・振動スペ
クトル； Ｈ（ω）：発破位置から観測点までの音・振動の伝達関
数（予め観測できる）； Ｑ（ω）：人または建物等の音・振動感度を表す伝達関
数（既存のデータから推定する）； である。この場合、問題となる周波数範囲に限って|T
（ω）｜≠０なる条件を付加しておく。 従って、あらゆる周波数において衝撃波を最小限に留
めるためには、上記（15）のＫ（ω）が近似的に満たさ
れるような時系列δｎを選ぶことに帰着する。 第１図はこのような思想に基づき、コンピュータで使
用して発破時系列を制御する本発明の一例のシステム構
成図である。先ず、Ｋ（ω）を除く上記ファクタＳ
（ω）,H（ω）,Q（ω）を予め決定するため、現場で単
一の発破をかけ、任意の観測点に設置した音・振動検出
部（マイクロホンまたは振動ピックアップ手段）１でそ
の伝播してきた音・振動を電気信号として検出し、これ
をA/Dコンバータ２によってデジタル信号に変換して取
り込む。そのデジタル量から、音・振動スペクトルＳ
（ω）測定部３で単一発破の音・振動スペクトルＳ
（ω）を測定するとともに、伝達関数Ｈ（ω）測定部４
で発破位置から観測点までの伝達関数Ｈ（ω）を測定す
る。なお、音・振動検出部１の伝達関数は便宜上理想的
で一定と仮定する。 測定したＳ（ω）とＨ（ω）に基づいて、人または建
物等の音・振動感度を表す伝達関数Ｑ（ω）を推定す
る。すなわち、既存のデータから数段階の伝達関数Q
1（ω）,Q2（ω）……QX（ω）を想定してこれをメモリ
５に予め記憶しておき、そのうちからＳ（ω）及びＨ
（ω）の測定値に適合する伝達関数Ｑ（ω）を伝達関数
Ｑ（ω）選択部６において選択する。 一方、キーボード７から入力した発破段数Ｎから、前
述のようにＭ系列発生部８においてＭ系列信号を発生す
る。そして、基本周期τを例えば段階的に代えたτ1,τ
２……τｘごとに、Ｍ系列による発破時系列データ（時
間間隔はそれぞれの基本周期の整数倍）を発破時系列生
成部９で生成するとともに、その各発破時系列によって
生ずるスペクトルデータK1（ω）,K2（ω）……K
x（ω）をスペクトル生成部10において上記（３）式に
基づいて生成し、これらのデータをメモリ５に予め記憶
しておく。 上記のように測定したＳ（ω）とＨ（ω）及び選択し
たＱ（ω）から単一発破の総合スペクトルＴ（ω） Ｔ（ω）＝Ｓ（ω）Ｈ（ω）Ｑ（ω） 及びその逆関数T-1（ω）をＴ（ω）及びT-1（ω）演算
部11において演算する。この後、上記（８）式を満足す
るＫ（ω）をＫ（ω）演算部12において上記（15）式に
従って求め、その求めたＫ（ω）と、メモリ５に記憶さ
れているスペクトルデータK1（ω）,K2（ω）……K
x（ω）とを照合し、これらスペクトルデータのうちか
らＫ（ω）に最も近似するスペクトルデータを探す。そ
して、そのスペクトルデータに対応する発破時系列デー
タを発破時系列選択部13によってメモリ５から選択す
る。最後に、その選択した発破時系列データをキーボー
ド７からの指示に従い出力し、D/Aコンバータ14によっ
てアナログ電気信号に変換して起爆装置15へ入力し、選
択した発破時系列δｎに従って起爆させる。 第９図及び第11図は、上記第３図に示した発破パター
ンにより基本周期を250msとしてＭ系列に従って起爆
（各発破点の火薬量0.8kg）させ、異なる観測点Ａ（ト
ンネル抗口外約5km）と観測点Ｂ（同約150km）で音圧を
観測し、その波形をフーリエ変換したパワースペクトル
の波形図である。なお、この場合はＳ（ω）,H（ω）,Q
（ω）については何等考慮していない。 第８図及び第10図は上記第２図の発破パターンにより
等時間間隔（250ms）で従来通りに起爆させ、同様に観
測したパワースペクトルの波形図である。 これらの図を比較すれば分かるように、いずれの場合
もＡ点、Ｂ点ともに約1.5Hzを基本周期として3.0,4.5,
6.0,7.5及び9.0（Hz）でピークが現れている。最大値は
従来法では3.0Hz前後、本発明の場合は4.5Hz前後にみら
れる。パワースペクトルの形は本発明の場合が鋭いが、
その数が従来法よりも増えている。また、一般に10Hz以
下の低周波領域で本発明の場合が従来法よりもレベルが
低く、Ａ点では3.0Hz付近で本発明の方がレベルが平均
約6dB低い。また、Ｂ点においては、本発明の方が1.5Hz
付近で5dB程度、3.0Hz付近で10dB程度低くなっている。 この実測データは上述した理論と一致するものであ
り、従って発破時間間隔をＭ系列に従って基本周期の整
数倍とし、さらに第１図に示した実施例のように上記Ｓ
（ω）,H（ω）,Q（ω）のファクタを加味した発破時系
列制御を行うことによって、特に周辺への被害が大きい
低周波領域で衝撃波を一層効果的に低減できるものであ
る。
本発明によれば次のような効果が期待される。 爆破工事において特に周辺住民や建物等に対し被害
の大きかった低周波騒音及び振動を発生源で軽減でき
る。従って、爆破工事の際に周辺住民の苦情から工事の
進歩に支障をきたすなどの場合があったが、発生源で有
効な対策を講じることができ、工期の短縮が図れ、また
被害補償等の問題もなくなる。 作業員に与える騒音・振動の影響も軽減でき、その
防止のための施設費を低減できる。 起爆装置自体は既製のものをそのまま使用できる。 コンピュータ制御によって被害の極力少ない有効な
発破制御が行える。 地質の調査や探査等にも応用できる。
の大きかった低周波騒音及び振動を発生源で軽減でき
る。従って、爆破工事の際に周辺住民の苦情から工事の
進歩に支障をきたすなどの場合があったが、発生源で有
効な対策を講じることができ、工期の短縮が図れ、また
被害補償等の問題もなくなる。 作業員に与える騒音・振動の影響も軽減でき、その
防止のための施設費を低減できる。 起爆装置自体は既製のものをそのまま使用できる。 コンピュータ制御によって被害の極力少ない有効な
発破制御が行える。 地質の調査や探査等にも応用できる。
第１図は本発明の発破方法を実施するシステム構成の一
例を示すブロック図、第２図は従来の等時間間隔発破法
（以下、従来法と記す）によるトンネル工事の際の発破
パターンの説明図、第３図は本発明による場合の発破パ
ターンの説明図、第４図は従来法の発破時系列説明図、
第５図は本発明の発破時系列説明図、第６図（ａ）〜
（ｋ）は従来法の場合の振動モデルのタイムチャート、
第７図（ａ）〜（ｋ）は本発明の場合の振動モデルのタ
イムチャート、第８図及び第10図は従来法によって起爆
させて観測した音圧のパワースペクトル波形図、第９図
及び第11図は本発明により起爆させて観測した音圧のパ
ワースペクトル波形図、第12図はＭ系列のつくる方法を
示す論理回路のブロック図である。
例を示すブロック図、第２図は従来の等時間間隔発破法
（以下、従来法と記す）によるトンネル工事の際の発破
パターンの説明図、第３図は本発明による場合の発破パ
ターンの説明図、第４図は従来法の発破時系列説明図、
第５図は本発明の発破時系列説明図、第６図（ａ）〜
（ｋ）は従来法の場合の振動モデルのタイムチャート、
第７図（ａ）〜（ｋ）は本発明の場合の振動モデルのタ
イムチャート、第８図及び第10図は従来法によって起爆
させて観測した音圧のパワースペクトル波形図、第９図
及び第11図は本発明により起爆させて観測した音圧のパ
ワースペクトル波形図、第12図はＭ系列のつくる方法を
示す論理回路のブロック図である。
Claims (1)
- 【請求項１】多数の爆薬を基本周期τに対してＭ系列に
従った整数倍の時間間隔で起爆することを特徴とする発
破方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1121420A JPH0814480B2 (ja) | 1989-05-17 | 1989-05-17 | 発破方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1121420A JPH0814480B2 (ja) | 1989-05-17 | 1989-05-17 | 発破方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH02302600A JPH02302600A (ja) | 1990-12-14 |
| JPH0814480B2 true JPH0814480B2 (ja) | 1996-02-14 |
Family
ID=14810711
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1121420A Expired - Lifetime JPH0814480B2 (ja) | 1989-05-17 | 1989-05-17 | 発破方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0814480B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1998021544A1 (fr) * | 1996-11-12 | 1998-05-22 | Asahi Kasei Kogyo Kabushiki Kaisya | Procede d'excavation par explosion |
-
1989
- 1989-05-17 JP JP1121420A patent/JPH0814480B2/ja not_active Expired - Lifetime
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1998021544A1 (fr) * | 1996-11-12 | 1998-05-22 | Asahi Kasei Kogyo Kabushiki Kaisya | Procede d'excavation par explosion |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH02302600A (ja) | 1990-12-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Haskell | Total energy and energy spectral density of elastic wave radiation from propagating faults: part II. A statistical source model | |
| CN102135445B (zh) | 爆破振动预测方法 | |
| Liu et al. | A one-dimensional meshfree particle formulation for simulating shock waves | |
| de Hoop et al. | Generation of acoustic waves by an impulsive line source in a fluid/solid configuration with a plane boundary | |
| Ragam et al. | Evaluation and prediction of blast-induced peak particle velocity using artificial neural network: A case study | |
| Tang | Predictive processing of array acoustic waveform data | |
| CA1059611A (en) | Method of seismic exploration | |
| JPWO1998021544A1 (ja) | 発破工法 | |
| CN109827481B (zh) | 一种微差干扰降振爆破中最优微差时间确定方法及装置 | |
| JP6484089B2 (ja) | 振動予測方法 | |
| JP6420054B2 (ja) | 弾性波速度測定方法 | |
| EP0323687A1 (en) | Method of measuring, analysing, predicting and controlling vibrations induced by explosive blasting in earth formations | |
| JPH0814480B2 (ja) | 発破方法 | |
| US20210183403A1 (en) | Frequency extraction method using dj transform | |
| Clark et al. | Plane and spherical waves in a Voigt medium | |
| JP6379261B2 (ja) | 弾性波速度測定方法および弾性波速度測定システム | |
| JPH11181753A (ja) | 発破による地盤振動の予測方法および発破方法 | |
| CN120339536B (zh) | 一种用于一次成井拉槽爆破数据分析处理方法以及系统 | |
| US3432807A (en) | Signal interpretation | |
| CN118962781B (zh) | 矿震震源反演方法以及矿震震源反演系统 | |
| SU584426A1 (ru) | Имитатор радиосигнала | |
| JP3127088B2 (ja) | 比抵抗トモグラフィー法測定装置 | |
| Bhagwat et al. | Comparison of some blast vibration predictors for blasting in underground drifts and some observations | |
| Makarov et al. | On the possibility of determining internal wave characteristics from the ray arrival time distribution in an underwater sound channel under conditions of ray chaos | |
| Vu et al. | Developing the combined Kalman-EMD algorithm to process blasting shockwave signals propagating in a water medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090214 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100214 Year of fee payment: 14 |
|
| EXPY | Cancellation because of completion of term | ||
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100214 Year of fee payment: 14 |