JP4633565B2 - River data measurement method and apparatus - Google Patents

River data measurement method and apparatus Download PDF

Info

Publication number
JP4633565B2
JP4633565B2 JP2005208406A JP2005208406A JP4633565B2 JP 4633565 B2 JP4633565 B2 JP 4633565B2 JP 2005208406 A JP2005208406 A JP 2005208406A JP 2005208406 A JP2005208406 A JP 2005208406A JP 4633565 B2 JP4633565 B2 JP 4633565B2
Authority
JP
Japan
Prior art keywords
ultrasonic
river
reflected
data
intensity
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 - Fee Related
Application number
JP2005208406A
Other languages
Japanese (ja)
Other versions
JP2007024715A (en
Inventor
由博 道口
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2005208406A priority Critical patent/JP4633565B2/en
Publication of JP2007024715A publication Critical patent/JP2007024715A/en
Application granted granted Critical
Publication of JP4633565B2 publication Critical patent/JP4633565B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Volume Flow (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Description

本発明は、超音波を用い河川などの川底形状(河床形状)や流量などの河川データを測定する方法及び装置に係り、特に、測定装置の移動量と移動方向の測定精度の向上に関する。   The present invention relates to a method and apparatus for measuring river data such as a riverbed shape (river bed shape) and a flow rate of a river using ultrasonic waves, and more particularly to improvement in measurement accuracy of a moving amount and a moving direction of a measuring device.

河川の河床形状は水面から河底までの距離、つまり、水深を河川断面で連続的に求めることで測定できる。流速分布は色々な深さを流れる反射体の移動速度から計測できる。流量はある位置での水深と水中の流速分布を検出すれば、河川断面の積分量として換算できる。これらの計測には超音波を用いることが多い。   The riverbed shape of a river can be measured by continuously obtaining the distance from the water surface to the riverbed, that is, the water depth in the river cross section. The flow velocity distribution can be measured from the moving speed of the reflector flowing through various depths. The flow rate can be converted as an integral value of the river cross section by detecting the water depth at a certain position and the flow velocity distribution in the water. Ultrasonic waves are often used for these measurements.

超音波を用いる検出方法は、例えば、特許文献1の「水深・流速・水温測定装置」に述べられている。この方法は、深さ方向に棒状の垂下材を設け、その垂下材に沿って複数の超音波トランスジューサを設置する。河底の超音波トランスジューサから超音波を発射し、垂下材各所のトランスジューサでこれを受信し、その伝搬時間を求めて各トランスジューサ間の音速・水深を測定するものである。この方法は、河川の断面全体の河床形状、流速分布を求めるためには、垂下材を多数設ける必要がある。   A detection method using ultrasonic waves is described, for example, in “Water depth / flow velocity / water temperature measuring device” of Patent Document 1. In this method, a rod-like hanging material is provided in the depth direction, and a plurality of ultrasonic transducers are installed along the hanging material. Ultrasonic waves are emitted from ultrasonic transducers at the bottom of the river, received by the transducers at various locations of the drooping material, and the propagation time is obtained to measure the speed of sound and the depth of water between the transducers. In this method, in order to obtain the riverbed shape and flow velocity distribution of the entire cross section of the river, it is necessary to provide many drooping materials.

超音波トランスジューサを移動させるため移動手段を備えたものに、特許文献2の「河川流量測定装置及び方法」がある。この方法も、超音波トランスジューサを複数個水面及び水中に配置する必要があり、さらに、河底付近に超音波反射体を沈める必要があるため、簡便な測定法ではない。   Patent Document 2 “River Flow Rate Measuring Device and Method” includes a moving means for moving an ultrasonic transducer. This method is also not a simple measurement method because it is necessary to place a plurality of ultrasonic transducers on the water surface and in the water, and it is necessary to sink an ultrasonic reflector near the riverbed.

特許文献3に記載された「河川状態計測方法および装置」の方法は、水のない部分はレーザー光線計測装置で形状を測定し、河川部分の形状計測には超音波を用いるものである。河床形状の計測に用いる超音波は100kHz〜500kHzを使っているが、超音波ビームの広がりによる精度低下がある。また、流速分布を測定することができない。   In the method of “river state measuring method and apparatus” described in Patent Document 3, the shape of a portion without water is measured with a laser beam measuring device, and ultrasonic waves are used for measuring the shape of the river portion. The ultrasonic wave used for measuring the riverbed shape is 100 kHz to 500 kHz, but there is a decrease in accuracy due to the spread of the ultrasonic beam. In addition, the flow velocity distribution cannot be measured.

複数の超音波ビームを使い海洋での水深、流速分布を測定し、流量を求める装置で、さらに計測器を搭載した水面上の浮揚体位置を求める公知例として特許文献4の方法がある。この方法は、例えば4本の超音波を垂直方向から設定角度離して斜め方向に、発信源から離れるにつれて互いに離れる方向に送信し、水中散乱体の反射を検出して流速を求める。さらに、海底の反射から深さを求めるとともに、測定装置を搭載した浮揚体の移動を検出する。   There is a method disclosed in Patent Document 4 as a known example for obtaining the position of a floating body on the water surface, which is a device for measuring the water depth and flow velocity distribution in the ocean using a plurality of ultrasonic beams and obtaining the flow rate, and further equipped with a measuring instrument. In this method, for example, four ultrasonic waves are transmitted obliquely at a set angle away from the vertical direction and in directions away from each other as they move away from the transmission source, and the reflection of underwater scatterers is detected to determine the flow velocity. Furthermore, the depth is obtained from the reflection of the seabed, and the movement of the floating body equipped with the measuring device is detected.

その方法は、複数回超音波を放射して反射波を検出し、その反射波検出時間の変化を求める。海底は静止しているので、超音波を照射した部分のもっとも反射強度が高い位置は不動と仮定すると、浮揚体の移動に伴って反射波検出時間は変化する。この変化量から不動である海底から相対的に浮揚体が移動した量と方向を求めるものである。   In this method, a reflected wave is detected by emitting an ultrasonic wave a plurality of times, and a change in the reflected wave detection time is obtained. Since the seabed is stationary, assuming that the position where the reflection intensity is highest in the portion irradiated with ultrasonic waves is immobile, the reflected wave detection time changes as the floating body moves. From this amount of change, the amount and direction of the relative movement of the buoyant body from the stationary seabed is obtained.

この方法では、ビーム間に角度があるため、水深が大きくなるほど超音波ビームは互いに離れて行き、離れた位置での海底状態が異なると検出時間差に誤差がでる。また、海底に凹凸があれば異なる超音波往復時間となって、大きな誤差を伴い、移動量の測定が不可能となる。海底に比べ川底の凹凸が激しい河川への適用はさらに難しい。   In this method, since there is an angle between the beams, the ultrasonic beams move away from each other as the water depth increases, and an error occurs in the detection time difference if the seabed state at the distant position is different. In addition, if there is unevenness on the seabed, the ultrasonic round-trip time will be different, which will cause a large error and make it impossible to measure the amount of movement. It is even more difficult to apply to rivers where the unevenness of the riverbed is more severe than the seabed.

超音波を使用せずに河床形状を求めるものに、特許文献5に記載された「河床探査具」がある。これは、河底から採取した河床材料に永久磁石を埋め込んだ探査具を多数河底に配置し、磁場の測定からこの探査具の位置や移動を測るものである。この方法は、測定範囲が限られること、水深測定に別の手段が必要なことなどの課題がある。   There exists a "bed search tool" described in patent document 5 as what calculates | requires a riverbed shape, without using an ultrasonic wave. In this method, a number of probe tools with permanent magnets embedded in riverbed material collected from the riverbed are placed on the riverbed, and the position and movement of the probe are measured from the measurement of the magnetic field. This method has problems such as a limited measurement range and the need for another means for water depth measurement.

特開2000−97738号公報(第6,7頁、図1)Japanese Unexamined Patent Publication No. 2000-97738 (6th and 7th pages, FIG. 1) 特許第2955920号公報(第8〜9頁、図1)Japanese Patent No. 2955920 (pages 8-9, FIG. 1) 特開平11−304484号公報(第3頁、図1)Japanese Patent Laid-Open No. 11-304484 (page 3, FIG. 1) USP5122990USP 5122990 特開2000−98046号公報(第3頁、図2)JP 2000-98046 A (page 3, FIG. 2) 応用物理学会誌、第46巻、第8号(1977年)、「スペックルの性質と応用」Journal of Japan Society of Applied Physics, Vol. 46, No. 8 (1977), “Properties and Applications of Speckle”

従来、浮揚体に搭載した発信源から発信される、中心軸線の異なる複数の超音波ビームを使って河川等の水深、流速分布を測定して、流量を求める装置においては、
1)河底の凹凸により浮揚体移動中に各超音波ビームの行程が変化し、河底での反射位置の水平方向移動量が、浮揚体の水平方向移動量と一致しない場合が生じるため移動量・移動方向の測定に大きな誤差が発生し、河床形状の測定精度が低下する、
2)ビーム間に角度があるため水深が大きくなるほど超音波ビームは離れて行き、異なる場所の流速を測ることになり流速、ひいては流量の誤差となる、
という惧れがあった。
Conventionally, in devices that measure the water depth and flow velocity distribution of rivers etc. using multiple ultrasonic beams with different central axes transmitted from the transmission source mounted on the levitation body,
1) The movement of each ultrasonic beam changes during the movement of the floating body due to the unevenness of the riverbed, and the amount of horizontal movement of the reflection position at the riverbed may not match the horizontal movement of the floating body. A large error occurs in the measurement in the quantity / movement direction, and the measurement accuracy of the riverbed shape decreases.
2) Since there is an angle between the beams, the ultrasonic beam moves away as the water depth increases, and the flow velocity at different locations is measured, resulting in an error in the flow velocity, and hence the flow rate.
There was a concern.

すなわち、本発明の第1の課題は、測定装置の移動量の計測精度低下に伴う河床形状の測定精度の低下を防ぐことであり、第2の課題は、流速分布、流量の測定精度の低下を防ぐことである。   That is, the first problem of the present invention is to prevent a decrease in the measurement accuracy of the riverbed shape due to a decrease in the measurement accuracy of the movement amount of the measuring device, and the second problem is a decrease in the measurement accuracy of the flow velocity distribution and the flow rate. Is to prevent.

本発明の第1の課題は、測定装置の移動量、移動方向の測定精度を向上させることにより解決される。具体的には、本発明の第1の課題は、測定線に沿って移動する浮揚体に搭載され、河底に向かって発信される超音波によって河川のデータを測定する測定装置であって、河底方向に所定の時間間隔で超音波を発信する超音波送信手段と、前記超音波送信手段から発信され、河底で散乱・反射された超音波を受信する2次元状に配置された2次元超音波受信手段と、前記超音波送信手段から時間をおいて発信された2回の超音波それぞれに対応する前記2次元超音波受信手段の出力に基づいて河底からの反射波の強度分布を前記2次元超音波受信手段の中心を基準位置として求める信号処理手段と、信号処理手段に接続され、信号処理手段から出力される二つの強度分布の間の相関の強度を、前記二つの強度分布のうちの先に発信された超音波に対応するものを他方に近づけつつ移動距離の関数として求め、相関の強度が最大となる移動距離及びそのときの移動方向を前記2回の超音波発信の間の超音波送信手段の移動ベクトルとして出力する信号解析手段と、前記移動ベクトルが格納される記憶手段とを有してなる河川データ測定装置により、解決される。 The first problem of the present invention is solved by improving the measurement accuracy of the moving amount and moving direction of the measuring apparatus. Specifically, a first problem of the present invention is a measuring device that is mounted on a floating body that moves along a measurement line, and that measures river data by ultrasonic waves transmitted toward the river bed, An ultrasonic transmission means for transmitting ultrasonic waves at a predetermined time interval in the riverbed direction, and two-dimensionally arranged to receive ultrasonic waves transmitted from the ultrasonic transmission means and scattered and reflected by the riverbed Based on the output of the two-dimensional ultrasonic wave reception unit corresponding to each of the two ultrasonic waves transmitted from the ultrasonic wave transmission unit with a time interval from the two-dimensional ultrasonic wave reception unit. and determined Mel signal processing means as a reference position the center of the fabric the 2-dimensional ultrasound receiving means is connected to the signal processing means, the strength of the correlation between the two intensity distribution output from the signal processing unit , ultrasound previously originated in one of the two intensity distribution The corresponding one is obtained as a function of the moving distance while approaching the other, and the moving distance at which the correlation strength is maximum and the moving direction at that time are output as the movement vector of the ultrasonic transmission means between the two ultrasonic transmissions. This is solved by a river data measuring device having signal analyzing means for storing and storage means for storing the movement vector.

上記構成によれば、超音波送信手段から河底に向け、一定周波数で位相のそろった超音波が発信される。発信された超音波は、河底の照射領域で反射される。照射領域の表面が超音波の波長以上のランダムな凹凸があると、超音波反射は散乱状態となる。通常、河底は砂や小石があるため、その凹凸は波長以上になり種々の方向に超音波が反射される散乱状態となる。散乱反射波は、空間的に互いに干渉しあい、超音波の重ね合わせで強め合う部分と打ち消しあう部分が発生し、強度の高い場所と低い場所をもった強度分布を形成する。この強度むらを形成する現象をスペックルといい、波長がそろったレーザー光や超音波で発生する。この強度むらであるスペックルパターンが、水面付近で2次元状に配置された2次元超音波受信手段で受信される。超音波送信手段と2次元超音波受信手段は浮揚体に搭載されているから、超音波は浮揚体と一体となって移動する。移動前後のスペックルパターンは、移動量が小さく超音波照射領域の重なりが大きいと、移動前後の強度分布の変化は少ない。スペックルパターンの強度分布を移動前後で求めると、移動により検出位置が動いた分、強度分布がシフトするため、移動前後の強度分布の相関を計算すると、相関のピーク位置から移動量と移動方向を求めることができる。   According to the above configuration, ultrasonic waves having a constant frequency are transmitted from the ultrasonic transmission means toward the riverbed. The transmitted ultrasonic wave is reflected by the irradiation area of the riverbed. If the surface of the irradiated region has random irregularities that are longer than the wavelength of the ultrasonic wave, the ultrasonic reflection is in a scattering state. Usually, the riverbed has sand and pebbles, so that the unevenness becomes longer than the wavelength, and it becomes a scattering state in which ultrasonic waves are reflected in various directions. The scattered reflected waves spatially interfere with each other, and a portion that reinforces and cancels out by superposition of ultrasonic waves is generated, and an intensity distribution having a high intensity place and a low place is formed. This phenomenon of unevenness in intensity is called speckle, and is generated by laser light or ultrasonic waves with the same wavelength. The speckle pattern, which is uneven in intensity, is received by a two-dimensional ultrasonic receiving means arranged two-dimensionally near the water surface. Since the ultrasonic transmission means and the two-dimensional ultrasonic reception means are mounted on the levitation body, the ultrasonic waves move together with the levitation body. The speckle pattern before and after the movement has a small change in the intensity distribution before and after the movement when the movement amount is small and the overlap of the ultrasonic irradiation regions is large. When the intensity distribution of the speckle pattern is calculated before and after movement, the intensity distribution shifts by the amount of movement of the detection position due to movement. Therefore, if the correlation of the intensity distribution before and after movement is calculated, the movement amount and movement direction from the correlation peak position. Can be requested.

このように、超音波を用いそのスペックルパターンの移動によるシフトから、超音波送信、受信部分、つまり、浮揚体の移動を精度よく検出できる。   In this way, it is possible to accurately detect the movement of the ultrasonic transmission and reception portion, that is, the floating body, from the shift due to the movement of the speckle pattern using ultrasonic waves.

前記超音波送信手段に代えて、河底方向に所定の時間間隔で超音波を発信するとともに反射波を受信する超音波送受信手段と、前記超音波送受信手段に接続され、超音波発信から河底の反射波受信までの時間と河底からの反射波の強度を測定する河底反射測定手段とが設けられ、前記信号解析手段は、前記河底反射測定手段の出力を入力として水深を算出して前記記憶手段に格納し、さらに、前記移動ベクトルを順次結合して前記超音波送受信手段の移動の軌跡を生成するとともに、前記移動の軌跡を構成する各移動ベクトルの終点位置に対応する水深のデータを前記記憶手段から読み出して前記終点位置のデータに組み合わせて前記記憶手段に格納するよう構成されていることを特徴とする河川データ測定装置としてもよい。   In place of the ultrasonic transmission means, ultrasonic transmission / reception means for transmitting ultrasonic waves at a predetermined time interval in the river bed direction and receiving reflected waves, and connected to the ultrasonic transmission / reception means, from the ultrasonic transmission to the river bed A bottom reflection measuring means for measuring the time until reception of the reflected wave and the intensity of the reflected wave from the bottom of the river is provided, and the signal analysis means calculates the water depth using the output of the bottom reflection measuring means as an input. Are stored in the storage means, and the movement vectors are sequentially combined to generate a movement trajectory of the ultrasonic transmission / reception means, and the water depth corresponding to the end point position of each movement vector constituting the movement trajectory is A river data measuring device may be configured to read data from the storage unit and store the data in the storage unit in combination with the data of the end point position.

このように構成することで、各移動ベクトルの終点位置における水深を、前記終点位置のデータに容易に組み合わせることができ、河床形状を精度よく測定することができる。   By comprising in this way, the water depth in the end point position of each movement vector can be easily combined with the data of the said end point position, and a riverbed shape can be measured accurately.

本発明の前記第2の課題は、前記河川データ測定装置において、前記超音波送信手段に代えて、河底方向に所定の時間間隔で超音波を発信するとともに反射波を受信する超音波送受信手段と、前記超音波送受信手段に接続され、超音波発信から河底の反射波受信までの時間と河底からの反射波の強度を測定する河底反射測定手段とが設けられ、前記2次元超音波受信手段は河底で散乱・反射された超音波に加え、水中に浮遊する水中反射体で反射された反射波を受信するよう構成されているとともに前記信号処理手段は、超音波発信から水中反射体で反射された反射波を検出するまでの時間と該反射波の強度分布及び位相分布をも測定するよう構成され、前記記憶手段には、それぞれ異なる深さでかつ前記超音波送受信手段との水平方向相対位置が所定の位置にある水中反射体によって前記超音波送受信手段から発信された超音波が反射されたときの反射波の強度分布及び位相分布が基準強度分布及び基準位相分布として格納され、前記信号解析手段は、前記河底反射測定手段の出力を入力として水深を算出して前記記憶手段に格納し、前記信号処理手段の出力に基づいて水中反射体が浮遊している深さを検出し、前記検出された深さの水中反射体の反射波の基準強度分布及び基準位相分布と前記測定された反射波の強度分布及び位相分布のいずれか若しくは双方の相関の強度を、測定された分布を基準分布に近づけつつ移動距離の関数として求め、相関の強度が最大となる移動距離及びそのときの移動方向に基づいて当該水中反射体の超音波送受信手段に対する水平方向位置を求め、所定の時間間隔で検出された水中反射体の超音波送受信手段に対する水平方向相対位置及び前記算出された超音波送受信手段の移動ベクトルに基づいて、前記移動ベクトルの終点位置における、前記検出された深さの流速を求めるように構成されていることを特徴とする河川データ測定装置により解決される。   The second problem of the present invention is that, in the river data measuring apparatus, in place of the ultrasonic transmission means, ultrasonic transmission / reception means for transmitting ultrasonic waves at a predetermined time interval and receiving reflected waves in the river bottom direction. And a bottom reflection measuring means connected to the ultrasonic transmission / reception means for measuring the time from transmission of ultrasonic waves to reception of reflected waves at the bottom of the river and the intensity of reflected waves from the bottom of the river. The sound wave receiving means is configured to receive the reflected wave reflected by the underwater reflector floating in the water in addition to the ultrasonic waves scattered and reflected at the riverbed, and the signal processing means is configured to receive the ultrasonic wave from the underwater transmission. The time until the reflected wave reflected by the reflector is detected and the intensity distribution and phase distribution of the reflected wave are also measured, and the storage means includes the ultrasonic transmission / reception means at different depths. Horizontal relative to The intensity distribution and phase distribution of the reflected wave when the ultrasonic wave transmitted from the ultrasonic transmitting / receiving means is reflected by the underwater reflector at a predetermined position are stored as the reference intensity distribution and the reference phase distribution, and the signal The analysis means calculates the water depth using the output of the river bed reflection measurement means as an input and stores it in the storage means, and detects the depth at which the underwater reflector floats based on the output of the signal processing means, The reference intensity distribution and reference phase distribution of the reflected wave of the underwater reflector of the detected depth and the intensity of the correlation between one or both of the intensity distribution and phase distribution of the measured reflected wave are measured distribution. The horizontal position of the underwater reflector with respect to the ultrasonic transmission / reception means is obtained based on the movement distance at which the correlation intensity is maximized and the movement direction at that time. The detected position at the end position of the movement vector based on the horizontal relative position of the underwater reflector detected with respect to the ultrasonic transmission / reception means and the calculated movement vector of the ultrasonic transmission / reception means detected at predetermined time intervals. Solved by a river data measuring device, characterized in that it is configured to determine the flow velocity of the depth.

上記構成によれば、浮揚体に搭載された超音波送受信手段から河底に向け、一定周波数で位相のそろった超音波が発信される。超音波は水中を伝搬するが、その水中には水中を浮遊しつつ水流とともに移動する反射体が多数存在する。発信された超音波ビーム中に前記反射体、すなわち水中反射体があると、超音波の一部はその水中反射体で反射され、水中反射体を中心とする球面波として上方に伝搬する。この反射球面波を2次元状に配置した2次元超音波受信手段で受信して2次元のデータとし、2次元超音波受信手段に接続された信号処理手段で前記2次元のデータに基づいて反射波の強度分布や位相分布を2次元のデータとして測定する。超音波を所定の時間間隔で複数回発信し、水中反射体の位置が連続的に動いていると球面波の2次元超音波受信手段に対する分布状態も移動し、この移動シフト量を検出することで水中反射体の移動、つまり、ある深さでの流速を検出できる。流速を検出する深さは、超音波送信からの時間で弁別し、深さ方向の流速分布を知ることができる。   According to the said structure, the ultrasonic wave with which the phase was uniform was transmitted to the riverbed from the ultrasonic transmission / reception means mounted in the floating body. Ultrasound propagates in water, but there are many reflectors that move with the water flow while floating in the water. If the reflector, that is, the underwater reflector, is present in the transmitted ultrasonic beam, a part of the ultrasonic wave is reflected by the underwater reflector and propagates upward as a spherical wave centered on the underwater reflector. The reflected spherical wave is received by a two-dimensional ultrasonic receiving means arranged two-dimensionally to obtain two-dimensional data, and reflected based on the two-dimensional data by a signal processing means connected to the two-dimensional ultrasonic receiving means. Wave intensity distribution and phase distribution are measured as two-dimensional data. When ultrasonic waves are transmitted a plurality of times at predetermined time intervals, and the position of the underwater reflector continuously moves, the distribution state of the spherical wave with respect to the two-dimensional ultrasonic wave receiving means also moves, and this movement shift amount is detected. The movement of the underwater reflector, that is, the flow velocity at a certain depth can be detected. The depth at which the flow velocity is detected can be discriminated by the time from ultrasonic transmission, and the flow velocity distribution in the depth direction can be known.

前記移動ベクトル終点位置となる超音波発信時に対応する深さ方向の流速分布を、各移動ベクトル終点位置となる超音波発信時について積分することにより、流量を得ることができる。   The flow rate can be obtained by integrating the flow velocity distribution in the depth direction corresponding to the time when ultrasonic waves are transmitted at the movement vector end positions with respect to the time when ultrasonic waves are transmitted at the movement vector end positions.

本発明の前記第1の課題は、測定線に沿って移動する浮揚体に搭載された測定装置から河底に向かって発信される超音波によって河川のデータを測定する河川データ測定方法であって、河底方向に所定の時間間隔で超音波を発信し、河底で散乱・反射された超音波を2次元状に配置された2次元超音波受信手段で受信し、時間をおいて発信された2回の超音波それぞれに対応する河底からの反射波の強度分布を前記2次元超音波受信手段の中心を基準位置としてそれぞれ求め、前記二つの強度分布の間の相関の強度を、前記二つの強度分布のうちの先に発信された超音波に対応するものを他方に近づけつつ移動距離の関数として求め、相関の強度が最大となる移動距離及びそのときの移動方向を前記2回の超音波発信の間の超音波発信源の移動距離及び移動方向として出力する手順を有してなる河川データ測定方法によっても、解決される。 The first problem of the present invention is a river data measurement method for measuring river data by ultrasonic waves transmitted toward a river bed from a measurement device mounted on a levitating body that moves along a measurement line. The ultrasonic waves are transmitted at predetermined time intervals in the riverbed direction, and the ultrasonic waves scattered / reflected at the riverbed are received by the two-dimensional ultrasonic receiving means arranged in a two-dimensional form, and transmitted at a time. twice the intensity distribution of the reflected wave from corresponding riverbed ultrasonic each respectively determined center of the two-dimensional ultrasonic receiving means as a reference position and the intensity of correlation between the two intensity distribution the movement direction of the two the one corresponding to the ultrasonic wave transmitted earlier of the intensity distribution while close to the other determined as a function of the moving distance, the moving distance and the time the strength of the correlation is the maximum Transfer of ultrasonic source between two ultrasonic transmissions Distance and by river data measuring method comprising a procedure of outputting a moving direction, is solved.

上記河川データ測定方法において、超音波発信から河底の反射波受信までの時間と河底からの反射波の強度を測定し、超音波発信から河底の反射波受信までの時間を入力として水深を算出して記憶手段に格納し、さらに、前記移動ベクトルを順次結合して測定装置の移動の軌跡を生成するとともに、前記移動の軌跡を構成する各移動ベクトルの終点位置に対応する水深のデータを前記記憶手段から読み出して前記終点位置のデータに組み合わせて前記記憶手段に格納する手順を有してなる河川データ測定方法としてもよい。   In the river data measurement method described above, the time from the transmission of ultrasonic waves to the reception of reflected waves at the bottom of the river and the intensity of reflected waves from the bottom of the river are measured. Is calculated and stored in the storage means, and further, the movement vectors are sequentially combined to generate a movement trajectory of the measuring apparatus, and water depth data corresponding to the end point position of each movement vector constituting the movement trajectory. May be read out from the storage means, combined with the data of the end point position, and stored in the storage means.

本発明の前記第2の課題は、上記河川データ測定方法において、超音波発信から河底の反射波受信までの時間と河底からの反射波の強度を測定し、河底で散乱・反射された超音波に加え、水中に浮遊する水中反射体で反射された反射波を前記2次元超音波受信手段で受信するとともに、超音波発信から水中反射体で反射された反射波を検出するまでの時間と該反射波の強度分布及び位相分布をも測定し、それぞれ異なる深さでかつ超音波発信源との水平方向相対位置が所定の位置にある水中反射体によって前記超音波送受信手段から発信された超音波が反射されたときの反射波の強度分布及び位相分布を基準強度分布及び基準位相分布として前記記憶手段に格納され、前記超音波発信から河底の反射波受信までの時間を入力として水深を算出して前記記憶手段に格納し、前記超音波発信から水中反射体の反射波受信までの時間に基づいて水中反射体が浮遊している深さを検出し、前記検出された深さの水中反射体の反射波の基準強度分布及び基準位相分布と前記測定された反射波の強度分布及び位相分布のいずれか若しくは双方の相関の強度を、測定された分布を基準分布に近づけつつ移動距離の関数として求め、相関の強度が最大となる移動距離及びそのときの移動方向に基づいて当該水中反射体の超音波発信源に対する水平方向位置を求め、所定の時間間隔で検出された水中反射体の超音波発信源に対する水平方向相対位置及び前記算出された超音波発信源の移動ベクトルに基づいて、前記移動ベクトルの終点位置における、前記検出された深さの流速を求める手順を有してなる河川データ測定方法によっても解決される。   The second problem of the present invention is that, in the river data measurement method, the time from ultrasonic transmission to reception of reflected waves on the riverbed and the intensity of reflected waves from the riverbed are measured and scattered and reflected by the riverbed. In addition to the ultrasonic wave, the reflected wave reflected by the underwater reflector floating in the water is received by the two-dimensional ultrasonic wave receiving means, and the reflected wave reflected by the underwater reflector is detected from the ultrasonic wave transmission. The time and intensity distribution and phase distribution of the reflected wave are also measured and transmitted from the ultrasonic transmission / reception means by an underwater reflector having different depths and a horizontal relative position to the ultrasonic transmission source at a predetermined position. The intensity distribution and the phase distribution of the reflected wave when the reflected ultrasonic wave is reflected are stored in the storage means as the reference intensity distribution and the reference phase distribution, and the time from the ultrasonic transmission to the reception of the reflected wave on the riverbed is input. Calculate water depth And detecting the depth at which the underwater reflector is floating based on the time from the transmission of the ultrasonic wave to the reception of the reflected wave of the underwater reflector, and the underwater reflector of the detected depth The intensity of the correlation between the reference intensity distribution and reference phase distribution of the reflected wave and the intensity distribution and / or phase distribution of the measured reflected wave is used as a function of the moving distance while bringing the measured distribution close to the reference distribution. Obtain the horizontal position of the underwater reflector with respect to the ultrasonic wave transmission source based on the moving distance and the moving direction at which the correlation intensity is maximum, and detect the ultrasonic wave of the underwater reflector detected at a predetermined time interval. Based on the horizontal relative position with respect to the transmission source and the calculated movement vector of the ultrasonic transmission source, a procedure for obtaining the flow velocity of the detected depth at the end point position of the movement vector is provided. Also solved by river data measurement method.

本発明によれば、測定装置の移動量と移動方向を十分な精度で測定することが可能になり、河床形状の測定精度が向上するとともに、流速分布、流量の測定精度が向上する効果が得られる。   According to the present invention, it becomes possible to measure the moving amount and moving direction of the measuring apparatus with sufficient accuracy, and the measurement accuracy of the riverbed shape is improved, and the effect of improving the measurement accuracy of the flow velocity distribution and the flow rate is obtained. It is done.

本発明では、水面付近に配置した超音波測定装置で河床形状、流速分布、流量などの河川データを測定する。測定線に沿った河底の形状は、指向性を小さくした超音波送受信器からの超音波を河底に向けて照射し、照射開始から反射波検出までの時間とあらかじめ知られている水中の音速から河底までの距離を算出する。流速分布は、上記の河底に向けて放射した超音波ビームのビーム領域中にある水中反射体がある時間中に移動する量を、反射体の球面反射波から求め、流速に換算する。これには、上記の超音波送受信器とは別個に2次元状に配置した超音波検出器を用い、この検出器で反射体の球面反射波の分布を検出する。   In the present invention, river data such as river bed shape, flow velocity distribution, and flow rate is measured by an ultrasonic measurement device disposed near the water surface. The shape of the riverbed along the measurement line irradiates the riverbed with ultrasonic waves from an ultrasonic transmitter / receiver with reduced directivity, the time from the start of irradiation until detection of the reflected wave, and the known underwater Calculate the distance from the speed of sound to the riverbed. The flow velocity distribution is obtained by calculating the amount of movement of the underwater reflector in the beam region of the ultrasonic beam radiated toward the riverbed during a certain time from the spherical reflected wave of the reflector and converting it into the flow velocity. For this purpose, an ultrasonic detector arranged two-dimensionally separately from the ultrasonic transmitter / receiver is used, and the distribution of the spherical reflected wave of the reflector is detected by this detector.

河川水面に浮かべた測定装置の位置は、河底の散乱で形成される超音波スペックルパターンの移動による変化を検出して求める。   The position of the measuring device floated on the river surface is obtained by detecting changes due to the movement of the ultrasonic speckle pattern formed by scattering of the riverbed.

河床形状、流速分布、流量の計測に必要な測定装置の位置はこの移動計測結果から算出する。
<第1の実施の形態>
本発明の第1の実施の形態を、図面を参照して説明する。図1に、河床形状、流速分布、流量の測定の基本概念を示す。測定は、図1に示すように、河川1に、河川1を横切る方向に、ある測定線2を設定し、測定線2に沿った河底の形状(河床形状)や河川断面の流速分布、流量を求める。このために、測定装置を搭載した浮揚体3を、自走或いは牽引などの方法で、測定線2に沿って移動させる。
The position of the measuring device necessary for measuring the riverbed shape, velocity distribution, and flow rate is calculated from the result of this movement measurement.
<First Embodiment>
A first embodiment of the present invention will be described with reference to the drawings. FIG. 1 shows the basic concept of riverbed shape, flow velocity distribution, and flow rate measurement. As shown in FIG. 1, a measurement line 2 is set in a direction crossing the river 1 in the river 1, and the shape of the river bed (river bed shape) along the measurement line 2 and the flow velocity distribution of the river cross section, Find the flow rate. For this purpose, the levitation body 3 equipped with the measuring device is moved along the measurement line 2 by a method such as self-propelled or towed.

図2に、本実施の形態に係る河川データ測定装置の要部構成を示す。図示の河川データ測定装置は、浮揚体3に搭載されており、浮揚体3底部外面に配置されて河底方向に超音波を放射する超音波送受信器301と、超音波送受信器301に接続された送受信回路302と、浮揚体3底部外面に配置された2次元超音波検出器303と、2次元超音波検出器303に接続された信号処理装置304と、前記送受信回路302及び信号処理装置304に接続された信号解析装置305と、信号解析装置305に接続された図示されていない記憶手段であるメモリ、図示されていない画面表示手段、入力手段及び電源とを含んで構成されている。   In FIG. 2, the principal part structure of the river data measuring apparatus which concerns on this Embodiment is shown. The river data measuring apparatus shown in the figure is mounted on the levitation body 3, and is connected to the ultrasonic transceiver 301, which is disposed on the outer surface of the bottom of the levitation body 3 and emits ultrasonic waves toward the riverbed. Transmission / reception circuit 302, a two-dimensional ultrasonic detector 303 disposed on the bottom outer surface of the levitation body 3, a signal processing device 304 connected to the two-dimensional ultrasonic detector 303, the transmission / reception circuit 302 and the signal processing device 304. The signal analysis device 305 is connected to the signal analysis device 305, the memory (not shown) is a storage unit (not shown), the screen display unit (not shown), the input unit, and the power source.

浮揚体3は、図2に示すように、測定装置を搭載して水面に浮いており、自走などにより水面を移動する。図2では、浮揚体3を水面上で移動させるための自走装置などは省略してあり、本発明に係わる測定装置部分のみを示してある。   As shown in FIG. 2, the levitation body 3 is equipped with a measuring device and floats on the water surface, and moves on the water surface by self-propelled or the like. In FIG. 2, a self-propelled device for moving the levitation body 3 on the surface of the water is omitted, and only the measurement device portion according to the present invention is shown.

河底反射測定手段である送受信回路302は、超音波送受信手段である超音波送受信器301から河底に向かって超音波を発信させ、反射波を超音波送受信器301を介して受信する。超音波送受信器301から送信された超音波は水中を伝搬し、河底面で反射してその反射波が超音波送受信器301で検出される。   The transmission / reception circuit 302 that is a riverbed reflection measurement unit transmits an ultrasonic wave toward the riverbed from an ultrasonic transmitter / receiver 301 that is an ultrasonic transmission / reception unit, and receives the reflected wave via the ultrasonic transmitter / receiver 301. The ultrasonic wave transmitted from the ultrasonic transmitter / receiver 301 propagates in water, is reflected by the riverbed, and the reflected wave is detected by the ultrasonic transmitter / receiver 301.

超音波送受信器301から河底に向かって発信された超音波が水中に浮遊している物体(水中反射体)に当たると、超音波は反射されて球面波となって浮揚体3の方向にも伝搬する。また、河底で反射・散乱した超音波も浮揚体3の方向に伝搬する。なお、本実施の形態では、超音波送受信器301は、超音波の送受信を行なうが、水深計測を別の装置で行う場合は、超音波送信手段として、送信専用としてもよい。   When the ultrasonic wave transmitted from the ultrasonic transmitter / receiver 301 toward the river bottom strikes an object floating in water (an underwater reflector), the ultrasonic wave is reflected and becomes a spherical wave in the direction of the levitation body 3. Propagate. In addition, ultrasonic waves reflected and scattered by the riverbed also propagate in the direction of the levitation body 3. In this embodiment, the ultrasonic transmitter / receiver 301 transmits / receives ultrasonic waves. However, when the water depth measurement is performed by another apparatus, the ultrasonic transmitter / receiver 301 may be dedicated to transmission.

2次元超音波受信手段である2次元超音波検出器303は2次元平面に密に配置された複数の検出器で構成され、個々の検出器はこれらの球面波、散乱波の強度、位相をそれぞれ検出し、信号処理手段である信号処理装置304に伝送する。   A two-dimensional ultrasonic detector 303, which is a two-dimensional ultrasonic receiver, is composed of a plurality of detectors arranged densely on a two-dimensional plane. Each detector detects the intensity and phase of these spherical waves, scattered waves. Each is detected and transmitted to a signal processing device 304 as signal processing means.

信号処理装置304は、2次元超音波検出器303で得た信号を処理する装置であり、入力された信号の増幅やフィルタリングを行ったあと強度、位相を抽出し、アナログ−デジタル変換して、その結果と検出タイミング信号を信号解析手段である信号解析装置305に転送する。   The signal processing device 304 is a device that processes the signal obtained by the two-dimensional ultrasonic detector 303, and after performing amplification and filtering of the input signal, extracts the intensity and phase, and performs analog-digital conversion, The result and the detection timing signal are transferred to the signal analysis device 305 which is signal analysis means.

信号解析装置305はメモリを備えた演算手段、いわゆるマイクロコンピュータであり、図示されていないメモリが接続されており、実行する処理内容はソフトウエアとして搭載されている。信号解析装置305は、図6に示す送信超音波波形信号W02を生成して送受信回路302に送るとともに、位相検出のための基準波形信号W01を生成して信号処理装置304に送る。また、前記メモリに、超音波送受信器301の真下で水面から予め定められた複数の深さレベルdにある反射体で反射された球面波、散乱波が前記2次元超音波検出器303で受信され、信号処理装置304で処理された時の強度、位相の分布状態のデータが予め算出されて基準強度分布、基準位相分布のパターンとして格納されている。前記複数の深さレベルdは、d=0を水面として、d=1,2,…i,…,mまで、例えば0.5m間隔で設定されている。前記ソフトウエアには、流速測定の際、演算対象となる深さレベルを指定するパラメータd=iが含まれており、初期値は、i=1に設定されている。本実施の形態では、1MHzの超音波を用いる。   The signal analysis device 305 is a calculation means including a memory, a so-called microcomputer, and is connected to a memory not shown, and the processing content to be executed is installed as software. The signal analysis device 305 generates a transmission ultrasonic waveform signal W02 shown in FIG. 6 and sends it to the transmission / reception circuit 302, and also generates a reference waveform signal W01 for phase detection and sends it to the signal processing device 304. In addition, the two-dimensional ultrasonic detector 303 receives spherical waves and scattered waves reflected by a reflector at a plurality of predetermined depth levels d from the water surface directly below the ultrasonic transmitter / receiver 301 in the memory. The intensity and phase distribution data when processed by the signal processing device 304 are calculated in advance and stored as a reference intensity distribution and reference phase distribution pattern. The depth levels d are set to d = 1, 2,..., I,..., M, for example, at intervals of 0.5 m, where d = 0 is the water surface. The software includes a parameter d = i for designating a depth level to be calculated in the flow velocity measurement, and the initial value is set to i = 1. In this embodiment, 1 MHz ultrasonic waves are used.

なお、本実施の形態の河川データ測定装置には、装置全体の電源を投入・遮断する電源スイッチのほかに、超音波の生成・発信を停止する超音波発信停止スイッチが設けられ、測定終了時、超音波の生成・発信を停止した状態(超音波発信停止スイッチオフ)で、データ処理を行なうことが出来るようになっている。   The river data measurement device of this embodiment is provided with an ultrasonic transmission stop switch for stopping the generation and transmission of ultrasonic waves, in addition to the power switch for turning on / off the power of the entire device. Data processing can be performed in a state where generation / transmission of ultrasonic waves is stopped (ultrasonic transmission stop switch is off).

以下、流速の検出について図3を参照して説明する。流速は、水中に浮遊して移動している物体(散乱・反射体)が反射する超音波を、時間をおいて複数回受信することで求められる。通常、河川の流水中には、超音波を反射するに足る大きさの微小な物体が種々の深さで浮遊し、流水とともに移動している。河川の流速は、それら物体が、超音波送受信器301から発信された超音波ビーム中にあって移動するとき、それら物体が反射する超音波(反射波)を、その移動の前後で検出して求められる。   Hereinafter, detection of the flow velocity will be described with reference to FIG. The flow velocity is obtained by receiving ultrasonic waves reflected by an object (scattering / reflecting body) floating and moving in water a plurality of times over time. Usually, in a flowing water of a river, a minute object large enough to reflect ultrasonic waves floats at various depths and moves together with the flowing water. As for the flow velocity of the river, when these objects move in the ultrasonic beam transmitted from the ultrasonic transceiver 301, the ultrasonic waves (reflected waves) reflected by these objects are detected before and after the movement. Desired.

図3に示すように、超音波ビーム中の位置Aに反射体があると、その反射波は球面波となって2次元超音波検出器303で検出される。球面波は2次元分布であるが、簡単のため図3では一方向のみの分布を示している。たとえば、2回の超音波送信で、1回目にある深さで位置Aに反射体が検出され、2回目に同じ深さで位置Aには反射体がなく、位置Bで反射体が検出されたら、反射体が、2回の超音波送信の間に位置Aから位置Bに移動したと考える。なお、2次元超音波検出器303は、超音波ビーム中の種々の深さにある反射体すべての反射波を受信するが、信号解析装置305は、前記予め定められた深さにある反射体からの反射波のデータを選別して解析する。つまり、超音波の発信後、前記定められた深さレベルそれぞれに対応する時間後に受信された反射波のデータを選別して解析する。   As shown in FIG. 3, if there is a reflector at position A in the ultrasonic beam, the reflected wave becomes a spherical wave and is detected by the two-dimensional ultrasonic detector 303. The spherical wave has a two-dimensional distribution, but for simplicity, FIG. 3 shows a distribution in only one direction. For example, in two ultrasonic transmissions, a reflector is detected at position A at the first depth, and at the second depth, there is no reflector at position A, and a reflector is detected at position B. Then, it is considered that the reflector has moved from position A to position B between the two ultrasonic transmissions. The two-dimensional ultrasonic detector 303 receives the reflected waves of all the reflectors at various depths in the ultrasonic beam, but the signal analysis device 305 receives the reflector at the predetermined depth. Select and analyze the reflected wave data from That is, after transmission of ultrasonic waves, reflected wave data received after a time corresponding to each of the predetermined depth levels is selected and analyzed.

反射・散乱波は、位置A、位置Bを中心とする球面波となって2次元超音波検出器303で検出される。図3に示すように、強度分布はほぼ位置A、位置Bの真上で大きく、真上から離れるに従って小さくなる特性となる。超音波送信から検出深さに応じた時間後の球面波位相は、位置A、位置Bの真上でほぼ平坦で、離れるのに従って遅れる位相となる。図3の位相は、±180度の範囲を示しているため、途中で切れた特性となっている。   The reflected / scattered wave is a spherical wave centered at position A and position B and is detected by the two-dimensional ultrasonic detector 303. As shown in FIG. 3, the intensity distribution has a characteristic that is almost right above position A and position B, and becomes smaller as it goes away from right above. The spherical wave phase after a time corresponding to the detection depth from the ultrasonic transmission is substantially flat immediately above the position A and the position B, and becomes a phase that is delayed as the distance increases. Since the phase in FIG. 3 indicates a range of ± 180 degrees, the phase is cut off halfway.

このように、強度、位相とも反射・散乱体の位置を仮定すれば分布特性を推測でき、あらかじめ計算し前記メモリに格納されている基準強度分布、基準位相分布のパターンと比較することによって逆に反射・散乱体が超音波送受信器301の真下の位置からどの方向にどれくらい離れた位置にあるかを算出できる。位置A、Bが求まると、超音波の送信間隔からその移動速度、つまり、その深さでの流速が算出できる。算出されたデータは、前記2回の超音波発信のうちの後の超音波発信の時刻データ、深さレベルデータとともにメモリに格納される。   In this way, the distribution characteristics can be estimated if the position of the reflector / scatterer is assumed for both the intensity and phase, and by comparing with the reference intensity distribution and reference phase distribution patterns stored in the memory in advance. It is possible to calculate how far and in what direction the reflector / scatterer is from the position directly below the ultrasonic transceiver 301. When the positions A and B are obtained, the moving speed, that is, the flow velocity at the depth can be calculated from the transmission interval of the ultrasonic waves. The calculated data is stored in the memory together with the time data and depth level data of the subsequent ultrasonic transmission of the two ultrasonic transmissions.

図3では一方向の反射体の移動について述べたが、2次元超音波検出器303は2次元面のデータを採取するため、強度、位相の2次元状の分布パターンをあらかじめ基準強度分布パターン、基準位相分布パターンとして求めて格納しておくことで、二次元平面での流速分布を検出できる。さらに、位相の算出時、超音波送信から検出深さに応じた時間を同じとすることで、A、B真上での位相の差から超音波の放射から検出までの時間を求めることができ、深さ方向の位置、つまり、深さ方向の流速も算出可能である。これにより、ある深さの3次元流速分布を検出できる。   Although the movement of the reflector in one direction is described in FIG. 3, the two-dimensional ultrasonic detector 303 collects two-dimensional surface data, so that a two-dimensional distribution pattern of intensity and phase is set in advance as a reference intensity distribution pattern, By obtaining and storing as a reference phase distribution pattern, a flow velocity distribution in a two-dimensional plane can be detected. Furthermore, when calculating the phase, the time from ultrasonic transmission to the detection depth can be made the same, so that the time from ultrasonic radiation to detection can be obtained from the phase difference directly above A and B. The position in the depth direction, that is, the flow velocity in the depth direction can also be calculated. Thereby, a three-dimensional flow velocity distribution of a certain depth can be detected.

次に、超音波のスペックルパターンを用いた浮揚体3、つまり測定装置の移動量と移動方向を求める方法について図4を参照して説明する。超音波送受信器301から河底に向かって送信された超音波は、河底に到達し、反射する。河底が鏡面であればその反射は入射角と反射角の関係で反射特性が決まる。しかし、通常、河底表面は砂や小石で覆われている。超音波の波長より砂や小石の凹凸が大きいと、超音波は鏡面反射とならず、凹凸により種々の方向に散乱される状態となる。水中を伝播する1MHzの超音波の波長は約1.5mmであり、河底表面の砂や小石の凹凸はそれよりも大きいから、反射波は種々の方向に散乱される。   Next, a method for obtaining the levitation body 3 using the ultrasonic speckle pattern, that is, the moving amount and moving direction of the measuring apparatus will be described with reference to FIG. The ultrasonic waves transmitted from the ultrasonic transceiver 301 toward the river bed reach the river bed and are reflected. If the riverbed is a mirror surface, its reflection characteristics are determined by the relationship between the incident angle and the reflection angle. However, the riverbed surface is usually covered with sand and pebbles. If the unevenness of sand or pebbles is larger than the wavelength of the ultrasonic wave, the ultrasonic wave is not mirror-reflected and is scattered in various directions by the unevenness. Since the wavelength of 1 MHz ultrasonic waves propagating in water is about 1.5 mm, and the unevenness of sand and pebbles on the riverbed surface is larger than that, reflected waves are scattered in various directions.

照射面の各点から散乱された超音波は水中空間で互いに干渉し合い、超音波の重ね合わせで強度の大きい部分と小さい部分の強度むら(スペックルパターン)を形成する。この強度むらを、浮揚体3の底面に設置した2次元超音波検出器303を介して信号処理装置304で検出する。2次元超音波検出器303では、浮揚体3の底面である受信面で、スペックルパターンである強度分布を受信することになる。   The ultrasonic waves scattered from each point on the irradiated surface interfere with each other in the underwater space, and form an intensity unevenness (speckle pattern) between the high intensity portion and the low intensity portion by superimposing the ultrasonic waves. This unevenness in intensity is detected by the signal processing device 304 via a two-dimensional ultrasonic detector 303 installed on the bottom surface of the levitation body 3. The two-dimensional ultrasonic detector 303 receives an intensity distribution that is a speckle pattern on the receiving surface that is the bottom surface of the levitation body 3.

超音波送受信器301は、ビームの指向角α度を持った円錐形に広がる超音波ビームを放射すると仮定する。2次元超音波検出器303を介して信号処理装置304で検出される強度むらには、種々の空間周波数が含まれるが、その強度むらの平均変化幅は、指向角α度と水中における超音波の波長できまり、深さに依存しない。これは、光の場合も同様であり、詳細は非特許文献1に述べられている。具体的な数値例をあげると、α=2.5度の超音波送信器を使い、周波数を1MHzにすると、水中の音速は1500m/s であるから、強度むらの平均変化幅は3.4cmとなる。2次元超音波検出器303の各検出器はこの平均変化幅より小さくする必要がある。また、α=2.5度 の超音波送信器で送信された超音波は、深さ5mで、直径22cmの円形の照射領域を持つ。   It is assumed that the ultrasonic transceiver 301 emits an ultrasonic beam spreading in a conical shape having a beam directivity angle α degrees. The intensity unevenness detected by the signal processing device 304 via the two-dimensional ultrasonic detector 303 includes various spatial frequencies. The average variation width of the intensity unevenness is determined by the directivity angle α degree and the ultrasonic wave in water. It depends on the wavelength, and does not depend on the depth. This is the same in the case of light, and details are described in Non-Patent Document 1. As a specific numerical example, if an ultrasonic transmitter with α = 2.5 degrees is used and the frequency is set to 1 MHz, the speed of sound in water is 1500 m / s. Therefore, the average variation width of the intensity unevenness is 3.4 cm. It becomes. Each detector of the two-dimensional ultrasonic detector 303 needs to be smaller than the average change width. Further, the ultrasonic wave transmitted by the ultrasonic transmitter with α = 2.5 degrees has a circular irradiation area with a depth of 5 m and a diameter of 22 cm.

図4に強度分布の例を示す。超音波送受信器301から発信された超音波ビームによる河底の散乱超音波は、2次元超音波検出器303で受信されて2次元の強度分布を示すが、図4では簡単のため一方向の散乱強度データを例として示した。浮揚体が実線で示すAFの位置にあると、スッペックルパターンの強度分布は、図中の<AFでの強度分布>となり、これが検出される。浮揚体が位置AFから破線で示す位置BFに移動した時を考える。位置AFから位置BFへの移動量が河底の超音波照射領域の大きさより十分小さければ、<BFでの強度分布>は、位置BFでの照射部分が位置AFでの照射部分と少し異なるため、<AFでの強度分布>と分布の形は若干違う。しかし、照射領域の大部分が位置AF、位置BFで同じであるため、分布の差は小さく、<BFでの強度分布>は、位置AFと検出位置が違うため、<AFでの強度分布>が位置AFから位置BFへの移動量だけシフトした強度分布波形に近い。なお、図5に示す<AFでの強度分布>と<BFでの強度分布>の横軸は、2次元超音波検出器303を構成する個々の検出器の位置に対応している。   FIG. 4 shows an example of the intensity distribution. Scattering ultrasonic waves in the riverbed by the ultrasonic beam transmitted from the ultrasonic transceiver 301 are received by the two-dimensional ultrasonic detector 303 and show a two-dimensional intensity distribution. In FIG. Scattering intensity data is shown as an example. When the levitation body is at the AF position indicated by the solid line, the intensity distribution of the speckle pattern becomes <intensity distribution at AF> in the figure, and this is detected. Consider a case where the levitation body has moved from position AF to position BF indicated by a broken line. If the amount of movement from the position AF to the position BF is sufficiently smaller than the size of the ultrasonic irradiation area on the riverbed, the <intensity distribution at BF> is slightly different from the irradiation part at the position AF. <Intensity distribution in AF> and the shape of the distribution are slightly different. However, since most of the irradiation area is the same at position AF and position BF, the difference in distribution is small, and <intensity distribution at BF> is different from the position AF and detection position, so <intensity distribution at AF> Is close to the intensity distribution waveform shifted by the amount of movement from the position AF to the position BF. Note that the horizontal axis of <intensity distribution at AF> and <intensity distribution at BF> shown in FIG. 5 corresponds to the position of each detector constituting the two-dimensional ultrasonic detector 303.

このシフト量の検出のため、図5に示すように、<AFでの強度分布>パターンと<BFでの強度分布>パターンの相関演算を行う。図5では、簡単のため一方向の相関計算を示してある。<BFでの強度分布>パターンを基準とし、<BFでの強度分布>パターンに対して<AFでの強度分布>パターンを予め定めた距離だけ移動させつつ二つの強度分布の相関演算を行う。図5の右側の図は、移動量を横軸に、相関の強度を縦軸にとって、相関演算の結果を示している。相関演算を行うと二つの強度分布が最も一致した移動量の位置で相関の強度が最大になる。このピーク位置を求めることで浮揚体3の位置AFと位置BFの距離差が検出される。   In order to detect this shift amount, as shown in FIG. 5, a correlation calculation is performed between the <intensity distribution at AF> pattern and the <intensity distribution at BF> pattern. In FIG. 5, one-way correlation calculation is shown for simplicity. The <intensity distribution at BF> pattern is used as a reference, and the correlation between the two intensity distributions is performed while the <intensity distribution at AF> pattern is moved by a predetermined distance with respect to the <intensity distribution at BF> pattern. The diagram on the right side of FIG. 5 shows the result of the correlation calculation with the movement amount on the horizontal axis and the correlation strength on the vertical axis. When the correlation calculation is performed, the correlation intensity becomes the maximum at the position of the movement amount where the two intensity distributions are most consistent. By obtaining this peak position, the distance difference between the position AF and the position BF of the levitation body 3 is detected.

図5は一方向の相関演算であるが、2次元超音波検出器303を介して信号処理装置304で2次元の強度分布を求め、2次元相関演算を行うと、浮揚体3がどちらの方向に、どれだけの距離を移動したか、つまり、方向と大きさを持ったベクトルとして浮揚体の移動を算出できる。このベクトルを順次連結することによりある基準点からの浮揚体3の水平方向移動軌跡を求めることが可能となる。得られた水平方向位置データは、位置BFにおいて超音波を発信したときの時刻データとともに、メモリに格納される。   FIG. 5 shows a correlation calculation in one direction. When a two-dimensional intensity calculation is performed by the signal processing device 304 via the two-dimensional ultrasonic detector 303 and the two-dimensional correlation calculation is performed, the levitation body 3 moves in which direction. In addition, it is possible to calculate how much the levitated body has moved, that is, as a vector having a direction and a size. By sequentially connecting these vectors, it is possible to obtain the horizontal movement locus of the levitation body 3 from a certain reference point. The obtained horizontal position data is stored in the memory together with time data when an ultrasonic wave is transmitted at the position BF.

ここで検出波形の例を示して、上記計測法を述べる。図6に、流速を求める時の波形の関係を示す。周波数fなる正弦波を基準波形W01とする。その一部の波形(送信超音波波形W02)を超音波送受信器301から放射する。2次元超音波検出器303のある任意の検出器で、超音波ビーム内で反射体位置AとBとで検出した波形(W03、W06)を得る。波形W03、W06は、送信超音波波形W02をある間隔で放射した時の、2次元超音波検出器303の同一検出器の検出波形であり、この間に反射体は点Aから点Bに移動したと仮定する。検出波形で強度W04、W07と位相W05、W08を算出する。強度は、例えば包絡線検波で求まる。位相は、基準波形と検出波形の積を低域ろ波することで算出できる。tを時間、ΔTを超音波放射から検出までの時間、検出信号の振幅をVとし、基準波形をsin(2πft)、検出波形をV*sin(2πf(t+ΔT))とすれば、その積の低域ろ波信号は、0.5V*cos(2πf*ΔT)となる。また、cos(2πft)なる波形と検出波形の積、低域ろ波処理で0.5V*sin(2πf*ΔT)が求まり、両者からtan(2πf*ΔT)を得て、位相を算出できる。   Here, the measurement method will be described with an example of a detected waveform. FIG. 6 shows the relationship of waveforms when obtaining the flow velocity. A sine wave having a frequency f is set as a reference waveform W01. A part of the waveform (transmission ultrasonic waveform W02) is radiated from the ultrasonic transceiver 301. Waveforms (W03, W06) detected at the reflector positions A and B in the ultrasonic beam by an arbitrary detector having the two-dimensional ultrasonic detector 303 are obtained. Waveforms W03 and W06 are detection waveforms of the same detector of the two-dimensional ultrasonic detector 303 when the transmission ultrasonic waveform W02 is emitted at a certain interval. During this time, the reflector has moved from point A to point B. Assume that Intensities W04 and W07 and phases W05 and W08 are calculated from the detected waveform. The intensity is obtained, for example, by envelope detection. The phase can be calculated by low-pass filtering the product of the reference waveform and the detected waveform. If t is time, ΔT is the time from ultrasonic radiation to detection, the amplitude of the detection signal is V, the reference waveform is sin (2πft), and the detection waveform is V * sin (2πf (t + ΔT)), the product of The low-pass filtered signal is 0.5 V * cos (2πf * ΔT). Further, 0.5 V * sin (2πf * ΔT) is obtained by the product of the waveform of cos (2πft) and the detected waveform, and low-pass filtering, and tan (2πf * ΔT) is obtained from both to obtain the phase.

位相の波形W05、W08で、意味をもつのは反射体の反射強度がある部分のみで、そのほかの部分は不定になる可能性がある。このような波形の処理手順で、反射体の移動前後の強度、位相が算出できる。この波形は、検出器1個の信号について示したが、複数の検出器の信号について同じ処理を行うことにより、反射体の移動前後で強度分布、位相分布のパターンを得ることができる。このパターンと、あらかじめ求めておいたパターンとの一致性から反射体が移動した距離を算定する。さらに、2回の超音波放射の時間間隔を考慮して、反射体の移動速度、つまり、流速を検出する。   In the phase waveforms W05 and W08, only the portion having the reflection intensity of the reflector is meaningful, and the other portions may be indefinite. With such a waveform processing procedure, the intensity and phase before and after the movement of the reflector can be calculated. This waveform is shown for the signal of one detector, but by performing the same processing on the signals of a plurality of detectors, patterns of intensity distribution and phase distribution can be obtained before and after the movement of the reflector. The distance traveled by the reflector is calculated from the coincidence between this pattern and the previously obtained pattern. Further, the moving speed of the reflector, that is, the flow velocity is detected in consideration of the time interval between the two ultrasonic radiations.

図7は、水深と浮揚体3の移動量を検出する際の波形例である。この場合、超音波の放射は1回(送信超音波波形W12)で水深とスペックルパターンが求まる。超音波送受信器301の受信波形(W13)で、河底の反射位置を検出し、これより水深がわかる。前記位置BFで発信された超音波で測定された水深データは、発信超音波の時刻データとともにメモリに格納される。2次元超音波検出器303で異なる2個の検出器で検出した波形(W14、W16)を検波して強度を求める。河底反射位置で反射、散乱された超音波を2次元超音波検出器303で受信した時の反射波強度分布がスペックルパターンとなる。この強度分布が、浮揚体3の移動で変化する。浮揚体3の移動前後でこのスペックルパターンを求め、両者の相関のピークから移動ベクトルを計算する。   FIG. 7 is a waveform example when detecting the water depth and the amount of movement of the levitation body 3. In this case, the depth of water and the speckle pattern can be obtained by emitting ultrasonic waves once (transmitted ultrasonic waveform W12). The reflection position of the riverbed is detected by the received waveform (W13) of the ultrasonic transmitter / receiver 301, and the water depth can be known from this. The water depth data measured by the ultrasonic wave transmitted at the position BF is stored in the memory together with the time data of the transmitted ultrasonic wave. The two-dimensional ultrasonic detector 303 detects waveforms (W14, W16) detected by two different detectors, and obtains the intensity. The reflected wave intensity distribution when the two-dimensional ultrasonic detector 303 receives the ultrasonic wave reflected and scattered at the riverbed reflection position becomes a speckle pattern. This intensity distribution changes as the levitation body 3 moves. This speckle pattern is obtained before and after the levitation body 3 moves, and a movement vector is calculated from the peak of the correlation between the two.

上記信号処理は主に信号処理装置304で実行される。信号処理装置304の詳細を図8によって説明する。2次元超音波検出器303は2次元平面状に配置された複数の検出器で構成され、前記検出器それぞれに個別に信号処理装置304の処理回路304iが接続されている。   The signal processing is mainly performed by the signal processing device 304. Details of the signal processing device 304 will be described with reference to FIG. The two-dimensional ultrasonic detector 303 includes a plurality of detectors arranged in a two-dimensional plane, and a processing circuit 304i of the signal processing device 304 is individually connected to each of the detectors.

各処理回路304iは、前記検出器に接続されたフィルターアンプ304i−1と、フィルターアンプ304i−1に接続された検波器304i−2及び位相検出器304i−3と、前記検波器304i−2に接続されたレベル判定器304i−4と、前記検波器304i−2及びレベル判定器304i−4に接続されたスイッチ素子304i−5と、スイッチ素子304i−5に接続されたアナログ−デジタル変換器304i−7と、前記レベル判定器304i−4及び位相検出器304i−3に接続されたスイッチ素子304i−6と、スイッチ素子304i−6に接続されたアナログ−デジタル変換器304i−8とを含んで構成されている。   Each processing circuit 304i includes a filter amplifier 304i-1 connected to the detector, a detector 304i-2 and a phase detector 304i-3 connected to the filter amplifier 304i-1, and a detector 304i-2. The level determiner 304i-4 connected, the switch element 304i-5 connected to the detector 304i-2 and the level determiner 304i-4, and the analog-digital converter 304i connected to the switch element 304i-5. -7, a switch element 304i-6 connected to the level determiner 304i-4 and the phase detector 304i-3, and an analog-digital converter 304i-8 connected to the switch element 304i-6. It is configured.

フィルターアンプ304i−1は、検出信号のノイズを低減し、信号レベルを増大させる。検波器304i−2は、検出信号の振幅波形を得る。位相検出器304i−3は、検出信号波形と基準信号とから上記の演算手順で位相を算出する。レベル判定器304i−4は、検波器304i−2の出力があらかじめ設定した振幅以上のとき信号を出し、スイッチ素子304i−5、6を導通させる。レベル判定器304i−4は、超音波を送信したタイミング信号を受けており、これを基準に超音波反射波を受信する時間範囲を規定する。アナログ−デジタル変換機304i−7、8は、スイッチ素子304i−5、6を通って来た信号、つまり、振幅と位相情報をデジタル信号に変換して出力する。   The filter amplifier 304i-1 reduces the noise of the detection signal and increases the signal level. The detector 304i-2 obtains the amplitude waveform of the detection signal. The phase detector 304i-3 calculates the phase by the above calculation procedure from the detection signal waveform and the reference signal. The level determiner 304i-4 outputs a signal when the output of the detector 304i-2 is equal to or larger than a preset amplitude, and makes the switch elements 304i-5 and 6 conductive. The level determiner 304i-4 receives a timing signal at which an ultrasonic wave is transmitted, and defines a time range for receiving an ultrasonic reflected wave based on this timing signal. The analog-to-digital converters 304i-7 and 8 convert the signals that have passed through the switch elements 304i-5 and 6, that is, amplitude and phase information, into digital signals and output them.

処理回路304iを構成する要素は、いずれもICにより製作可能であり、小型にできる。各処理回路304iでこれらの強度、位相信号を得て、信号解析装置305に送信し、信号解析装置305が強度、位相の分布情報としてまとめる。   Any of the elements constituting the processing circuit 304i can be manufactured by an IC and can be made compact. Each processing circuit 304i obtains these intensity and phase signals and transmits them to the signal analysis device 305. The signal analysis device 305 collects the information as intensity and phase distribution information.

送受信回路302は、図9に示すように、超音波送受信器301に接続されたスイッチ素子302−1と、スイッチ素子302−1に接続されたフィルターアンプ302−3及び増幅器302−2と、フィルターアンプ302−3に接続された検波器302-4と、検波器302-4に接続されたレベル判定器302−5と、を含んで構成されている。   As shown in FIG. 9, the transmission / reception circuit 302 includes a switch element 302-1 connected to the ultrasonic transmitter / receiver 301, a filter amplifier 302-3 and an amplifier 302-2 connected to the switch element 302-1, and a filter. A detector 302-4 connected to the amplifier 302-3 and a level determiner 302-5 connected to the detector 302-4 are included.

増幅器302−2は、信号解析装置305から入力された送信超音波波形W02を増幅してスイッチ素子302−1に伝送する。スイッチ素子302−1は、超音波送受信器301の送信、受信を切り替えるもので、超音波送信時は、増幅器302−2で増幅された送信超音波波形W02が超音波送受信器301に伝送されるように、また、受信時は、超音波送受信器301とフィルターアンプ302−3がつながるように切替動作を行う。フィルターアンプ302−3は、スイッチ素子302−1を介して超音波送受信器301から入力された検出信号のノイズを低減し、信号レベルを増大させる。検波器302−4は、検出信号の振幅波形を得る。レベル判定器302−5は、検波器302-4の出力があらかじめ設定したレベルを超えたとき、河底の反射と判断してその検出タイミング信号を信号解析装置305に送る。   The amplifier 302-2 amplifies the transmission ultrasonic waveform W02 input from the signal analysis device 305 and transmits the amplified waveform to the switch element 302-1. The switch element 302-1 switches between transmission and reception of the ultrasonic transceiver 301. At the time of ultrasonic transmission, the transmission ultrasonic waveform W02 amplified by the amplifier 302-2 is transmitted to the ultrasonic transceiver 301. In addition, at the time of reception, the switching operation is performed so that the ultrasonic transceiver 301 and the filter amplifier 302-3 are connected. The filter amplifier 302-3 reduces noise in the detection signal input from the ultrasonic transceiver 301 via the switch element 302-1, and increases the signal level. The detector 302-4 obtains the amplitude waveform of the detection signal. When the output of the detector 302-4 exceeds a preset level, the level determiner 302-5 determines that the reflection is from the riverbed and sends the detection timing signal to the signal analyzer 305.

信号解析装置305の処理フローを図10によって説明する。計測処理の実施前に、流速測定に必要な、種々の深さの反射体から反射された反射球面波の強度パターン、位相パターンのデータをあらかじめ計算しておき、基準強度分布、基準位相分布のデータベースとして信号解析装置305のメモリに格納しておく。   The processing flow of the signal analyzer 305 will be described with reference to FIG. Before performing the measurement process, the intensity pattern and phase pattern data of the reflected spherical waves reflected from the reflectors of various depths necessary for the flow velocity measurement are calculated in advance, and the reference intensity distribution and the reference phase distribution are calculated. The data is stored in the memory of the signal analyzer 305 as a database.

計測を開始すると、まず、パラメータ設定が行なわれる(305−01)。これは、観測する深さ(最大深さ何メートルまで測定するか、つまりd=1〜n)や浮揚体3の移動情報(浮揚体3の移動速度を設定)などの計測条件の設定である。浮揚体3の移動速度を設定することにより、浮揚体3の移動量を測定する場合のデータ採取間隔が決まる。また、流速測定のためのデータ採取間隔を設定する。ここでいうデータ採取間隔は、移動量測定あるいは流速測定の演算を行なうとき、異なる時点で発信された2度の超音波のデータを用いるが、定められた時間間隔で発信される超音波に対し、何回おきのデータを用いるかということである。また、その時点における浮揚体3の位置を特定するデータ、つまり測定始点が記録される。例えば、GPSシステムにより決定される位置、あるいは、川岸に設定された基線との相対関係で決定される位置である。   When measurement is started, parameter setting is first performed (305-01). This is the setting of measurement conditions such as the depth to be observed (how many meters the maximum depth is measured, that is, d = 1 to n) and the movement information of the levitation body 3 (the movement speed of the levitation body 3 is set). . By setting the moving speed of the levitation body 3, the data collection interval when measuring the movement amount of the levitation body 3 is determined. Moreover, the data collection interval for the flow velocity measurement is set. The data collection interval here uses the data of two ultrasonic waves transmitted at different times when calculating the movement amount measurement or the flow velocity measurement, but for the ultrasonic waves transmitted at a predetermined time interval. How many times to use the data. Further, data for specifying the position of the levitation body 3 at that time, that is, the measurement start point is recorded. For example, the position determined by the GPS system or the position determined by the relative relationship with the base line set on the riverbank.

次に超音波信号を送信するタイミング信号、基準波形W01を信号処理装置304に送信するとともに、送信超音波波形W02を送受信回路302に送信する(305−02)。この段階で超音波送受信器301から超音波が送信される。ここでいうタイミング信号は超音波を発信する時間間隔を指示する信号であるが、この信号を受ける都度、超音波送受信器301から超音波が送信されるものであってもよい。   Next, a timing signal for transmitting an ultrasonic signal and a reference waveform W01 are transmitted to the signal processing device 304, and a transmission ultrasonic waveform W02 is transmitted to the transmission / reception circuit 302 (305-02). At this stage, ultrasonic waves are transmitted from the ultrasonic transceiver 301. The timing signal here is a signal that indicates a time interval for transmitting an ultrasonic wave, but an ultrasonic wave may be transmitted from the ultrasonic transceiver 301 each time this signal is received.

次いで、超音波送受信器301で受信した超音波信号を、送受信回路302で処理した結果(河底検出タイミング)と、2次元超音波検出器303で受信した超音波信号を、信号処理装置304で処理した結果(種々の深さの反射体からの反射強度パターン、位相パターン及びタイミング信号)を受信し、発信時刻データとともに格納する(305−03)。このとき、反射体からの反射強度パターン、位相パターンについて、同時に入力されるタイミング信号で反射体の深さレベルを判定し、深さレベルdが前記d=1〜nのいずれかに合致するものを、合致する深さレベルを示すd=1〜nのいずれかとともに格納する。   Next, the result of processing the ultrasonic signal received by the ultrasonic transmitter / receiver 301 by the transmitter / receiver circuit 302 (river bottom detection timing) and the ultrasonic signal received by the two-dimensional ultrasonic detector 303 are processed by the signal processing device 304. The processed results (reflection intensity patterns, phase patterns and timing signals from reflectors of various depths) are received and stored together with the transmission time data (305-03). At this time, with respect to the reflection intensity pattern and the phase pattern from the reflector, the depth level of the reflector is determined by a timing signal input at the same time, and the depth level d matches any of d = 1 to n. Are stored with any of d = 1 to n indicating the matching depth level.

さらに、受信した河底検出タイミングに基づいて河底深さ情報を算出して超音波発信時刻データとともに格納し(305−04)、河底からの反射超音波のスペックルパターンである強度分布を、同様に、超音波発信時刻データとともに格納する(305−05)。   Further, riverbed depth information is calculated based on the received riverbed detection timing and stored together with ultrasonic transmission time data (305-04), and an intensity distribution that is a speckle pattern of reflected ultrasonic waves from the riverbed is calculated. Similarly, it is stored together with the ultrasonic transmission time data (305-05).

次いで手順305−06に進み、今回得られた強度分布と予め設定された前記データ採取間隔だけ以前に格納された強度分布の相関演算を行ない、相関の強度が最大となる位置を求めて、浮揚体3の移動ベクトルを算出し、移動ベクトルの始点、終点の超音波発信時刻データとともに格納する。手順305−05、305−06の処理が、浮揚体の移動検出処理である。   Next, proceeding to step 305-06, the correlation between the intensity distribution obtained this time and the intensity distribution stored previously by the preset data collection interval is performed, and the position where the correlation intensity is maximized is obtained, The movement vector of the body 3 is calculated and stored together with the ultrasonic transmission time data of the start point and end point of the movement vector. The processes of procedures 305-05 and 305-06 are the floating body movement detection process.

次に、305−07に進んで流速測定処理を行う。深さレベルiは、初期値1にセットされている。手順305−03で、種々の深さで水中に浮遊している反射体の反射球面波の強度パターン、位相パターンの測定データが反射体の深さレベルとともに格納されており、また、定められた複数の深さにある反射体による反射球面波の強度パターンと位相パターンはデータベースとしてメモリに格納されている。   Next, it progresses to 305-07 and performs a flow velocity measurement process. The depth level i is set to an initial value of 1. In step 305-03, the reflected spherical wave intensity pattern and phase pattern measurement data of the reflector floating in the water at various depths are stored together with the depth level of the reflector. The intensity pattern and phase pattern of the reflected spherical wave by the reflectors at a plurality of depths are stored in the memory as a database.

まず、305−07で、前記格納されている測定データのなかの最新のもののうち、深さレベル(d=i)の反射体の反射球面波の強度パターン、位相パターン及び、データベースとしてメモリに格納されている同じ深さレベルの反射体による反射球面波の強度パターンと位相パターンを取り出し、測定データの強度パターンとデータベース内強度パターン、測定データの位相パターンとデータベース内位相パターン間の相関演算をそれぞれ行い、相関強度のピーク位置を検出して反射体の最新位置を求める。このとき、強度パターンの相関の強度のピークから算出した反射体の位置と位相パターンの相関の強度のピークから算出した反射体の位置との平均の位置を、反射体の位置とする。次いで、前記格納されている測定データのなかの前記データ採取間隔だけ以前のもので、同じ深さレベルの反射体の反射球面波の強度パターン、位相パターン及び、データベースとしてメモリに格納されている同じ深さレベルの反射体による反射球面波の強度パターンと位相パターンを取り出し、測定データの強度パターンとデータベース内強度パターン、測定データの位相パターンとデータベース内位相パターン間の相関演算をそれぞれ行い、相関強度のピーク位置を検出して反射体の前回位置を求める。求めた前記最新位置と前回位置から反射体の移動ベクトルを算出する。   First, at 305-07, the intensity pattern, phase pattern, and database of the reflected spherical wave of the reflector of the depth level (d = i) out of the latest stored measurement data are stored in the memory as a database. The intensity pattern and phase pattern of the reflected spherical wave from the reflector of the same depth level are extracted, and the intensity pattern of the measurement data and the intensity pattern in the database, and the correlation calculation between the phase pattern of the measurement data and the phase pattern in the database, respectively And the peak position of the correlation strength is detected to obtain the latest position of the reflector. At this time, an average position of the position of the reflector calculated from the intensity peak of the correlation of the intensity pattern and the position of the reflector calculated from the intensity peak of the correlation of the phase pattern is set as the position of the reflector. Next, the intensity pattern and phase pattern of the reflected spherical wave of the reflector of the same depth level, which is the previous one of the stored measurement data, and the same stored in the memory as a database. Extract the intensity pattern and phase pattern of the reflected spherical wave by the reflector at the depth level, and perform correlation calculation between the intensity pattern of the measurement data and the intensity pattern in the database, and between the phase pattern of the measurement data and the phase pattern in the database The previous position of the reflector is obtained by detecting the peak position of the reflector. A movement vector of the reflector is calculated from the obtained latest position and previous position.

算出された反射体の移動ベクトルには、浮揚体3の移動ベクトル成分が入っているため、差を取って真の流速ベクトルを求め、移動ベクトル終点の時刻データとともにメモリに格納する(305−08)。この操作を、深さレベルdについて深さレベル1から設定された最大深さレベルnに達するまで繰り返す(305−09、305−10)。   Since the calculated movement vector of the reflector includes the movement vector component of the levitation body 3, a true flow velocity vector is obtained by taking the difference and stored in the memory together with the time data of the movement vector end point (305-08). ). This operation is repeated for the depth level d until reaching the maximum depth level n set from the depth level 1 (305-09, 305-10).

305−09で、最大深さレベルnでの流速算出終了が確認されたら、さらに次の浮揚体位置での測定を行なうのかどうか、言い換えると超音波発信停止スイッチがオフかオンかを判断する(305−11)。超音波発信停止スイッチがオンのとき、計測終了ではないと判断して深さレベルを指定するパラメータiを初期値i=1として手順305−02に戻り、所定の時間間隔で上記手順を繰り返す。   When it is confirmed in 305-09 that the calculation of the flow velocity at the maximum depth level n has been completed, it is determined whether or not to perform measurement at the next floating body position, in other words, whether the ultrasonic wave transmission stop switch is off or on ( 305-11). When the ultrasonic wave transmission stop switch is on, it is determined that the measurement is not finished, the parameter i for designating the depth level is set to the initial value i = 1, and the procedure returns to the procedure 305-02, and the above procedure is repeated at a predetermined time interval.

超音波発信停止スイッチがオフのとき、計測終了と判断して手順305−12に進み、浮揚体の移動ベクトルの和として浮揚体位置を求め、流速分布、河床形状、流量などを画面表示手段に表示する。   When the ultrasonic transmission stop switch is off, it is determined that the measurement is finished, and the process proceeds to step 305-12. The floating body position is obtained as the sum of the floating body moving vectors, and the flow velocity distribution, river bed shape, flow rate, etc. indicate.

浮揚体の移動の軌跡は、前記移動ベクトルの和として求められ、各移動ベクトルの終点位置が、水深、流速の測定点となる。したがって、河床形状は、各移動ベクトルの終点位置の時刻データと同じ時刻データを持つ水深データを各測定点の水深データとして前記浮揚体の移動の軌跡上に配列することで得られる。また、流量は、各深さレベルでの流速が検出されているため、各移動ベクトルの終点位置の時刻データと同じ時刻データを持つ流速データをメモリから順に取り出し、これらを積分計算して求める。   The trajectory of the movement of the levitation body is obtained as the sum of the movement vectors, and the end point position of each movement vector becomes the measurement point of the water depth and the flow velocity. Therefore, the riverbed shape can be obtained by arranging the water depth data having the same time data as the time data of the end point position of each movement vector on the trajectory of the movement of the floating body as the water depth data of each measurement point. Further, since the flow rate at each depth level is detected, the flow rate is obtained by sequentially extracting flow rate data having the same time data as the time data at the end point position of each movement vector, and integrating these.

本実施の形態によれば、河川水面を移動する浮揚体、すなわち超音波を発信する測定装置の移動速度・方向を精度よく求めることができる。また、超音波ビームを用いて水深、流速分布を検出し、前記測定装置の移動速度・方向から測定装置移動軌跡の河床形状、流速分布、流量を精度良く算出できる。
<第2の実施の形態>
本発明の第2の実施の形態が前記第1の実施の形態と異なるのは、超音波送受信器301は、発信する超音波ビームの指向性を制御することが可能になっている点である。他の構成は前記第1の実施の形態と同じであるので、説明を省略する。
According to the present embodiment, it is possible to accurately obtain the moving speed and direction of a floating body that moves on the river surface, that is, a measuring device that transmits ultrasonic waves. In addition, it is possible to detect water depth and flow velocity distribution using an ultrasonic beam, and to accurately calculate the riverbed shape, flow velocity distribution, and flow rate of the measurement device movement locus from the movement speed and direction of the measurement device.
<Second Embodiment>
The second embodiment of the present invention is different from the first embodiment in that the ultrasonic transceiver 301 can control the directivity of the transmitted ultrasonic beam. . Since other configurations are the same as those of the first embodiment, description thereof is omitted.

第1の実施の形態の説明で述べたように、2次元超音波検出器303で検出される河底からの反射波の強度むらの平均的な変化幅は、超音波ビームの指向角αと水中における超音波の波長できまり、深さに依存しない。本実施の形態では、流れの方向と垂直な方向に浮揚体3を移動させる場合、浮揚体3の移動方向にはビーム角度を狭く、流れの方向にはビーム角度が広くなるように超音波送受信器301を制御する。ビーム指向角度の制御には、超音波レンズの使用、アレイ状センサの使用など公知の方法が使える。この制御により、河底反射によるスペックルパターン(強度分布)の平均的な変化幅は浮揚体の移動方向で大きくなり、移動量が大きい場合有効である。また流れの方向にビーム角が広いため、水中に浮遊する反射体の移動幅も大きくすることができ、測定範囲が拡大する。   As described in the description of the first embodiment, the average change width of the intensity unevenness of the reflected wave from the riverbed detected by the two-dimensional ultrasonic detector 303 is the directivity angle α of the ultrasonic beam. It depends on the wavelength of the ultrasonic wave in water and does not depend on the depth. In the present embodiment, when the levitation body 3 is moved in a direction perpendicular to the flow direction, ultrasonic waves are transmitted and received so that the beam angle is narrow in the movement direction of the levitation body 3 and the beam angle is wide in the flow direction. The device 301 is controlled. For controlling the beam directing angle, a known method such as the use of an ultrasonic lens or an array sensor can be used. By this control, the average change width of the speckle pattern (intensity distribution) due to riverbed reflection increases in the moving direction of the levitating body, which is effective when the moving amount is large. In addition, since the beam angle is wide in the flow direction, the moving width of the reflector floating in water can be increased, and the measurement range is expanded.

本実施の形態によれば、前記第1の実施の形態におけると同様の効果が得られるとともに、測定装置を搭載した浮揚体の移動速度が大きい場合や、流速が大きい場合の対応が可能になるという効果がある。
<第3の実施の形態>
本発明の第3の実施の形態が前記第1の実施の形態と異なるのは、浮揚体3に傾きセンサを設置し、超音波放射時の浮揚体の傾きを求めて河底までの距離、反射体の位置の補正を行うようにした点である。他の構成は前記第1の実施の形態と同じであるので、説明を省略する。
According to the present embodiment, the same effect as in the first embodiment can be obtained, and it is possible to cope with a case where the moving speed of the floating body on which the measuring device is mounted is large or a case where the flow velocity is large. There is an effect.
<Third Embodiment>
The third embodiment of the present invention differs from the first embodiment in that a tilt sensor is installed in the levitation body 3, and the distance to the river bed is obtained by obtaining the tilt of the levitation body during ultrasonic radiation. The point is that the position of the reflector is corrected. Since other configurations are the same as those of the first embodiment, description thereof is omitted.

超音波の放射時、浮揚体3が傾くと河底や反射体深さまでの超音波経路が長くなる。また、超音波照射ごとに浮揚体3の傾きが異なると反射体位置が見かけ上ずれ、流速測定の誤差となる。このため、浮揚体3に傾きセンサを設置し、超音波放射時の浮揚体3の傾きを求めて河底までの距離、反射体の位置の補正を行う。   When the levitation body 3 is tilted during the emission of ultrasonic waves, the ultrasonic path to the riverbed or the depth of the reflector becomes longer. Moreover, if the inclination of the levitation body 3 is different for each ultrasonic irradiation, the reflector position is apparently shifted, resulting in an error in flow velocity measurement. For this reason, an inclination sensor is installed in the levitation body 3, the inclination of the levitation body 3 at the time of ultrasonic radiation is obtained, and the distance to the riverbed and the position of the reflector are corrected.

本実施の形態によれば、前記第1の実施の形態におけると同様の効果が得られるとともに、河床形状、流速分布の計測がより正確になる効果がある。   According to the present embodiment, the same effect as in the first embodiment can be obtained, and the riverbed shape and the flow velocity distribution can be more accurately measured.

なお、河床形状の測定のみを行なう場合には、流速を測定する部分は除去し形状測定部分のみを使用してもよい。   When only the riverbed shape measurement is performed, the portion for measuring the flow velocity may be removed and only the shape measurement portion may be used.

河川データ測定の基本的概念を示す斜視図である。It is a perspective view which shows the basic concept of river data measurement. 本発明の第1の実施の形態に係る河川データ測定装置の要部構成を示すブロック図である。It is a block diagram which shows the principal part structure of the river data measuring apparatus which concerns on the 1st Embodiment of this invention. 本発明の第1の実施の形態による流速分布の測定を模式的に説明する概念図である。It is a conceptual diagram which illustrates typically the measurement of the flow-velocity distribution by the 1st Embodiment of this invention. 本発明の第1の実施の形態による浮揚体の移動測定を模式的に説明する概念図である。It is a conceptual diagram which illustrates typically the movement measurement of the floating body by the 1st Embodiment of this invention. 本発明の第1の実施の形態による浮揚体の移動測定における相関演算を説明する概念図である。It is a conceptual diagram explaining the correlation calculation in the movement measurement of the levitating body by the 1st Embodiment of this invention. 本発明の第1の実施の形態における超音波の波形、強度、位相について説明する概念図である。It is a conceptual diagram explaining the waveform, intensity | strength, and phase of an ultrasonic wave in the 1st Embodiment of this invention. 本発明による河床形状、流速分布、流量の測定装置に関し、河底などの波形関係と強度の検出を説明する図である。It is a figure explaining the waveform relationship and intensity | strength detection, such as a river bed, regarding the riverbed shape by this invention, flow velocity distribution, and the measuring apparatus of flow volume. 図2に示す実施の形態の部分の詳細を示すブロック図である。It is a block diagram which shows the detail of the part of embodiment shown in FIG. 図2に示す実施の形態の他の部分の詳細を示すブロック図である。It is a block diagram which shows the detail of the other part of embodiment shown in FIG. 本発明の第1の実施の形態の信号解析装置305の処理の流れを示す手順図である。It is a procedure figure showing the flow of processing of signal analysis device 305 of a 1st embodiment of the present invention.

符号の説明Explanation of symbols

1 河川
2 測定線
3 浮揚体
301 超音波送受信器
302 送受信回路
302−1 スイッチ素子
302−2 増幅器
302−3 フィルターアンプ
302−4 検波器
302−5 レベル判定器
303 2次元超音波検出器
304 信号処理装置
304i 信号処理装置304のひとつの処理回路
304i−1 フィルターアンプ
304i−2 検波器
304i−3 位相検出器
304i−4 レベル判定器
304i−5、6 スイッチ素子
304i−7、8 アナログ−デジタル変換器
305 信号解析装置
DESCRIPTION OF SYMBOLS 1 River 2 Measuring line 3 Floating body 301 Ultrasonic transmitter / receiver 302 Transmission / reception circuit 302-1 Switch element 302-2 Amplifier 302-3 Filter amplifier 302-4 Detector 302-5 Level judgment device 303 Two-dimensional ultrasonic detector 304 Signal Processing device 304i One processing circuit of the signal processing device 304i-1 Filter amplifier 304i-2 Detector 304i-3 Phase detector 304i-4 Level determination device 304i-5, 6 Switch element 304i-7, 8 Analog-digital conversion 305 Signal Analyzer

Claims (9)

測定線に沿って移動する浮揚体に搭載され、河底に向かって発信される超音波によって河川のデータを測定する測定装置であって、河底方向に所定の時間間隔で超音波を発信する超音波送信手段と、前記超音波送信手段から発信され、河底で散乱・反射された超音波を受信する2次元状に配置された2次元超音波受信手段と、前記超音波送信手段から時間をおいて発信された2回の超音波それぞれに対応する前記2次元超音波受信手段の出力に基づいて河底からの反射波の強度分布を前記2次元超音波受信手段の中心を基準位置として求める信号処理手段と、信号処理手段に接続され、信号処理手段から出力される二つの強度分布の間の相関の強度を、前記二つの強度分布のうちの先に発信された超音波に対応するものを他方に近づけつつ移動距離の関数として求め、相関の強度が最大となる移動距離及びそのときの移動方向を前記2回の超音波発信の間の超音波送信手段の移動ベクトルとして出力する信号解析手段と、前記移動ベクトルが格納される記憶手段とを有してなる河川データ測定装置。 A measuring device that is mounted on a levitating body that moves along a measurement line and measures river data using ultrasonic waves transmitted toward the bottom of the river, and emits ultrasonic waves at predetermined time intervals in the direction of the bottom of the river. Ultrasonic transmission means, two-dimensional ultrasonic reception means arranged in a two-dimensional manner for receiving ultrasonic waves transmitted from the ultrasonic transmission means and scattered / reflected at the river bed, and time from the ultrasonic transmission means reference position the center of the two-dimensional ultrasonic receiving means the intensity distribution of the reflected waves from the riverbed on the basis of the output of two of the two-dimensional ultrasonic receiving means corresponding to the ultrasonic each originating at a and determined Mel signal processing means as a, is connected to the signal processing means, the strength of the correlation between the two intensity distribution output from the signal processing means, originating earlier of said two intensity distribution Move the one corresponding to the ultrasonic wave close to the other. A signal analysis means for obtaining a movement distance and a moving direction at which the correlation intensity is maximized as a function of the distance and a movement direction of the ultrasonic transmission means between the two ultrasonic transmissions; A river data measuring device comprising storage means for storing the data. 請求項1記載の河川データ測定装置において、前記超音波送信手段に代えて、河底方向に所定の時間間隔で超音波を発信するとともに反射波を受信する超音波送受信手段と、前記超音波送受信手段に接続され、超音波発信から河底の反射波受信までの時間と河底からの反射波の強度を測定する河底反射測定手段とが設けられ、前記信号解析手段は、前記河底反射測定手段の出力を入力として水深を算出して前記記憶手段に格納し、さらに、前記移動ベクトルを順次結合して前記超音波送受信手段の移動の軌跡を生成するとともに、前記移動の軌跡を構成する各移動ベクトルの終点位置に対応する水深のデータを前記記憶手段から読み出して前記終点位置のデータに組み合わせて前記記憶手段に格納するよう構成されていることを特徴とする河川データ測定装置。   2. The river data measuring apparatus according to claim 1, wherein, in place of the ultrasonic transmission means, ultrasonic transmission / reception means for transmitting ultrasonic waves at a predetermined time interval in the river bottom direction and receiving reflected waves, and the ultrasonic transmission / reception. And a bottom reflection measuring means for measuring the time from ultrasonic transmission to reception of reflected waves at the bottom of the river and the intensity of reflected waves from the bottom of the river, and the signal analyzing means is provided with the bottom reflection. The water depth is calculated by using the output of the measuring means as input and stored in the storage means. Further, the movement vectors of the ultrasonic transmitting / receiving means are generated by sequentially combining the movement vectors, and the movement locus is configured. A water depth data corresponding to the end point position of each movement vector is read from the storage means and combined with the end position data and stored in the storage means. Data measuring device. 請求項1記載の河川データ測定装置において、前記超音波送信手段に代えて、河底方向に所定の時間間隔で超音波を発信するとともに反射波を受信する超音波送受信手段と、前記超音波送受信手段に接続され、超音波発信から河底の反射波受信までの時間と河底からの反射波の強度を測定する河底反射測定手段とが設けられ、前記2次元超音波受信手段は河底で散乱・反射された超音波に加え、水中に浮遊する水中反射体で反射された反射波を受信するよう構成されているとともに前記信号処理手段は、超音波発信から水中反射体で反射された反射波を検出するまでの時間と該反射波の強度分布及び位相分布をも測定するよう構成され、前記記憶手段には、それぞれ異なる深さでかつ前記超音波送受信手段との水平方向相対位置が所定の位置にある水中反射体によって前記超音波送受信手段から発信された超音波が反射されたときの反射波の強度分布及び位相分布が基準強度分布及び基準位相分布として格納され、前記信号解析手段は、前記河底反射測定手段の出力を入力として水深を算出して前記記憶手段に格納し、前記信号処理手段の出力に基づいて水中反射体が浮遊している深さを検出し、前記検出された深さの水中反射体の反射波の基準強度分布及び基準位相分布と前記測定された反射波の強度分布及び位相分布のいずれか若しくは双方の相関の強度を、測定された分布を基準分布に近づけつつ移動距離の関数として求め、相関の強度が最大となる移動距離及びそのときの移動方向に基づいて当該水中反射体の超音波送受信手段に対する水平方向位置を求め、所定の時間間隔で検出された水中反射体の超音波送受信手段に対する水平方向相対位置及び前記算出された超音波送受信手段の移動ベクトルに基づいて、前記移動ベクトルの終点位置における、前記検出された深さの流速を求めるように構成されていることを特徴とする河川データ測定装置。   2. The river data measuring apparatus according to claim 1, wherein, in place of the ultrasonic transmission means, ultrasonic transmission / reception means for transmitting ultrasonic waves at a predetermined time interval in the river bottom direction and receiving reflected waves, and the ultrasonic transmission / reception. And a bottom reflection measuring means for measuring the time from transmission of ultrasonic waves to reception of reflected waves from the bottom of the river and the intensity of reflected waves from the bottom of the river. In addition to the ultrasonic waves scattered / reflected at, the reflected wave reflected by the underwater reflector floating in the water is received and the signal processing means is reflected by the underwater reflector from the ultrasonic transmission The time until the reflected wave is detected and the intensity distribution and phase distribution of the reflected wave are also measured, and the storage means has different relative depths in the horizontal direction with respect to the ultrasonic transmission / reception means. In place The intensity distribution and phase distribution of the reflected wave when the ultrasonic wave transmitted from the ultrasonic transmission / reception means is reflected by a certain underwater reflector are stored as a reference intensity distribution and a reference phase distribution, and the signal analysis means The water depth is calculated using the output of the bottom reflection measuring means as an input and stored in the storage means, and the depth at which the underwater reflector is floating is detected based on the output of the signal processing means, and the detected depth The intensity of the correlation between the reference intensity distribution and reference phase distribution of the reflected wave of the underwater reflector and the intensity distribution and / or phase distribution of the measured reflected wave is moved while bringing the measured distribution close to the reference distribution. As a function of distance, the horizontal position of the underwater reflector with respect to the ultrasonic transmission / reception means is obtained based on the moving distance and the moving direction at which the correlation intensity is maximum, and for a predetermined time Based on the horizontal relative position of the underwater reflector to the ultrasonic transmission / reception means detected in step 1 and the calculated movement vector of the ultrasonic transmission / reception means, the flow velocity of the detected depth at the end position of the movement vector is calculated. A river data measuring device, characterized by being configured to obtain. 請求項3記載の河川データ測定装置において、前記信号解析手段は、前記移動ベクトルを順次結合して前記超音波送受信手段の移動の軌跡を生成し、前記移動の軌跡を構成する各移動ベクトルの終点位置に対応する水深のデータを前記記憶手段から読み出して前記終点位置のデータに組み合わせて前記記憶手段に格納するよう構成されていることを特徴とする河川データ測定装置。   4. The river data measuring apparatus according to claim 3, wherein the signal analysis means sequentially combines the movement vectors to generate a movement trajectory of the ultrasonic transmission / reception means, and an end point of each movement vector constituting the movement trajectory. A river data measuring apparatus configured to read out water depth data corresponding to a position from the storage means and combine the data with the data of the end point position and store the data in the storage means. 請求項1乃至4のいずれかに記載の河川データ測定装置において、前記超音波送信手段若しくは超音波送受信手段は、発信する超音波の指向特性を制御する手段を備えてなることを特徴とする河川データ測定装置。   5. The river data measuring apparatus according to claim 1, wherein the ultrasonic transmission means or the ultrasonic transmission / reception means includes means for controlling directivity characteristics of ultrasonic waves to be transmitted. Data measuring device. 請求項2乃至5のいずれかに記載の河川データ測定装置において、前記浮揚体の傾きを検出する傾き検出手段を備え、前記信号解析手段は、傾き検出手段の出力に基づいて算出する水深を補正するように構成されていることを特徴とする河川データ測定装置。   6. The river data measuring apparatus according to claim 2, further comprising an inclination detection unit that detects an inclination of the floating body, wherein the signal analysis unit corrects a water depth calculated based on an output of the inclination detection unit. A river data measuring device characterized by being configured to perform. 測定線に沿って移動する浮揚体に搭載された測定装置から河底に向かって発信される超音波によって河川のデータを測定する河川データ測定方法であって、河底方向に所定の時間間隔で超音波を発信し、河底で散乱・反射された超音波を2次元状に配置された2次元超音波受信手段で受信し、時間をおいて発信された2回の超音波それぞれに対応する河底からの反射波の強度分布を前記2次元超音波受信手段の中心を基準位置として求め、前記二つの強度分布の間の相関の強度を、前記二つの強度分布のうちの先に発信された超音波に対応するものを他方に近づけつつ移動距離の関数として求め、相関の強度が最大となる移動距離及びそのときの移動方向を前記2回の超音波発信の間の超音波発信源の移動距離及び移動方向として出力する手順を有してなる河川データ測定方法。 A river data measurement method for measuring river data by ultrasonic waves transmitted toward a river bed from a measurement device mounted on a levitating body that moves along a measurement line, at predetermined time intervals in the river bed direction. The ultrasonic wave is transmitted, and the ultrasonic wave scattered and reflected at the riverbed is received by the two-dimensional ultrasonic wave receiving means arranged in a two-dimensional manner, corresponding to each of the two ultrasonic waves transmitted with time. the intensity distribution of the reflected waves from the riverbed as a reference position the center of the two-dimensional ultrasonic receiving means determined Me, the strength of the correlation between the two intensity distribution among the two intensity distribution The one corresponding to the ultrasonic wave transmitted earlier is obtained as a function of the moving distance while approaching the other, and the moving distance and the moving direction at which the intensity of the correlation is the maximum between the two ultrasonic transmissions are determined. Hands that output as the moving distance and moving direction of the ultrasonic transmission source River data measurement method comprising a. 請求項7記載の河川データ測定方法において、超音波発信から河底の反射波受信までの時間と河底からの反射波の強度を測定し、超音波発信から河底の反射波受信までの時間を入力として水深を算出して記憶手段に格納し、さらに、前記移動ベクトルを順次結合して測定装置の移動の軌跡を生成するとともに、前記移動の軌跡を構成する各移動ベクトルの終点位置に対応する水深のデータを前記記憶手段から読み出して前記終点位置のデータに組み合わせて前記記憶手段に格納する手順を有してなる河川データ測定方法。   8. The river data measurement method according to claim 7, wherein the time from ultrasonic transmission to reception of reflected waves on the riverbed and the intensity of reflected waves from the riverbed are measured, and the time from transmission of ultrasonic waves to reception of reflected waves on the riverbed. The water depth is calculated and stored in the storage means, and the movement vectors are sequentially combined to generate a movement trajectory of the measuring device and correspond to the end point position of each movement vector constituting the movement trajectory. A river data measuring method comprising a step of reading out data of water depth to be read from the storage means and combining the data with the data of the end point position and storing the data in the storage means. 請求項7記載の河川データ測定方法において、超音波発信から河底の反射波受信までの時間と河底からの反射波の強度を測定し、河底で散乱・反射された超音波に加え、水中に浮遊する水中反射体で反射された反射波を前記2次元超音波受信手段で受信するとともに、超音波発信から水中反射体で反射された反射波を検出するまでの時間と該反射波の強度分布及び位相分布をも測定し、それぞれ異なる深さでかつ超音波発信源との水平方向相対位置が所定の位置にある水中反射体によって前記超音波送受信手段から発信された超音波が反射されたときの反射波の強度分布及び位相分布を基準強度分布及び基準位相分布として前記記憶手段に格納され、前記超音波発信から河底の反射波受信までの時間を入力として水深を算出して前記記憶手段に格納し、前記超音波発信から水中反射体の反射波受信までの時間に基づいて水中反射体が浮遊している深さを検出し、前記検出された深さの水中反射体の反射波の基準強度分布及び基準位相分布と前記測定された反射波の強度分布及び位相分布のいずれか若しくは双方の相関の強度を、測定された分布を基準分布に近づけつつ移動距離の関数として求め、相関の強度が最大となる移動距離及びそのときの移動方向に基づいて当該水中反射体の超音波発信源に対する水平方向位置を求め、所定の時間間隔で検出された水中反射体の超音波発信源に対する水平方向相対位置及び前記算出された超音波発信源の移動ベクトルに基づいて、前記移動ベクトルの終点位置における、前記検出された深さの流速を求める手順を有してなることを特徴とする河川データ測定方法。   The river data measurement method according to claim 7, wherein the time from the transmission of ultrasonic waves to the reception of reflected waves at the bottom of the river and the intensity of reflected waves from the bottom of the river are measured, and in addition to the ultrasonic waves scattered and reflected at the bottom of the river, The reflected wave reflected by the underwater reflector floating in the water is received by the two-dimensional ultrasonic receiving means, and the time from the ultrasonic transmission until the reflected wave reflected by the underwater reflector is detected and the reflected wave The intensity distribution and the phase distribution are also measured, and the ultrasonic wave transmitted from the ultrasonic wave transmitting / receiving means is reflected by an underwater reflector having a different depth and a horizontal relative position to the ultrasonic wave transmission source at a predetermined position. The intensity distribution and the phase distribution of the reflected wave at that time are stored in the storage means as a reference intensity distribution and a reference phase distribution, and the water depth is calculated by inputting the time from the transmission of the ultrasonic wave to the reception of the reflected wave of the riverbed as the input. Storage means Storing and detecting the floating depth of the underwater reflector based on the time from the transmission of the ultrasonic wave to the reception of the reflected wave of the underwater reflector, and the reference of the reflected wave of the underwater reflector of the detected depth The intensity of correlation between the intensity distribution and the reference phase distribution and / or the intensity distribution and / or phase distribution of the measured reflected wave is obtained as a function of the moving distance while bringing the measured distribution close to the reference distribution. The horizontal position of the underwater reflector with respect to the ultrasonic transmission source is obtained based on the moving distance at which the maximum is obtained and the moving direction at that time, and the horizontal direction with respect to the ultrasonic transmission source of the underwater reflector detected at a predetermined time interval. A step of obtaining a flow velocity of the detected depth at an end point position of the movement vector based on the relative position and the calculated movement vector of the ultrasonic wave transmission source; River data measurement method.
JP2005208406A 2005-07-19 2005-07-19 River data measurement method and apparatus Expired - Fee Related JP4633565B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005208406A JP4633565B2 (en) 2005-07-19 2005-07-19 River data measurement method and apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005208406A JP4633565B2 (en) 2005-07-19 2005-07-19 River data measurement method and apparatus

Publications (2)

Publication Number Publication Date
JP2007024715A JP2007024715A (en) 2007-02-01
JP4633565B2 true JP4633565B2 (en) 2011-02-16

Family

ID=37785672

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005208406A Expired - Fee Related JP4633565B2 (en) 2005-07-19 2005-07-19 River data measurement method and apparatus

Country Status (1)

Country Link
JP (1) JP4633565B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110608720A (en) * 2019-10-16 2019-12-24 北京交通大学 Riverbed elevation detection method and system based on ultrasonic reflection waveform
CN110645964A (en) * 2019-09-30 2020-01-03 浙江海洋大学 Device and method for evaluating rising of offshore sea level based on river runoff

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1757956A1 (en) * 2005-08-24 2007-02-28 Leica Geosystems AG Multiple target capable ranging method according to the phase measuring method
KR20130137005A (en) 2010-11-16 2013-12-13 퀄컴 인코포레이티드 System and method for object position estimation based on ultrasonic reflected signals
CN106525192A (en) * 2016-11-16 2017-03-22 北京中船信息科技有限公司 Real-time river section measurement device and method
JP2019035677A (en) * 2017-08-17 2019-03-07 株式会社 拓和 Underwater distance observation device
JP7054507B2 (en) * 2017-11-17 2022-04-14 国立大学法人金沢大学 Measuring device, mobile body, and measuring method
CN109444894B (en) * 2018-12-27 2023-11-28 镇江市高等专科学校 Road ponding depth rapid movement measuring device and method
CN110068822B (en) * 2019-06-11 2024-05-14 水利部交通运输部国家能源局南京水利科学研究院 Two-dimensional accurate positioning auxiliary device of portable sonar sounder
CN111596376B (en) * 2020-05-07 2023-02-10 中国石油天然气股份有限公司 Method and device for determining river form conversion of river facies stratum
CN115930918B (en) * 2023-03-10 2023-05-30 山东省地质矿产勘查开发局八〇一水文地质工程地质大队(山东省地矿工程勘察院) Intelligent observation equipment for hydrologic dynamics and use method
CN117250366B (en) * 2023-11-17 2024-02-02 威海多鱼海洋科技有限公司 Underwater monitoring device and method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS50142072A (en) * 1974-05-02 1975-11-15
JPS5582067A (en) * 1978-11-06 1980-06-20 Gen Electric Correlating soner apparatus for measuring speed
JPS57149908A (en) * 1981-03-11 1982-09-16 Ishikawajima Harima Heavy Ind Co Ltd Sounding machine
JPS60242383A (en) * 1984-05-16 1985-12-02 Furuno Electric Co Ltd Three-dimensional underwater detecting device
JPS61262674A (en) * 1985-05-16 1986-11-20 Nec Corp Apparatus for measuring position in water
JP2003315455A (en) * 2002-04-18 2003-11-06 Koden Electronics Co Ltd Ultrasonic surveying apparatus

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS50142072A (en) * 1974-05-02 1975-11-15
JPS5582067A (en) * 1978-11-06 1980-06-20 Gen Electric Correlating soner apparatus for measuring speed
JPS57149908A (en) * 1981-03-11 1982-09-16 Ishikawajima Harima Heavy Ind Co Ltd Sounding machine
JPS60242383A (en) * 1984-05-16 1985-12-02 Furuno Electric Co Ltd Three-dimensional underwater detecting device
JPS61262674A (en) * 1985-05-16 1986-11-20 Nec Corp Apparatus for measuring position in water
JP2003315455A (en) * 2002-04-18 2003-11-06 Koden Electronics Co Ltd Ultrasonic surveying apparatus

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110645964A (en) * 2019-09-30 2020-01-03 浙江海洋大学 Device and method for evaluating rising of offshore sea level based on river runoff
CN110608720A (en) * 2019-10-16 2019-12-24 北京交通大学 Riverbed elevation detection method and system based on ultrasonic reflection waveform

Also Published As

Publication number Publication date
JP2007024715A (en) 2007-02-01

Similar Documents

Publication Publication Date Title
JP4633565B2 (en) River data measurement method and apparatus
JP4415192B2 (en) Riverbed measuring device
JP5389267B2 (en) Method and apparatus for measuring seabed contours
US20100142326A1 (en) Method For Identifying Gas Leaks Using A Stationary Seabed Placed Steered Beamformed Acoustic Antenna With Active Narrow Beam Transmitter Interrogation Capability
WO2017177967A1 (en) Underwater detection system and underwater detection method
KR101899079B1 (en) System for detecting location of underground objects
RU2343502C2 (en) Method and system of positional analysis of object under observation by depth in aqueous medium
WO2020097390A1 (en) Horizontal acoustic sediment and current profiler apparatus and methods
RU2515957C1 (en) Complex for ultrasonic control of products and optical measurement device of complex
JP2010112895A (en) Target object detector, target object detection control program, and target object detection method
CN103471998B (en) Thermoplastic material reflection and transmission coefficients laser measurement system
JP2009264965A (en) Underwater sailing body and obstacle detection apparatus
RU2350983C2 (en) Method for determination of object submersion depth
US20220236437A1 (en) Method and system for determining top and bottom depth of an under water mud layer
KR101331333B1 (en) Method and device for measuring a profile of the ground
RU2541435C1 (en) Method of determining iceberg immersion
RU2125278C1 (en) Method measuring distance to controlled object ( its versions )
RU2563317C1 (en) Method of acoustic monitoring of offshore parameters variability
EA035249B1 (en) Detection system and method to check the position of a pipeline in a bed of a body of water
RU2452977C1 (en) Method of measuring distance to monitored facility
JP2010190682A (en) Method and apparatus for exploring and classifying object in the bottom of water
RU2376612C1 (en) Method of hydrometeorological monitoring water body of sea test site and device to this end
Wanis Design and applications of a vertical beam in acoustic Doppler current profilers
RU2625041C1 (en) Method for measuring object immersion depth
RU2280266C2 (en) Method of detection of objects lying in thickness of bottom ground

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070717

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20100520

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100525

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100726

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20101102

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20101117

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131126

Year of fee payment: 3

LAPS Cancellation because of no payment of annual fees