JP2007232621A - Electromagnetic underground prospecting method and system for the same - Google Patents
Electromagnetic underground prospecting method and system for the same Download PDFInfo
- Publication number
- JP2007232621A JP2007232621A JP2006056226A JP2006056226A JP2007232621A JP 2007232621 A JP2007232621 A JP 2007232621A JP 2006056226 A JP2006056226 A JP 2006056226A JP 2006056226 A JP2006056226 A JP 2006056226A JP 2007232621 A JP2007232621 A JP 2007232621A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- measurement
- amplitude modulation
- information
- wave
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000005259 measurement Methods 0.000 claims abstract description 119
- 230000005540 biological transmission Effects 0.000 claims abstract description 34
- 239000013598 vector Substances 0.000 claims description 31
- 230000004044 response Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 10
- 238000001514 detection method Methods 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims 1
- 230000014509 gene expression Effects 0.000 description 25
- 238000006243 chemical reaction Methods 0.000 description 14
- 230000000875 corresponding effect Effects 0.000 description 12
- 238000004891 communication Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000009466 transformation Effects 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000003786 synthesis reaction Methods 0.000 description 3
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 2
- 230000002411 adverse Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 210000003462 vein Anatomy 0.000 description 2
- 239000004020 conductor Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000005191 phase separation Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
この発明は電磁的地下探査技術に関するもので、特に、所定の発信局から空中に放射された電波によって生じる1次磁場と当該電波に対する地下構造の応答によって生じる2次磁場との合成に相当する変動磁場を測定することにより地下構造の情報を得る技術に関する。 The present invention relates to an electromagnetic underground exploration technique, and in particular, a variation corresponding to a combination of a primary magnetic field generated by a radio wave radiated from a predetermined transmitting station into the air and a secondary magnetic field generated by a response of the underground structure to the radio wave. The present invention relates to a technique for obtaining information on an underground structure by measuring a magnetic field.
地層構造や地中の鉱脈などの物理探査法として電磁的地下探査方法が知られている。これは、鉱脈などの探査対象物の導電的性質が周囲の地層とは異なることを利用して、地表側から電磁波を送信するとともに地下構造の電磁的応答を測定し、その測定結果から探査対象物の存在場所や分布形態を推定する方法である。 Electromagnetic underground exploration methods are known as geophysical exploration methods such as the stratum structure and underground veins. This utilizes the fact that the conductive properties of exploration objects such as mineral veins are different from the surrounding geological formations, and transmits electromagnetic waves from the ground surface and measures the electromagnetic response of the underground structure, and from the measurement results, This is a method for estimating the location and distribution form of an object.
可搬の送信機と受信機との組合せを用いた電磁的地下探査方法についての先行特許文献としては、たとえば特許文献1や特許文献2がある。
For example, there are
一方、電磁波の送信機能をインフラによって代用するという発想から、海上自衛隊および米軍が対潜水艦通信用に全地球的に設置しているVLF(超長波)通信局から発信されているVLF電波の搬送波を受動的に利用する方法(VLF法)が知られている。VLF法は、電場と磁場とを測定するVLF−MT法(VLF−MT:Very Low Frequency-Magnetotelluri:地磁気地電流法)と、磁場のみを測定するVLF−EM法(Very Low Frequency-Electromagnetic:VLF−電磁法)とに分類されるが、後者は地表と非接触で測定が可能という利点がある。 On the other hand, from the idea of substituting the electromagnetic wave transmission function by infrastructure, the VLF radio wave carrier transmitted from the VLF (very long wave) communication stations installed by the JMSDF and the US Army globally for anti-submarine communication A method (VLF method) that passively uses the above is known. The VLF method includes a VLF-MT method (VLF-MT: Very Low Frequency-Magnetotelluri) that measures an electric field and a magnetic field, and a VLF-EM method (Very Low Frequency-Electromagnetic: VLF) that measures only a magnetic field. -Electromagnetic method), but the latter has the advantage that it can be measured without contact with the ground surface.
ところで従来のVLF法では、広域にわたって探査を行うためには測定エリア内に仮想的に設定した測線上を等間隔で移動しつつ測定を繰り返す必要がある。このため、この方法では測定機材とともに長距離にわたって地上を移動せねばならないだけでなく、地形的な制約があったり地上に人工障害物が存在したりすると、測定場所への移動が不可能ないしは困難になる。また移動そのものは可能であっても広域にわたって地上を移動しつつ測定を繰り返すためには多大の時間と労力を要する。 By the way, in the conventional VLF method, in order to search over a wide area, it is necessary to repeat measurement while moving at equal intervals on a survey line virtually set in the measurement area. For this reason, this method not only has to move over the ground for a long distance together with the measurement equipment, but it is impossible or difficult to move to the measurement location if there are topographical restrictions or there are artificial obstacles on the ground. become. Even if movement is possible, it takes a lot of time and labor to repeat measurements while moving on the ground over a wide area.
また、従来のLVF法では測定精度も不十分になりがちである。すなわち、上記のVLF通信局は本来の目的が軍事通信であるため、実際の通信中には搬送波が変調される。そして、これによってチューニングオフを生じることがあり、安定した測定ができないことがある。 In addition, the measurement accuracy tends to be insufficient in the conventional LVF method. That is, since the above-mentioned VLF communication station is originally intended for military communication, the carrier wave is modulated during actual communication. As a result, tuning off may occur, and stable measurement may not be possible.
さらに、広域調査の場合には測定点の間隔をあまり狭くすると多大な労力を要するため、測定点の間隔はかなり広くならざるを得ない。すると、それぞれの部分では測線密度の疎な測線データとなり、測定結果の解釈が難しくなる場合が少なくない。 Furthermore, in the case of a wide-area survey, if the interval between measurement points is too narrow, a great deal of labor is required, and thus the interval between measurement points must be considerably wide. Then, in each part, it becomes survey line data with a sparse line density, and it is often difficult to interpret the measurement result.
[発明の目的]
この発明は、従来技術における上記の問題の克服を意図しており、所定の発信局からの電波を用いて受動的に地下構造の情報を得るにあたって、安定的に高速かつ精度よく測定できる技術を提供することを目的とする。
[Object of invention]
The present invention is intended to overcome the above-mentioned problems in the prior art, and in passively obtaining information on the underground structure using radio waves from a predetermined transmission station, a technique capable of stably measuring at high speed and with high accuracy. The purpose is to provide.
上述の目的を達成するため、請求項1の発明では、所定の発信局から空中に放射された電波によって生じる1次磁場と前記電波に対する地下構造の応答によって生じる2次磁場との合成に相当する変動磁場を測定することにより地下構造の情報を得る方法であって、前記電波として、時刻報知のために所定の時間周期ごとに所定の規則で振幅変調を行いつつ発信局から常時に公衆送信される送信波を利用する。
In order to achieve the above object, the invention of
そして前記方法が、独立な3方向についての磁場測定手段を保持した飛行体を測定対象エリアの上空を移動させつつ、前記3方向についての変動磁場を前記磁場測定手段によって測定する測定工程と、前記時間周期ごとに生じる前記送信波の振幅変調期間のそれぞれにつき、前記変動磁場の測定の結果から、
1) 前記3方向の変動磁場のそれぞれの振幅の情報と、
2) 前記3方向の変動磁場の相対位相の情報と、
を求める特定工程と、前記振幅の情報と前記相対位相の情報とから、前記3方向の変動磁場を成分とする3次元磁場ベクトルが3次元磁場空間において動く扁平楕円面の特徴情報を、前記振幅変調期間のそれぞれについて求める算出工程とを備え、それぞれの振幅変調期間において前記飛行体が飛行していた各測定場所付近の地下構造についての情報を、各振幅変調期間に対応して得られた扁平楕円面の特徴情報に基づいて得ることを特徴とする。
The method includes a measuring step of measuring a magnetic field measuring in three directions by the magnetic field measuring means while moving the flying object holding the magnetic field measuring means in three independent directions over the measurement target area; For each of the amplitude modulation periods of the transmitted wave that occurs every time period, from the results of the measurement of the varying magnetic field,
1) Information on the amplitude of each of the three-direction fluctuating magnetic fields;
2) Information on the relative phase of the three-direction varying magnetic field;
The characteristic information of the flat ellipsoidal surface in which a three-dimensional magnetic field vector having the three-direction varying magnetic field as a component is obtained from the specific step of obtaining the amplitude information and the relative phase information. A calculation step for each of the modulation periods, and information about the underground structure near each measurement place where the flying object flew in each amplitude modulation period is obtained according to each amplitude modulation period. It is obtained based on the feature information of the ellipsoid.
請求項2の発明は、請求項1の方法において、前記測定工程が、前記送信波を受信して、前記送信波のそれぞれの振幅変調期間の開始タイミングを検出するタイミング検出工程と、それぞれの振幅変調期間の開始タイミングに同期して前記磁場測定手段の測定結果の取り込みを開始させる工程とを備えることと特徴とする。 According to a second aspect of the present invention, in the method of the first aspect, the measurement step receives the transmission wave and detects a start timing of each amplitude modulation period of the transmission wave, and each amplitude. And a step of starting capturing of the measurement result of the magnetic field measuring means in synchronization with the start timing of the modulation period.
請求項3の発明は、請求項1または請求項2の方法において、前記算出工程が、各測定場所での変動磁場測定結果に基づいてそれぞれ実行されることにより、各測定場所ごとに前記扁平楕円面の楕円軸に関する特徴情報が得られることを特徴とする。 According to a third aspect of the present invention, in the method of the first or second aspect, the calculation step is performed based on a measurement result of a varying magnetic field at each measurement location, whereby the flat ellipse is measured at each measurement location. Characteristic information about the elliptical axis of the surface is obtained.
請求項4の発明は、請求項1ないし請求項3のいずれかの方法において、前記磁場測定手段に相対的に固定された姿勢関係で3軸傾斜計が前記飛行体に保持されており、前記3軸傾斜計の検出結果に基づいて、前記磁場測定手段に相対的に固定された局所座標系と地上に固定された絶対座標系との姿勢関係を補正する工程を、前記算出工程に付随して実行することを特徴とする。 According to a fourth aspect of the present invention, in the method according to any one of the first to third aspects, a three-axis inclinometer is held on the flying body in a posture relationship relatively fixed to the magnetic field measuring means, A step of correcting the posture relationship between the local coordinate system fixed relative to the magnetic field measuring means and the absolute coordinate system fixed on the ground based on the detection result of the three-axis inclinometer is attached to the calculation step. It is characterized by executing.
請求項5の発明は、所定の発信局から空中に放射された電波によって生じる1次磁場と前記電波に対する地下構造の応答によって生じる2次磁場との合成に相当する変動磁場を測定することにより地下構造の情報を得る電磁的地下探査システムであって、前記電波として、時刻報知のために所定の時間周期ごとに所定の規則で振幅変調を行いつつ発信局から常時に公衆送信される送信波を利用するものである。 According to a fifth aspect of the present invention, an underground magnetic field is measured by measuring a variable magnetic field corresponding to a combination of a primary magnetic field generated by a radio wave radiated from a predetermined transmitting station into the air and a secondary magnetic field generated by a response of the underground structure to the radio wave. An electromagnetic underground exploration system for obtaining structural information, wherein the radio wave is a transmission wave that is transmitted from the transmitting station to the public at all times while performing amplitude modulation according to a predetermined rule every predetermined time period for time notification. It is what you use.
そして、前記システムが、独立な3方向についての変動磁場を空中において測定する磁場測定手段と、前記磁場測定測定手段の出力から、前記振幅変調期間における波形を抽出する抽出手段と、前記波形に基づいて
1) 前記3方向の変動磁場のそれぞれの振幅の情報と、
2) 前記3方向の変動磁場の相対位相の情報と、
を求める振幅位相特定手段と、前記振幅の情報と前記相対位相の情報とから、前記3方向の変動磁場を成分とする3次元磁場ベクトルが3次元磁場空間において動く扁平楕円面の特徴情報を、前記振幅変調期間のそれぞれについて求める算出手段とを備え、それぞれの振幅変調期間において前記飛行体が飛行していた各測定場所付近の地下構造についての情報を、各振幅変調期間に対応して得られた扁平楕円面の特徴情報に基づいて得ることを特徴とする。
Then, the system is based on the magnetic field measuring means for measuring the variable magnetic field in three independent directions in the air, the extracting means for extracting the waveform in the amplitude modulation period from the output of the magnetic field measuring and measuring means, and the waveform. The
1) Information on the amplitude of each of the three-direction fluctuating magnetic fields;
2) Information on the relative phase of the three-direction varying magnetic field;
From the amplitude phase specifying means for obtaining the characteristic information of the flat ellipsoidal surface in which the three-dimensional magnetic field vector having the three-direction variable magnetic field as a component moves from the amplitude information and the relative phase information. Calculation means for obtaining each of the amplitude modulation periods, and information on the underground structure near each measurement location where the flying object flew in each amplitude modulation period can be obtained corresponding to each amplitude modulation period. It is obtained based on the feature information of the flat ellipsoid.
請求項6の発明は、請求項5記載のシステムにおいて、前記送信波を受信して、前記時間周期ごとに生じる前記送信波の振幅変調期間を検知する検知手段をさらに備え、前記抽出手段は、前記検知手段によって検知された振幅変調期間から前記波形を抽出することを特徴とする。 The invention according to claim 6 is the system according to claim 5, further comprising detection means for receiving the transmission wave and detecting an amplitude modulation period of the transmission wave that occurs every time period, and the extraction means includes: The waveform is extracted from an amplitude modulation period detected by the detecting means.
請求項7の発明は、請求項5または請求項6のシステムにおいて、前記算出手段が、各測定場所での変動磁場測定結果に基づいて、各測定場所ごとに前記扁平楕円面の楕円軸に関する特徴情報を得ることを特徴とする。 According to a seventh aspect of the present invention, in the system according to the fifth or sixth aspect, the calculating means relates to the elliptical axis of the flat elliptical surface for each measurement location based on the measurement result of the varying magnetic field at each measurement location. It is characterized by obtaining information.
請求項8の発明は、請求項5ないし請求項7のいずれかのシステムにおいて、前記磁場測定手段に相対的に固定された姿勢関係で3軸傾斜計が前記飛行体に保持されており、前記3軸傾斜計の検出結果に基づいて、前記磁場測定手段に相対的に固定された局所座標系と地上に固定された絶対座標系との姿勢関係を補正する補正手段をさらに備えることを特徴とする。 According to an eighth aspect of the present invention, in the system according to any one of the fifth to seventh aspects, a three-axis inclinometer is held on the flying body in a posture relationship relatively fixed to the magnetic field measuring means, And a correction means for correcting the posture relationship between the local coordinate system fixed relative to the magnetic field measurement means and the absolute coordinate system fixed on the ground based on the detection result of the triaxial inclinometer. To do.
請求項1〜請求項8の発明によれば、時刻報知のために所定の時間周期ごとに所定の規則で振幅変調を行いつつ発信局から常時に公衆送信される送信波のうちで、周期的な振幅変調期間の部分を利用し、その期間の発信電波によって生じる変動磁場を3方向につき測定し、その測定結果に基づいて測定場所付近の地下構造の情報を得るため、チューニングオフを生じることもなく、安定的で精度よく測定を行うことが可能となる。 According to the first to eighth aspects of the present invention, among the transmission waves constantly transmitted from the source station to the public while performing amplitude modulation according to a predetermined rule every predetermined time period for time notification, A variable amplitude modulation period is used to measure the variable magnetic field generated by the transmitted radio waves in that period in three directions, and to obtain information on the underground structure near the measurement location based on the measurement results, tuning off may occur. Therefore, stable and accurate measurement can be performed.
また、空中を飛行しつつこのような3方向の磁場測定を行うことにより、地形上の制約や地上の人工障害物の影響を受けることなく高速に測定を進めることが可能となる。 Further, by performing such three-direction magnetic field measurement while flying in the air, it becomes possible to proceed with the measurement at high speed without being affected by topographical constraints or artificial obstacles on the ground.
特に請求項2および請求項6の発明では、時刻報知のための前記送信波を受信し、その送信波のそれぞれの振幅変調期間の開始タイミングを検出するとともに、その開始タイミングに同期して磁場測定手段の測定結果の取り込みを開始させるため、特に安定かつ高精度な測定が可能となる。 In particular, in the second and sixth aspects of the invention, the transmission wave for time notification is received, the start timing of each amplitude modulation period of the transmission wave is detected, and the magnetic field measurement is synchronized with the start timing. Since the acquisition of the measurement result of the means is started, particularly stable and highly accurate measurement is possible.
さらに特に請求項4および請求項8の発明では、磁場測定手段に相対的に固定された姿勢関係で3軸傾斜計を用いて座標系の姿勢関係を補正することにより、飛行体の姿勢の変動の影響を補償して特に高精度な測定が可能となる。 More particularly, in the inventions of claims 4 and 8, the attitude of the flying object is changed by correcting the attitude relation of the coordinate system using a triaxial inclinometer with the attitude relation fixed relative to the magnetic field measuring means. In particular, it is possible to perform highly accurate measurement by compensating for the influence of the above.
<●1.システムの概要>
図1は、この発明の実施形態である電磁的地下探査方法を実現するシステム100の全体構成を示す模式図である。このシステム100では、固定された発信局200から公衆送信される送信波SWを利用する。この送信波SWは、時刻報知のために所定の時間周期ごとに所定の規則で振幅変調されて常時に公衆送信される電波であり、この実施形態で利用する送信波SWは、我が国の電波時計に使用されているJJY標準電波である。
<● 1. System overview>
FIG. 1 is a schematic diagram showing an overall configuration of a
JJY標準電波は、独立行政法人情報通信研究機構日本標準時グループ(東京都小金井市)の原子時計と同期して大鷹鳥谷山(福島県)および羽金山(佐賀県)の両発信局から日本標準時を公衆に報知するために常時に送信される40kHzおよび60kHzの電波(長波)であり、受信範囲は日本全国をカバーしている。上記の電波名称における「JJY」は、これらの発信局のコールサインである。これと同様な標準時間の電波報知システムはいずれの国・地域においても設置可能であり、我が国以外でもそれぞれの標準電波が運用されている。 JJY standard time signals are transmitted from both Otaka Toriyayama (Fukushima Prefecture) and Haganeyama (Saga Prefecture) transmitting stations in sync with the atomic clock of the Japan Standard Time Group (Koganei City, Tokyo). 40kHz and 60kHz radio waves (long waves) transmitted at all times to notify the public, and the reception range covers the whole of Japan. “JJY” in the above radio wave name is a call sign of these transmitting stations. The same standard time radio wave notification system can be installed in any country / region, and each standard radio wave is operated outside of Japan.
なお、この明細書においては、地上の固定施設としての無線発信所を特に「発信局」と呼ぶことによって、可搬の送信機と区別する。 In this specification, a radio transmitting station as a fixed facility on the ground is particularly called a “transmitting station” to distinguish it from a portable transmitter.
一方、飛行体としての無線制御の無人ヘリコプター(以下「ラジコンへリ」)1には後述する測定部2が吊り下げられないしは搭載されており、ラジコンへリ1は、地上のラジオコントローラ1aによる遠隔操作によって、地下構造の探査対象となる測定エリアARの上空を飛行することができる。一般にラジコンヘリのコントロール電波はメガヘルツ(MHz)帯であり、JJY電波は長波帯であるため、ラジコンヘリ1のコントロール電波が地下構造の測定に悪影響を及ぼすことはない。また、この実施形態のシステム100では、ラジオコントローラ1aを操作するオペレータ1Pは必要であるが、ラジコンヘリ1そのものについての搭乗員は不要である。
On the other hand, a radio controlled unmanned helicopter (hereinafter referred to as “radio control helix”) 1 as a flying object is suspended or mounted with a
地上側にはまた、測定部2によって測定した結果を解析する解析部3が準備される。この実施形態では解析部3は可搬とされて測定エリアARの近傍に設置されているが、測定エリアARとは異なる場所に固定的に設けられた解析センターに解析部3が置かれていてもよい。また、飛行体そのものに空間的余裕がある場合は、飛行体の内部に解析部3を搭載していてもかまわない。
On the ground side, an
発信局200から空中に放射された送信波SWは、その直接波として空中に1次磁場を生じさせるほか、この送信波SWに対する地下構造の応答によって空中に2次磁場が生じる。そしてこれらの合成に相当する変動磁場を測定部2によって測定する。解析部2を用いた解析によって地下構造の情報を得ることができる。
The transmission wave SW radiated from the transmitting
測定部2は、後述するように独立した3方法についての変動磁場を空中で測定可能とされており、ここで得られた測定結果はデジタルデータの形で解析部3へ受け渡される。このようなデータの授受は、測定部2で得られた測定結果を記憶メディアに記憶させておいて、全ての測定が完了した後、あるいは測定エリアARのある程度の部分の測定が完了した都度にラジコンヘリ1を着陸させ、その記憶メディアをオフラインで解析部3に与えることによって行ってもよい。また、ラジコンヘリ1の測定部2と地上の解析部3とを無線データ通信によってオンラインで結び、それによってデータ伝送を行ってもよい。デジタルデータ伝送に使用される電波も一般に高い周波数を持つから、このようなオンライン構成としても、地下構造の測定に悪影響を及ぼすことはない。
As will be described later, the
そして、ラジコンヘリ1は、測定エリアARに想定された測線PLに沿って交互飛行しながら測定エリアARを走査し、測線PL上の多数の離散的な測定場所Pi(i=1,2,3,…)上で測定を実行する。測線PLはラジコンヘリ1が飛行すべき軌跡としておおまかに想定されたものであり、それぞれの測定場所Piの位置も、そのような飛行を行いつつ一連の測定を行った結果として特定される地点である。
Then, the radio controlled
<●2.標準電波の構成>
図2は、JJY標準電波のタイムコード図である。JJY標準電波においては、40kHzまたは60kHzの搬送波が常時発信されており、1秒(1sec)おきにこの搬送波の振幅変調が開始されて振幅変調パルスとなる。ただし、ここで言う「パルス」は振幅の変調に用いられる波形がパルスであるという意味であり、JJY標準電波の搬送波の周波数はこれらの振幅変調期間であるか否かにかかわらず、40kHzまたは60kHに固定されている。この事情が図2中の部分CAに模式的に描かれている。
<● 2. Standard radio wave configuration>
FIG. 2 is a time code diagram of the JJY standard radio wave. In the JJY standard radio wave, a carrier wave of 40 kHz or 60 kHz is constantly transmitted, and amplitude modulation of this carrier wave is started every 1 second (1 sec) to become an amplitude modulation pulse. However, “pulse” here means that the waveform used for amplitude modulation is a pulse, and the frequency of the carrier wave of the JJY standard radio wave is 40 kHz or 60 kHz regardless of whether or not these amplitude modulation periods. It is fixed to. This situation is schematically depicted in a portion CA in FIG.
1回の振幅変調が持続している時間幅(単位変調期間)は、0.2sec、0.5sec、0.8secの3種類があり、これらのうち0.5sec、0.8secの単位変調期間を1secごとに所定の規則で組み合わせて配置することにより、1日のうちのどの時刻を指示しているかがコード化されている。残りの1種類である0.2secの単位変調期間は上記2種類の単位変調期間の組合せコードを区切るマーカとして機能する。 There are three types of time width (unit modulation period) in which one amplitude modulation lasts: 0.2 sec, 0.5 sec, and 0.8 sec. Of these, 0.5 sec and 0.8 sec unit modulation periods are predetermined every 1 sec. By combining them according to the above rules, it is coded which time of the day is indicated. The remaining one type of 0.2 sec unit modulation period functions as a marker for separating the combination codes of the two types of unit modulation periods.
送信波としてのJJY標準電波においては、各単位変調期間の組合せや順序は所定の時刻表示コード化規則に従って変化するが、「1secおきに搬送波の振幅変調期間が開始される」という周期性は、どの時間帯でも正確に保たれている。したがって、1secおきの各振幅変調パルス内で、40kHzまたは60kHzの振幅変調信号を用いた短時間の電磁法を実施することが可能である。 In the JJY standard radio wave as a transmission wave, the combination and order of the unit modulation periods change according to a predetermined time display coding rule, but the periodicity that “the amplitude modulation period of the carrier wave starts every 1 second” is It is kept accurate at all times. Therefore, it is possible to carry out a short-time electromagnetic method using an amplitude modulation signal of 40 kHz or 60 kHz within each amplitude modulation pulse every 1 sec.
このようなことから、この発明の実施形態ではラジコンヘリ1で飛行を行いつつJJY標準電波の振幅変調期間の開始に同期して1secおきに3方向の磁場測定を行い、それぞれの測定場所での地下構造の比抵抗を特定するように構成される。
For this reason, in the embodiment of the present invention, the magnetic field measurement in three directions is performed every 1 second in synchronization with the start of the amplitude modulation period of the JJY standard radio wave while flying with the radio controlled
以下ではまずその原理について説明し、具体的なシステム構成および動作についてはは後述する。 In the following, the principle will be described first, and the specific system configuration and operation will be described later.
<●3.測定および解析原理>
一般に、時刻tで測定された周波数fの磁場のx,y,zの各成分Hx(t),Hy(t),Hz(t) は、式1のように書ける。ただし、x,y,zは3次元直交座標であり、φkはk成分(k=x,y,z)についての初期位相である。ここでは水平面をxy面、鉛直軸をz軸としている。
<● 3. Principle of measurement and analysis>
In general, the components H x (t), H y (t), and H z (t) of x, y, and z of the magnetic field of frequency f measured at time t can be written as
また、角周波数ω=2πfを定義し、さらに式2の瞬時位相θを定義すれば、式1は
式3のように書き換えることができる。
Further, if the angular frequency ω = 2πf is defined and the instantaneous phase θ of
これらHk(θ)を3成分とする式4の3次元磁場ベクトルH(θ)#の終点は、時刻tの進行に伴って、3次元磁場空間中の扁平楕円面Ep(図3)の周Ec上を動く。 The end point of the three-dimensional magnetic field vector H (θ) # of Equation 4 having these three components of H k (θ) is the flat elliptical surface Ep (FIG. 3) in the three-dimensional magnetic field space as time t progresses. Move on Zhou Ec.
ただし、記号「T」は行と列の転置を表す。また、この明細書において右肩に記号「#」が付いた量はベクトルを、右肩に記号「##」が付いた量はマトリクスをそれぞれ示すものとする。 However, the symbol “T” represents transposition of rows and columns. In this specification, a quantity with a symbol “#” on the right shoulder indicates a vector, and a quantity with a symbol “##” on the right shoulder indicates a matrix.
扁平楕円面Epは、3次元磁場空間中で定義される扁平楕円体を、座標原点(扁平楕円体の中心)Oを通る面で切断したときの断面に相当する。便宜上、以下では周Ecを「扁平楕円周」とも呼ぶ。 The flat ellipsoidal surface Ep corresponds to a cross section when a flat ellipsoid defined in the three-dimensional magnetic field space is cut by a plane passing through the coordinate origin (center of the flat ellipsoid) O. For convenience, the circumference Ec is also referred to as a “flat elliptical circumference” below.
3次元磁場ベクトルH(θ)#は、
(1) 発信局200から送信された送信波SW(JJY標準電波)を空中で直接に受信することによって生じる1次磁場と、
(2) 送信波SWに応答して地下構造中に生じた渦電流などによって空中に生じる2次磁場と、
の合成に相当する変動磁場である。
The three-dimensional magnetic field vector H (θ) # is
(1) a primary magnetic field generated by directly receiving a transmission wave SW (JJY standard radio wave) transmitted from the
(2) a secondary magnetic field generated in the air by an eddy current generated in the underground structure in response to the transmission wave SW;
It is a fluctuating magnetic field corresponding to the synthesis of
すなわち、測定場所(たとえば測定場所Pi)の上空において磁場を測定すると、送信波SWに関して同相成分と離相成分(90度ずれた成分)との合成磁場が測定される。そして、測定場所付近の地下構造に導電性物質が多く存在したり、構造的な特異部分が多く存在するほど離相成分が増大し、扁平楕円面Epの長軸Amaxに対する短軸Aminの相対的長さが増大するとともに、これらの軸が座標軸から傾いてくる。このため、扁平楕円面Epの長軸Amaxおよび短軸Aminの長さや傾きを知ることによって地下構造の情報を得ることができる。 That is, when the magnetic field is measured above the measurement location (for example, measurement location Pi), the combined magnetic field of the in-phase component and the phase-separated component (component shifted by 90 degrees) is measured with respect to the transmission wave SW. The more the conductive material is present in the underground structure near the measurement site or the more the structural peculiar part is present, the larger the phase separation component is, and the minor axis A min with respect to the major axis A max of the flat elliptical surface Ep is increased. As the relative length increases, these axes tilt from the coordinate axes. For this reason, information on the underground structure can be obtained by knowing the length and inclination of the major axis A max and the minor axis A min of the flat elliptical surface Ep.
そのような情報の典型的なものは地下構造の比抵抗である。扁平楕円面Epの楕円率e(長軸Amaxと短軸Aminとの長さの比)と傾き角(伏角)αとの比の値(e/α)が、その付近の地下構造の誘導定数ないしは比抵抗と相関しているという事実が電磁的地下探査の分野で公知であり、この関係から地下構造の比抵抗を知ることができる。 A typical example of such information is the resistivity of underground structures. The value (e / α) of the ratio (e / α) of the ellipticity e (the ratio of the length between the major axis A max and the minor axis A min ) and the inclination angle (the depression angle) α of the flat elliptical surface Ep is The fact that it is correlated with the induction constant or resistivity is known in the field of electromagnetic underground exploration, and the resistivity of the underground structure can be known from this relationship.
ところで、扁平楕円周Ecは磁場の3成分の初期位相φkのすべてに個別に依存するのではなく、式5a〜式5cに示すような、それらの2成分ごとの位相差φxy,φyz,φzxにのみ依存する。これは、磁場の各成分の位相差が回転磁場の軌跡を決定するという事実に対応する。ただし、式6の関係があるから、これら3つの位相差φxy,φyz,φzxのうち独立な成分は2つである。 By the way, the flat elliptical circumference Ec does not depend on all of the initial phases φ k of the three components of the magnetic field individually, but the phase differences φ xy , φ yz of these two components as shown in the equations 5a to 5c. , Φ zx only. This corresponds to the fact that the phase difference of each component of the magnetic field determines the locus of the rotating magnetic field. However, because of the relationship of Equation 6, there are two independent components of these three phase differences φ xy , φ yz , and φ zx .
したがって、扁平楕円面Epの長軸Amaxおよび短軸Aminの長さおよび方向を求めるにあたっては、独立な2つの位相差としてたとえば式7b,式7cのΔφy,Δφzを使用し、残りの1つである式7aのΔφxはゼロと置くことができる(これはx成分の初期位相を基準にして考えていることに相当する)。 Therefore, in obtaining the length and direction of the major axis A max and the minor axis A min of the flat elliptical surface Ep, for example, Δφ y and Δφ z in the equations 7b and 7c are used as the independent two phase differences, and the rest Δφ x in equation 7a, which is one of the above, can be set to zero (this corresponds to the fact that the initial phase of the x component is considered as a reference).
これらの量を用いると、式3は、位相差Δφx,Δφy,Δφzを用いて式8のように統一的に表現される。
When these quantities are used,
そして式8で各成分が規定される磁場ベクトルH(θ)#の大きさHn(θ)は、式9のようになる。ただし、量A,B,Cは式9a〜式9cで定義され、これらの式および以後の式での記号Σは「 k=x,y,z 」の和を示す。 The magnitude H n (θ) of the magnetic field vector H (θ) # in which each component is defined by Expression 8 is expressed by Expression 9. However, the quantities A, B, and C are defined by Equations 9a to 9c, and the symbol Σ in these equations and the following equations represents the sum of “k = x, y, z”.
この量Hn(θ)2の最大値は扁平楕円面Epの長軸Amaxの長さRmaxの2乗に相当し、量Hn(θ)2の最小値は短軸Aminの長さRminの2乗に相当する(図3参照)。 The maximum value of this quantity H n (θ) 2 is equivalent to the square of the length R max of the long axis A max of the flat ellipsoid Ep, the minimum value of the amount H n (θ) 2 is the minor axis A min long This corresponds to the square of R min (see FIG. 3).
この最大最小問題は解析的に求めることが可能であって、式10a〜式10cが得られる。 This maximum and minimum problem can be obtained analytically, and equations 10a to 10c are obtained.
また、磁場ベクトルが長軸Amaxに一致するとき、すなわち式11が成立するときのθをθmaxと書き、磁場ベクトルが短軸Aminに一致するとき、すなわち式12が成立するときのθをθminと書けは、これらは式13a〜式13dのように与えられる。 Further, θ when the magnetic field vector coincides with the long axis A max , that is, when Expression 11 is satisfied is written as θ max, and θ when the magnetic field vector coincides with the short axis A min , that is, when Expression 12 is satisfied. Can be written as θ min , these are given by Equations 13a to 13d.
そして、長軸Amaxの方向余弦umax(k) および短軸Aminの方向余弦umin(k) は、k=x,y,zにつき、式14aおよび式14bのように得られる。 The major axis A max direction cosine u max (k) and a minor axis A direction cosines of min u min (k) is, k = x, y, per z, obtained by the equation 14a and Formula 14b.
以上をまとめると、適当な時間幅にわたって磁場の3方向成分Hk(t)の値を測定することにより、これらの3成分のそれぞれの振幅hkと、2軸間の位相差Δφx ,Δφy ,Δφz (実質的にはそのうちの独立な2つ)とを求めれば、式9a〜式9cから式14a〜式14bまでの関係によって、扁平楕円面Epの長軸Amaxおよび短軸Aminの長さおよび方向を特定することができるということになる。 To summarize the above, by measuring the value of the three-direction component H k (t) of the magnetic field over an appropriate time width, the amplitude h k of each of these three components and the phase differences Δφ x , Δφ between the two axes. If y and Δφ z (substantially two of them) are obtained, the major axis A max and the minor axis A of the elliptical elliptical surface Ep are determined according to the relationship from the equations 9a to 9c to the equations 14a to 14b. This means that the length and direction of min can be specified.
そして、扁平楕円面Epの楕円率eは式15aの関係から定まり、傾き角αは方向余弦umax(k),umin(k)の値から特定できる(その詳細については後述する)。 The ellipticity e of the flat elliptical surface Ep is determined from the relationship of the expression 15a, and the inclination angle α can be specified from the values of the direction cosines u max (k) and u min (k) (details will be described later).
その結果、探査対象物の電気的比抵抗ρは、実験的あるいは経験的にあらかじめ特定されている既知関数fを用いて式15bから求めるか、あるいは2変数e,αをρに変換するためにあらかじめ作成しておいた数値テーブルから特定することができる。 As a result, the electrical resistivity ρ of the object to be searched is obtained from the expression 15b using a known function f specified experimentally or empirically in advance, or in order to convert the two variables e and α into ρ. It can be specified from a numerical table created in advance.
なお、以上のプロセス中で位相差Δφy、Δφzを求めるには、式17a,式17bの関係により、Hk(θ)の時間変動成分のうちのcos 成分ak とsin 成分bk とを求め、これらの変形である式18a,式18bと、三角関数の加法定理に相当する式18cを用いればよい。 In order to obtain the phase differences Δφ y and Δφ z in the above process, the cos component a k and the sin component b k of the time-varying components of H k (θ) are obtained from the relationship of the equations 17a and 17b. And these expressions, Expressions 18a and 18b, and Expression 18c corresponding to the addition theorem of trigonometric functions may be used.
ところで、上記の原理は、直交座標系(x,y,z)が大地(および地下構造)に対して固定された方向関係にあるという前提でのものであるが、ラジコンヘリ1は飛行中にその姿勢が変動することが多いから、その影響を考慮した補正を行うことが好ましい。 By the way, the above principle is based on the premise that the Cartesian coordinate system (x, y, z) is in a fixed directional relationship with the ground (and the underground structure). Since the posture often fluctuates, it is preferable to perform correction in consideration of the influence.
そこでまず、ラジコンヘリ1の測定部2に固定された座標系を局所座標系S(x,y,z)と呼び、大地に固定された座標系を絶対座標系S0(X,Y,Z)と呼んで、これらを相互に区別することにする。
First, the coordinate system fixed to the
後述するようにこの実施形態のシステムでは、測定部2に固定された直交3軸コイルを備えており、それによって局所座標系S(x,y,z)と絶対座標系S0(X,Y,Z)との角度関係を検出することができるようになっている。
As will be described later, the system of this embodiment includes an orthogonal three-axis coil fixed to the
このような座標関係を用いて局所座標系S(x,y,z)から絶対座標系S0(X,Y,Z)へのデータ変換を行う処理は、以下のように行うことができる。 The process of performing data conversion from the local coordinate system S (x, y, z) to the absolute coordinate system S 0 (X, Y, Z) using such a coordinate relationship can be performed as follows.
まず、局所座標系S(x,y,z)における長軸Amax方向の方向余弦ベクトルUmax(S)#と、短軸Amin方向の方向余弦ベクトルUmin(S)#とを、式14aおよび式14bの各成分を用いて式19a,式19bのように定義する。すると、3軸傾斜計の出力から得られる3行3列の座標変換マトリクスをM##としたとき、絶対座標系での長軸Amax方向および短軸Amin方向のそれぞれの方向余弦ベクトルUmax(0)#、Umin(0)#は、それぞれ式20a、式20bのように定まることになり、この姿勢補正によって、扁平楕円面Epの長軸Amaxの方向および短軸Aminの方向が、絶対座標系S0(X,Y,Z)で特定されることになる。 First, the local coordinate system S (x, y, z) and direction cosine vector U max (S) # of major axis A max direction of, the short axis A min direction of direction cosine vector U min (S) #, wherein Using the components of 14a and 14b, definitions are made as in equations 19a and 19b. Then, when a 3-by-3 coordinate conversion matrix obtained from the output of the 3-axis inclinometer is M ## , the respective direction cosine vectors U in the major axis A max direction and the minor axis A min direction in the absolute coordinate system. max (0) # and U min (0) # are determined as shown in Expression 20a and Expression 20b, respectively. By this posture correction, the direction of the major axis A max and the minor axis A min of the flat ellipsoid Ep The direction is specified by the absolute coordinate system S 0 (X, Y, Z).
絶対座標系S0(X,Y,Z)において、方向余弦ベクトルUmax(0)#,Umin(0)#から傾き角αを特定するには、たとえば以下のようにすることができる。 In the absolute coordinate system S 0 (X, Y, Z), the inclination angle α can be specified from the direction cosine vectors U max (0) # , U min (0) # , for example, as follows.
まず、図4に示すように、扁平楕円面EpとXY平面(水平面に相当)との交線上に単位ベクトルW#を定義する(以下これを「水平基準ベクトル」と呼ぶ)。方向余弦ベクトルUmax(0)#,Umin(0)#は、扁平楕円面Epが属する2次元傾斜平面上で互いに直交する基底ベクトルとなっているから、水平基準ベクトルW#は、定数係数rmax ,rminとこれらの方向余弦ベクトルUmax(0)#,Umin(0)#を用いて式21-1,式21-2のように表現できる。 First, as shown in FIG. 4, a unit vector W # is defined on the intersection line between the flat elliptical surface Ep and the XY plane (corresponding to a horizontal plane) (hereinafter referred to as “horizontal reference vector”). Since the direction cosine vectors U max (0) # and U min (0) # are basis vectors orthogonal to each other on the two-dimensional inclined plane to which the flat elliptical surface Ep belongs, the horizontal reference vector W # is a constant coefficient. Using r max and r min and their direction cosine vectors U max (0) # and U min (0) # , they can be expressed as shown in equations 21-1 and 21-2.
方向余弦ベクトルUmax(0)#,Umin(0)#を式21-3,式21-4のように成分表示し、絶対座標系S0(X,Y,Z)におけるZ方向の単位ベクトルZ#を式21-5のように成分表示すると、水平基準ベクトルW#がXY平面上にあることから、式21-6が成立する(式21-6中の記号「・」はベクトルの内積)。 The direction cosine vectors U max (0) # , U min (0) # are displayed as components as shown in Equation 21-3 and Equation 21-4, and the unit in the Z direction in the absolute coordinate system S 0 (X, Y, Z) When the vector Z # is displayed as a component as shown in Equation 21-5, since the horizontal reference vector W # is on the XY plane, Equation 21-6 is established (the symbol “•” in Equation 21-6 represents the vector inner product).
式21-1および式21-3〜式21-5を用いると、式21-6は式21-7となる。そして式21-2を考慮してこれを解くと、式21-8,式21-9が得られる。 Using Expression 21-1 and Expression 21-3 to Expression 21-5, Expression 21-6 becomes Expression 21-7. Solving this in consideration of Equation 21-2 yields Equation 21-8 and Equation 21-9.
一方、図4からわかるように、水平基準ベクトルW#と、長軸Amax方向の方向余弦ベクトルUmin(0)#との間の角度が、求めるべき傾き角度αであるから、これらのベクトルの内積として式22-1が成立する。これらのベクトルはいずれも単位ベクトルである(式22-2)。また、式21-1の両辺とベクトルW#との内積をとると式22-3が得られる。 On the other hand, as can be seen from FIG. 4, the angle between the horizontal reference vector W # and the direction cosine vector U min (0) # in the major axis A max direction is the inclination angle α to be obtained. Equation 22-1 is established as the inner product of. These vectors are all unit vectors (Equation 22-2). Further, when the inner product of both sides of Expression 21-1 and the vector W # is taken, Expression 22-3 is obtained.
式22-1〜式22-3から式22-4が得られるが、これが絶対座標系S0(X,Y,Z)における長軸Amaxの方向余弦ベクトルUmax(0)#のZ成分から傾き角αを求める式である。 Expression 22-4 is obtained from Expression 22-1 to Expression 22-3, which is the Z component of the direction cosine vector U max (0) # of the major axis A max in the absolute coordinate system S 0 (X, Y, Z). Is an expression for obtaining the inclination angle α from
この式22-4から得た傾き角αの値を式15bの傾き角αの値として使用することにより、比抵抗ρを得ることができる。 The resistivity ρ can be obtained by using the value of the inclination angle α obtained from the equation 22-4 as the value of the inclination angle α in the equation 15b.
なお、姿勢補正の必要がないような状況の場合には、図4に対応する上記の原理を局所座標系S(x,y,z)で得た方向余弦ベクトルそのものについて適用すればよい。 In the case where there is no need for posture correction, the above principle corresponding to FIG. 4 may be applied to the direction cosine vector itself obtained by the local coordinate system S (x, y, z).
<●4A.測定部2の構成と動作>
以上の原理に沿った処理を実行する具体的なシステム構成の説明に戻る。
<● 4A. Configuration and operation of measuring
Returning to the description of a specific system configuration for executing processing in accordance with the above principle.
図5および図6はそれぞれ、図1のシステムにおける測定部2および解析部3の内部の機能構成を示すブロック図である。
FIG. 5 and FIG. 6 are block diagrams showing functional configurations inside the
図5に示す測定部2は、直交3軸方向のそれぞれの磁場を測定する手段の感知部として3方向磁場センサ21を有している。この3方向磁場センサ21は、直交座標系の3軸方向に3つの単位コイル21x,21y,21zをそれぞれ配向させた直交3軸コイルCLを用いて構成されている。この直交3軸コイルCLの3軸方向は、既述した局所座標系S(x,y,z)で定義される方向である。
The measuring
3方向磁場センサ21の3つの出力信号は3チャネル増幅器22でそれぞれ増幅され、AD変換器23においてデジタル信号へと変換されるが、このAD変換器23は、AD変換トリガー信号生成回路25から与えられるトリガー信号に同期してAD変化を行う。
The three output signals of the three-direction
測定部2はまた、受信アンテナ24aを持つ電波時計24を備えている。電波時計24は図2のJJY標準電波に相当する送信波SWを受信するが、AD変換トリガー信号生成回路25は電波時計24の受信回路の出力から振幅変調開始時の立ち上がりを検出することによって、1secごとの単位変調期間の開始タイミングを検知し、そのタイミングに応答してAD変換器23にAD変換トリガー信号を出力する。AD変換器23はこのトリガー信号に応答して3チャネル増幅器22からの3入力それぞれのAD変換を開始し、所定間隔(たとえば1μsec)ごとにサンプリングとAD変換とを行う。AD変換器23は、AD変換開始から所定期間(たとえば128μsecまたは256μsec)経過すると一連のAD変換をいったん中止し、次のトリガー信号の入力を待つ待機状態となる。
The measuring
したがって、図7に示すように、3方向についての磁場の測定結果を反映した一連のデジタル信号(以下「3方向磁場データ列」)が1secおきに生成されるが、これらはトリガー信号に同期するタイミングで図5のデータ収録装置28に取り込まれ、その中に設けた記憶メディア28aにデジタル記録される。これは、既述した原理における各成分Hx(t),Hy(t),Hz(t)に相当する。
Therefore, as shown in FIG. 7, a series of digital signals (hereinafter referred to as “three-direction magnetic field data string”) reflecting the measurement results of the magnetic fields in the three directions are generated every 1 second, and these are synchronized with the trigger signal. It is taken into the
測定部2にはまた、3方向磁場センサ21に相対的に固定されてこの3方向磁場センサ21の傾きを検出するための3軸傾斜計26と、その時点での測定部2の経緯度位置を特定するためのGPSセンサ27とが設けられている。これらの3軸傾斜計26およびGPSセンサ27のそれぞれの出力はインターフェイス回路(I/F)を介してデータ収録装置28に入力される。このインターフェイス回路にもAD変換トリガー信号生成回路25からのトリガー信号が与えられており、このトリガー信号に応答して3軸傾斜計26からの傾斜情報信号とGPSセンサ27からの経緯度位置出力とがデータ収録装置28に取り込まれ、その時点で取り込まれた3方向磁場データ列に対応させて記憶メディア28aにデジタル記録される。図7中に、各単位変調期間ごとに姿勢情報やGPS位置情報が取り込まれることが模式的に示されている。
The measuring
このため、記憶メディア28aには、(3方向磁場データ列、傾斜情報、経緯度位置情報)が組となる形で、1secごとに記録されてゆくことになる。
For this reason, (three-direction magnetic field data string, gradient information, longitude and latitude position information) are recorded in a
GPSセンサ27の出力は3方向磁場データ列のそれぞれがどの経緯度位置で得られたデータであるかを特定するために使用され、3軸傾斜計26は絶対座標系S0(X,Y,Z)に対する局所座標系S(x,y,z)の傾きを検出するために使用される。後者は、既述した姿勢補正のために利用される。
The output of the
ラジコンヘリ1自身がその飛行制御のために電波時計やGPSセンサを備えているときには、それらを上記の電波時計24やGPSセンサ27として流用することもできる。また、ラジコンヘリ1と測定部2との相対的な姿勢関係が固定されており、かつラジコンヘリ1自身が飛行制御のために傾斜計を備えている場合には、それを3軸傾斜計26として流用してもよい。
When the radio controlled
また、ラジコンヘリ1の飛行記録(速度、高度など)をデータ収録装置28に入力させ、これらも記憶メディア28aにあわせて記録しておいてもよい。
Further, flight records (speed, altitude, etc.) of the radio controlled
ラジコンヘリ1を飛行させつつ以上の構成を用いた処理を繰り返して実行することにより、JJY標準信号における1secごとの周期的な振幅変調の開始に同期して、3方向磁場データ列とされた磁場波形が、局所座標の姿勢情報や経緯度位置とともに記憶メディア28aに記憶されてゆくことになる。
By repeatedly executing the process using the above configuration while flying the radio controlled
<●4B.解析部3の構成と動作>
図6の解析部3の機能はパーソナルコンピュータなどを用いてソフト的に実現することもできるが、図6には機能ブロックの形で表現されている。
<● 4B. Configuration and operation of
The functions of the
記憶メディア28aに記録された3方向磁場データ列などの情報は、オンラインまたはオフラインで解析部3に伝送される。
Information such as a three-direction magnetic field data string recorded in the
スペクトル解析信号抽出部31では、1単位変調期間分の3方向磁場データ列(すなわち128μsecまたは256μsecの磁場波形データ)をスペクトル解析して、この期間でのJJY標準電波の基準周波数(たとえば40Hz)成分を抽出する。そして、特徴抽出部32に含まれる振幅位相差特定部32aにおいて、3軸方向の磁場の基準周波数での振幅hkとそれぞれの2軸間の位相差Δφx,Δφy,Δφzを特定する。
The spectrum analysis
そして、扁平楕円特定部32bにおいては、3方向それぞれについての振幅hkと2軸位相差Δφx,Δφy,Δφzとに基づいて、扁平楕円面Epの特徴情報ないしは楕円パラメータとしての、長軸の長さRmax、短軸の長さRmin、および長軸短軸それぞれの方向余弦を決定する。これは、既述した式14までの算出原理に基づいて行われる。
Then, in the flat
また局所座標系S(x,y,z)と絶対座標系S0(X,Y,Z)との相対的な角度関係に基づいて補正を行う(式20a,式20b)。図6ではこの補正を行う機能手段として姿勢補正手段33が示されている。そして、これによって長軸短軸それぞれの方向余弦は、絶対座標系S0(X,Y,Z)で特定される。 Further, correction is performed based on the relative angular relationship between the local coordinate system S (x, y, z) and the absolute coordinate system S 0 (X, Y, Z) (Equation 20a, Equation 20b). In FIG. 6, the posture correcting means 33 is shown as a functional means for performing this correction. Thus, the direction cosine of each of the long axis and the short axis is specified by the absolute coordinate system S 0 (X, Y, Z).
3軸傾斜計26の出力を用いたこの補正処理は、特徴抽出部32で行われる各処理のどの段階で行ってもよい。すなわち、このような姿勢補正に対応する座標変換は線形変換であり、また扁平楕円面の決定も3次元線形空間内での処理であるから、座標変換処理と扁平楕円面Epの決定処理とは数学的に可換である。このため、局所座標系S(x,y,z)での3方向磁場データ列を座標変換マトリクスによって絶対座標系S0(X,Y,Z)での成分に変換した後に扁平楕円面Epを決定してもよく、また、局所座標系S(x,y,z)で扁平楕円面を決定した後に座標変換を行ってもよい。扁平楕円面Epの長軸の長さRmaxおよび短軸の長さRminはいずれもスカラー量であるから座標変換による影響は受けない。長軸短軸それぞれの方向余弦だけが座標変換によって変化する。
This correction process using the output of the
そして、これらが終わると式22-3,式22-4に基づいて傾き角αが決定され、また、扁平楕円の楕円率e(式15a)が特定される。式15bないしはこれに代わる評価式(または数値テーブル)によって、各測定場所Pi付近の地下構造の比抵抗ρが求まる。比抵抗ρは各測定場所Piごとに特定されるから、記号的には、測定場所Piに対応させて比抵抗ρiと書くことができる。 When these are finished, the inclination angle α is determined based on the equations 22-3 and 22-4, and the ellipticity e (formula 15a) of the flat ellipse is specified. The specific resistance ρ of the underground structure in the vicinity of each measurement location P i is obtained by the equation 15b or an evaluation equation (or a numerical table) instead. Since the specific resistance ρ is specified for each measurement location P i , symbolically, it can be written as the specific resistance ρ i corresponding to the measurement location P i .
それぞれの測定時におけるラジコンヘリ1の緯度経度位置(GPS27から得た位置情報)は記憶メディア28aから読出され、対応する各測定場所Piでの比抵抗ρiの値と関連づけて出力装置35(ディスプレイ、プリンタ、送信装置など)から出力される。
The latitude / longitude position (position information obtained from the GPS 27) of the radio controlled
また、測定処理を繰り返しながらラジコンヘリ1が測定エリアAR全体にわたる測線PL上を飛んで全測定位置での測定が完了すると、その解析によって比抵抗ρiの2次元的ないしは3次元的分布を地図と対応させてグラフィカルに表現することが可能になり、ディスプレイまたはプリンタによってそれを出力することができる。
In addition, when the
<●5.本方法の利点>
以上の方法においては、測定場所1地点についての測定は126μsecまたは256μsec(さらに長く測定時間をとるとしてもたとえば1024μsec程度)の期間でのサンプリングデータを用いて実行可能である一方、図2に示すようにどの単位変調期間も0.2sec(200000μsec)以上の長さを持つから、いずれの単位変調期間においても測定場所1地点分の測定を十分に行うことができる。
<● 5. Advantages of this method>
In the above method, measurement at one measurement location can be performed using sampling data in a period of 126 μsec or 256 μsec (for example, about 1024 μsec even if taking a longer measurement time), as shown in FIG. In addition, since any unit modulation period has a length of 0.2 sec (200000 μsec) or more, measurement for one measurement place can be sufficiently performed in any unit modulation period.
また、ラジコンヘリ1は約10m/sec程度の速さで飛行可能であるから、1secおきに測定を行うとすれば、水平距離で約10mおきに測定をすることが可能になる。時間的観点から見ると、たとえば延べ距離にして1000mの測線距離があっても、それに沿った約100点の測定が100sec(=約2分)程度で完了することになる。
In addition, since the
すなわち、この実施形態の方法においては、ミクロに見れば測定は離散的・間欠的に行われるが、上記のような空間的緻密さで測定を行って測定エリアARの地図上にその結果をマッピングすることにより、測定場所Piの全体が測定エリアARをほぼ覆い尽くすような分布になる。このため、マクロに見ると実質的に2次元的な「面走査」に近い測定が可能となる。 That is, in the method of this embodiment, the measurement is performed discretely or intermittently when viewed microscopically, but the measurement is performed with the spatial precision as described above and the result is mapped on the map of the measurement area AR. By doing so, the distribution is such that the entire measurement location Pi substantially covers the measurement area AR. For this reason, when viewed macroscopically, it is possible to measure substantially similar to a two-dimensional “surface scan”.
また、この発明の実施形態では、測定場所のひとつひとつについて測定磁場の扁平楕円面を特定し、比抵抗などの地下構造情報を得ることが可能である。従来のように地表を移動しつつ測定する場合は、既述した地形などの障害によって一つの測点場所での測定結果の信頼性が乏しいため、扁平楕円面の傾き方向を経験的に推定したり、いくつもの測定場所での測定結果を総合的に評価しなければ有意な情報が得られない場合も多い。これに対してこの発明の実施形態では測定場所のひとつひとつについて、高密度でかつ信頼性が高い測定結果を高速に得ることが可能になる。 Further, in the embodiment of the present invention, it is possible to specify a flat ellipsoidal surface of a measurement magnetic field for each measurement place and obtain underground structure information such as specific resistance. When measuring while moving on the ground as in the past, the reliability of the measurement results at one station location is poor due to the above-mentioned obstacles such as topography, so the inclination direction of the flat ellipsoid is estimated empirically. In many cases, significant information cannot be obtained without comprehensively evaluating the measurement results at several measurement locations. On the other hand, in the embodiment of the present invention, it is possible to obtain high-density and highly reliable measurement results at high speed for each measurement place.
<●6.変形例>
上記の実施形態では測定部3に電波時計を設け、その電波時計によってJJY電波の振幅変調開始タイミングをリアルタイムで検出し、それに基づいて測定磁場の1secごとのAD変換トリガーをかけているが、JJY標準電波によって事前に時刻調整された高精度の非電波時計を用いて1secごとのAD変換トリガーをかけてもよい。
<● 6. Modification>
In the above embodiment, the
測定磁場のデータは連続的にすべて取得して記憶しておき、扁平楕円面の決定に必要な1secごとのデータ部分だけを後から抽出して解析してもよい。 All of the measurement magnetic field data may be continuously acquired and stored, and only the data portion of every 1 second necessary for determining the oblate ellipsoid may be extracted and analyzed later.
独立な3方向の取り方や、独立な2つの位相差の取り方は任意である。たとえば円筒座標系で磁場測定を行ってもかまわない。 The method of taking three independent directions and the method of taking two independent phase differences are arbitrary. For example, magnetic field measurement may be performed in a cylindrical coordinate system.
飛行体としては、ラジコンヘリ1ではなく、有人のヘリコプターを使用することも可能である。
It is also possible to use a manned helicopter instead of the
100 電磁的地下探査システム
200 時刻報知用のJJY標準電波の発信局
1 ラジコンヘリ
2 測定部
3 解析部
Ep 扁平楕円面
Ec 扁平楕円面の外周
Amax 扁平楕円面の長軸
Amin 扁平楕円面の短軸
Rmax 扁平楕円面の長軸の長さ
Rmin 扁平楕円面の短軸の長さ
α 長軸の傾き角(伏角)
100 Electromagnetic
Claims (8)
前記電波として、時刻報知のために所定の時間周期ごとに所定の規則で振幅変調を行いつつ発信局から常時に公衆送信される送信波を利用し、
前記方法が、
独立な3方向についての磁場測定手段を保持した飛行体を測定対象エリアの上空を移動させつつ、前記3方向についての変動磁場を前記磁場測定手段によって測定する測定工程と、
前記時間周期ごとに生じる前記送信波の振幅変調期間のそれぞれにつき、前記変動磁場の測定の結果から、
1) 前記3方向の変動磁場のそれぞれの振幅の情報と、
2) 前記3方向の変動磁場の相対位相の情報と、
を求める特定工程と、
前記振幅の情報と前記相対位相の情報とから、前記3方向の変動磁場を成分とする3次元磁場ベクトルが3次元磁場空間において動く扁平楕円面の特徴情報を、前記振幅変調期間のそれぞれについて求める算出工程と、
を備え、
それぞれの振幅変調期間において前記飛行体が飛行していた各測定場所付近の地下構造についての情報を、各振幅変調期間に対応して得られた扁平楕円面の特徴情報に基づいて得ることを特徴とする電磁的地下探査方法。 A method of obtaining information on the underground structure by measuring a variable magnetic field corresponding to a combination of a primary magnetic field generated by a radio wave radiated from a predetermined transmitting station into the air and a secondary magnetic field generated by a response of the underground structure to the radio wave. There,
As the radio wave, using a transmission wave that is constantly transmitted from the source station to the public while performing amplitude modulation with a predetermined rule every predetermined time period for time notification,
Said method comprises
A measuring step of measuring a varying magnetic field in the three directions by the magnetic field measuring means while moving the flying object holding the magnetic field measuring means in three independent directions over the measurement target area;
For each amplitude modulation period of the transmitted wave that occurs in each time period, from the measurement result of the varying magnetic field,
1) Information on the amplitude of each of the three-direction fluctuating magnetic fields;
2) Information on the relative phase of the three-direction varying magnetic field;
Specific process for
From the amplitude information and the relative phase information, feature information of a flat ellipsoid in which a three-dimensional magnetic field vector having the three-direction varying magnetic field as a component moves in a three-dimensional magnetic field space is obtained for each of the amplitude modulation periods. A calculation process;
With
Obtaining information about the underground structure near each measurement location where the flying object flew during each amplitude modulation period based on the feature information of the flat ellipsoid obtained corresponding to each amplitude modulation period Electromagnetic underground exploration method.
前記測定工程が、
前記送信波を受信して、前記送信波のそれぞれの振幅変調期間の開始タイミングを検出するタイミング検出工程と、
それぞれの振幅変調期間の開始タイミングに同期して前記磁場測定手段の測定結果の取り込みを開始させる工程と、
を備えることと特徴とする電磁的地下探査方法。 The method of claim 1, wherein
The measuring step is
A timing detection step of receiving the transmission wave and detecting a start timing of each amplitude modulation period of the transmission wave;
Starting to capture measurement results of the magnetic field measurement means in synchronization with the start timing of each amplitude modulation period;
An electromagnetic underground exploration method characterized by comprising:
前記算出工程が、各測定場所での変動磁場測定結果に基づいてそれぞれ実行されることにより、各測定場所ごとに前記扁平楕円面の楕円軸に関する特徴情報が得られることを特徴とする電磁的地下探査方法。 The method of claim 1 or claim 2, wherein
The calculation process is performed based on the measurement results of the varying magnetic field at each measurement location, whereby characteristic information about the elliptical axis of the flat elliptical surface is obtained for each measurement location. Exploration method.
前記磁場測定手段に相対的に固定された姿勢関係で3軸傾斜計が前記飛行体に保持されており、
前記3軸傾斜計の検出結果に基づいて、前記磁場測定手段に相対的に固定された局所座標系と地上に固定された絶対座標系との姿勢関係を補正する工程、
を、前記算出工程に付随して実行することを特徴とする電磁的地下探査方法。 A method according to any one of claims 1 to 3,
A triaxial inclinometer is held on the flying body in a posture relationship relatively fixed to the magnetic field measuring means;
Correcting the posture relationship between the local coordinate system fixed relative to the magnetic field measuring means and the absolute coordinate system fixed on the ground based on the detection result of the three-axis inclinometer;
Is performed in association with the calculation step.
前記電波として、時刻報知のために所定の時間周期ごとに所定の規則で振幅変調を行いつつ発信局から常時に公衆送信される送信波を利用するものであり、
前記システムが、
独立な3方向についての変動磁場を空中において測定する磁場測定手段と、
前記磁場測定測定手段の出力から、前記振幅変調期間における波形を抽出する抽出手段と、
前記波形に基づいて
1) 前記3方向の変動磁場のそれぞれの振幅の情報と、
2) 前記3方向の変動磁場の相対位相の情報と、
を求める振幅位相特定手段と、
前記振幅の情報と前記相対位相の情報とから、前記3方向の変動磁場を成分とする3次元磁場ベクトルが3次元磁場空間において動く扁平楕円面の特徴情報を、前記振幅変調期間のそれぞれについて求める算出手段と、
を備え、
それぞれの振幅変調期間において前記飛行体が飛行していた各測定場所付近の地下構造についての情報を、各振幅変調期間に対応して得られた扁平楕円面の特徴情報に基づいて得ることを特徴とする電磁的地下探査システム。 Electromagnetic information for obtaining information on the underground structure by measuring a variable magnetic field corresponding to a combination of a primary magnetic field generated by a radio wave radiated from a predetermined transmitting station into the air and a secondary magnetic field generated by a response of the underground structure to the radio wave. An underground exploration system,
As the radio wave, a transmission wave that is always transmitted from the source station to the public while performing amplitude modulation according to a predetermined rule for a predetermined time period for time notification,
The system is
Magnetic field measuring means for measuring in the air a variable magnetic field in three independent directions;
Extraction means for extracting a waveform in the amplitude modulation period from the output of the magnetic field measurement measurement means;
Based on the waveform
1) Information on the amplitude of each of the three-direction fluctuating magnetic fields;
2) Information on the relative phase of the three-direction varying magnetic field;
Amplitude phase identification means for obtaining
From the amplitude information and the relative phase information, feature information of a flat ellipsoid in which a three-dimensional magnetic field vector having the three-direction varying magnetic field as a component moves in a three-dimensional magnetic field space is obtained for each of the amplitude modulation periods. A calculation means;
With
Obtaining information about the underground structure near each measurement location where the flying object flew during each amplitude modulation period based on the feature information of the flat ellipsoid obtained corresponding to each amplitude modulation period Electromagnetic underground exploration system.
前記送信波を受信して、前記時間周期ごとに生じる前記送信波の振幅変調期間を検知する検知手段、
をさらに備え、
前記抽出手段は、前記検知手段によって検知された振幅変調期間から前記波形を抽出することを特徴とする電磁的地下探査システム。 The system of claim 5, wherein
Detecting means for receiving the transmission wave and detecting an amplitude modulation period of the transmission wave generated at each time period;
Further comprising
The electromagnetic underground exploration system, wherein the extraction means extracts the waveform from an amplitude modulation period detected by the detection means.
前記算出手段が、各測定場所での変動磁場測定結果に基づいて、各測定場所ごとに前記扁平楕円面の楕円軸に関する特徴情報を得ることを特徴とする電磁的地下探査システム。 The system of claim 5 or claim 6,
The electromagnetic subsurface exploration system, wherein the calculation means obtains characteristic information about the elliptical axis of the flat ellipsoid for each measurement location based on the measurement result of the varying magnetic field at each measurement location.
前記磁場測定手段に相対的に固定された姿勢関係で3軸傾斜計が前記飛行体に保持されており、
前記3軸傾斜計の検出結果に基づいて、前記磁場測定手段に相対的に固定された局所座標系と地上に固定された絶対座標系との姿勢関係を補正する補正手段、
をさらに備えることを特徴とする電磁的地下探査システム。 The system according to any one of claims 5 to 7,
A triaxial inclinometer is held on the flying body in a posture relationship relatively fixed to the magnetic field measuring means;
Correction means for correcting the posture relationship between the local coordinate system fixed relative to the magnetic field measurement means and the absolute coordinate system fixed on the ground based on the detection result of the three-axis inclinometer;
An electromagnetic underground exploration system characterized by further comprising:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006056226A JP2007232621A (en) | 2006-03-02 | 2006-03-02 | Electromagnetic underground prospecting method and system for the same |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006056226A JP2007232621A (en) | 2006-03-02 | 2006-03-02 | Electromagnetic underground prospecting method and system for the same |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2007232621A true JP2007232621A (en) | 2007-09-13 |
Family
ID=38553333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006056226A Pending JP2007232621A (en) | 2006-03-02 | 2006-03-02 | Electromagnetic underground prospecting method and system for the same |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2007232621A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015088205A1 (en) * | 2013-12-09 | 2015-06-18 | 한국지질자원연구원 | 3d aeromagnetic surveying system and 3d aeromagnetic surveying method using same |
WO2016017245A1 (en) * | 2014-07-31 | 2016-02-04 | ソニー株式会社 | Information processing device, information processing method, and image display system |
JP2018194454A (en) * | 2017-05-18 | 2018-12-06 | 国立研究開発法人産業技術総合研究所 | Underground electromagnetic survey device |
-
2006
- 2006-03-02 JP JP2006056226A patent/JP2007232621A/en active Pending
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015088205A1 (en) * | 2013-12-09 | 2015-06-18 | 한국지질자원연구원 | 3d aeromagnetic surveying system and 3d aeromagnetic surveying method using same |
US10185048B2 (en) | 2013-12-09 | 2019-01-22 | Korea Institute Of Geoscience And Mineral Resources | 3-dimensional airborne magnetic survey system and 3-dimensional airborne magnetic survey method using the same |
WO2016017245A1 (en) * | 2014-07-31 | 2016-02-04 | ソニー株式会社 | Information processing device, information processing method, and image display system |
JPWO2016017245A1 (en) * | 2014-07-31 | 2017-04-27 | ソニー株式会社 | Information processing apparatus, information processing method, and image display system |
US20170278262A1 (en) | 2014-07-31 | 2017-09-28 | Sony Corporation | Information processing device, method of information processing, and image display system |
US10269132B2 (en) | 2014-07-31 | 2019-04-23 | Sony Corporation | Displaying images according to head posture and camera posture |
JP2018194454A (en) * | 2017-05-18 | 2018-12-06 | 国立研究開発法人産業技術総合研究所 | Underground electromagnetic survey device |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2269086B1 (en) | Positioning, detection and communication system and method | |
US9151822B2 (en) | Precise positioning using a distributed sensor network | |
US20110066379A1 (en) | survey system for locating geophysical anomalies | |
US8018382B2 (en) | Positioning system and method | |
EP3472650A1 (en) | Systems and methods for uniquely identifying buried utilities in a multi-utility environment | |
US10048073B2 (en) | Beacon-based geolocation using a low frequency electromagnetic field | |
Li et al. | Using geomagnetic field for indoor positioning | |
US20100026304A1 (en) | Method and Apparatus for Analysing Geological Features | |
US20190011592A1 (en) | Tracked distance measuring devices, systems, and methods | |
CN102918426B (en) | Use being accurately positioned for distributed sensor networks | |
JP2005098746A (en) | Continuous cable position survey device, method, and program | |
WO2015085426A1 (en) | Electromagnetic surveying at low frequencies using an airborne transmitter with receivers on the ground | |
JP2007232621A (en) | Electromagnetic underground prospecting method and system for the same | |
CN112505787B (en) | Coal seam roof water electromagnetic method perspective exploration system and method | |
CN102252674A (en) | Underwater geomagnetic positioning and navigation device | |
US20230176233A1 (en) | Gnss positioning methods and devices using ppp-rtk, rtk, ssr, or like correction data | |
JP2018194454A (en) | Underground electromagnetic survey device | |
JP3740539B2 (en) | Electromagnetic exploration apparatus and method | |
CN109782222A (en) | A kind of the underground utilities localization method and system of integrated RTLS | |
CN217467206U (en) | Three-component measurement system for magnetic source ground-air transient electromagnetism | |
JP2009276254A (en) | Buried object locating system | |
JP4059516B1 (en) | Electromagnetic wave measurement system | |
EP4185898A2 (en) | Vehicle-based utility locating using principal components | |
Helle et al. | Indoor Localization Solutions for a Marine Industry Augmented Reality Tool | |
CN114002747A (en) | Geologic body resistivity detection method and system for ground transmission and air reception |