CN101532825A - Method for measuring thickness of sea surface spilled oil film based on differential laser triangulation method - Google Patents

Method for measuring thickness of sea surface spilled oil film based on differential laser triangulation method Download PDF

Info

Publication number
CN101532825A
CN101532825A CN200910068627A CN200910068627A CN101532825A CN 101532825 A CN101532825 A CN 101532825A CN 200910068627 A CN200910068627 A CN 200910068627A CN 200910068627 A CN200910068627 A CN 200910068627A CN 101532825 A CN101532825 A CN 101532825A
Authority
CN
China
Prior art keywords
msubsup
mrow
msub
oil film
laser
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN200910068627A
Other languages
Chinese (zh)
Other versions
CN101532825B (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN200910068627XA priority Critical patent/CN101532825B/en
Publication of CN101532825A publication Critical patent/CN101532825A/en
Application granted granted Critical
Publication of CN101532825B publication Critical patent/CN101532825B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention relates to a method for measuring thickness of sea surface spilled oil film based on differential laser triangulation. The method uses dual laser triangulation ranging system, and laser beams emitted by two lasers are respectively irradiated to the upper and the lower surfaces of an oil film to be measured in any vertical direction; the off-plane displacements of the upper and the lower surfaces can be obtained simultaneously according to the size of image drift and configuration parameters of imaging system; the thickness of the measured oil film can be determined by subtracting; and the thickness distribution of oil film in spilled oil area can be obtained by multi-point sampling. The non-contact measuring method has the characteristics of simple theory, convenient measurement, low cost, high precision, small volume and the like, and can be applied to dynamic on-line measurement of oil film thickness at the surface of aqueous solution.

Description

Sea surface oil spill oil film thickness measuring method based on differential laser triangulation method
[ technical field ]:
the invention belongs to the technical field of marine environment three-dimensional monitoring, and particularly relates to a technology for measuring the thickness of a sea surface spilled oil film.
[ background Art ] A method of:
the measurement of the thickness of the oil film on the sea surface is one of indispensable works of marine environment management, and is an important problem in the fields of environmental protection, chemical analysis, oil film dynamics research and the like. The ocean monitoring department and experts have always sought a method for accurately measuring the thickness of the water surface oil film, but until now, a credible method for measuring the thickness of the water surface oil film cannot be found. The oil film thickness provides effective and reasonable guidance for oil spill treatment (such as correct and reasonable use of a dispersing agent, field incineration and the like), provides important scientific basis for understanding the basic theory of oil layer extended dynamics, and has important significance for weather, climate and marine environment treatment. The following methods have been proposed for measuring the thickness of different oil and oil films on the surface of aqueous solutions:
the method can simultaneously realize oil product identification and oil film thickness measurement by analyzing the fluorescence life of different oil products and the spectral intensity under different oil film thicknesses, and is suitable for measuring the thickness of a thin oil film.
The microwave radiometer obtains the oil layer thickness information by comparing the measured brightness temperature, but the method needs to determine some environmental parameters and oil property parameters in advance.
In the multi-beam interference method, the oil film forms two parallel plane layers of oil/gas and oil/water on the sea surface. Light is incident on the oil-gas interface, part of the light is reflected to the air, and part of the light is transmitted to enter the oil layer. The light beam entering the oil layer is reflected/transmitted at the oil-water interface, and the reflected light beam is reflected back to the oil layer through the oil-water interface and then is transmitted back to the air. The light waves in the air are superposed to generate interference, and the thickness of the oil film is calculated according to the interference phase difference.
The method is characterized in that oil quality identification is carried out by analyzing the peak ratio of interference envelopes formed by different oil products on a water surface and an air interface, and the oil film thickness is calculated according to the distance between two groups of interference fringe peaks formed by an oil-gas interface and an oil-water interface. The technology based on the interference method is suitable for measuring the thickness of a transparent oil layer with a flat oil-water interface.
The basic mechanism of the laser acoustics/laser ultrasonic technology is based on laser radar and ultrasonic principle, and the thickness of an oil layer is calculated by utilizing the time of sound pulse transmitted through the oil layer and then returned to an oil-gas interface and the sound wave transmission speed in the oil layer. A laser ultrasonic remote sensing oil film thickness measuring sensor has been developed, and the sensor has high measuring precision and great potential in the field of aerial oil spill remote sensing measurement, but three lasers are needed: a CO2Pulsed lasers, a Nd: YAG laser and a He-Ne laser, and are large in volume.
[ summary of the invention ]:
the invention aims to provide a method for measuring the thickness of a sea surface spilled oil film based on a differential laser triangulation method, which is used for solving the problem of measuring the thickness of the sea surface spilled oil film, filling the blank of the technical research in the aspect of marine oil monitoring in China, promoting the development of the existing oil monitoring level in China and providing a technical basis for the measurement of the thickness of the sea surface spilled oil film, the oil layer extension dynamics and the marine environment science.
The invention firstly proposes that the thickness measuring technology of the laser triangulation method is applied to the measurement of the thickness of the oil film of the spilled oil, the fluctuation of the seawater and the oil film is considered, and a double-optical-path differential system, namely a symmetrical oblique-incidence double-laser triangulation distance measuring system, is adopted. Laser beams emitted by the laser device are converged on the upper surface and the lower surface of the oil film by the converging lens, and an elliptical laser spot is formed at an incidence point on the upper surface and the lower surface of the oil film. The imaging system images the light spot on the light receiver CCD to form an image point. By using a double-laser triangulation system of the upper and lower surfaces, according to the image motion size h and the structural parameter k of the imaging system, the out-of-plane displacement of the upper and lower surfaces at the point can be obtained simultaneously, and the oil film thickness of the measured point can be determined by subtraction. In the optical path system, the laser beam, the optical axis of the imaging lens and the CCD are ensured to be in the same plane, the included angle between the optical axis of the laser beam and the optical axis of the imaging lens is strictly ensured to be 90 degrees, the included angle between the CCD and the optical axis of the imaging lens is 90 degrees, the object distance is ensured to be unchanged after the position of the measured surface is changed, the measured point in a certain depth of field can be ensured to be focused and imaged on the detector, and the measurement precision is ensured.
The specific measuring process of the measuring method comprises the following steps:
first, light path System building
Building a symmetrical oblique double-laser triangular distance measuring system according to the figure 1, wherein each laser triangular distance measuring system comprises a light source part consisting of a laser and a converging lens, and a receiving imaging system consisting of an imaging lens, a detector CCD (charge coupled device) or a position detector PSD (position sensitive detector) and a data acquisition and processing module; like parts are designated by like reference numerals in the drawings;
second, light spot image acquisition
The laser beam emitted by the first laser device is converged and obliquely incident on the upper surface of the oil film to be detected through the converging lens, an elliptical laser spot is formed at an incident point on the upper surface of the oil film (the position of the lens is adjusted, so that the focus is positioned near the upper surface of the oil film to be detected, and the spot has strong brightness and small radius to form a bright point), and the imaging lens in the receiving imaging system images the spot on the CCD; the laser beam emitted by the second laser is converged by the convergent lens and obliquely enters the seawater through the air, the seawater is reflected and refracted at the air-water interface, the refracted light is obliquely incident on the lower surface of the oil film, an elliptical laser spot is also formed at the incident point on the lower surface of the oil film (the position of the lens is adjusted, the focus is positioned near the lower surface of the oil film to be detected, the spot has very strong brightness and very small radius at the moment, a bright spot is formed), the imaging system images the spot on the CCD, and the upper illuminating point and the lower illuminating point are made to correspond to the upper surface and the lower surface of any point of the object to be detected in. And collecting light spot image points of the two systems on the CCD. The reference plane is set as a virtual reference plane, and the selection can be arbitrary theoretically;
third, image displacement measurement
The collected light spot image is processed by adopting filtering technologies such as morphological filtering and the like, the centroid of two image points on the CCD is obtained by utilizing methods such as a centroid method and Hough transformation, and the image displacement h of the upper and lower light spot image points relative to the reference plane light spot image point can be obtained by utilizing a distance formula between the two points1And <math> <mrow> <msub> <mi>h</mi> <mn>2</mn> </msub> <mo>=</mo> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>-</mo> <msubsup> <mi>h</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>.</mo> </mrow></math>
fourthly, establishing a mathematical model for calculating out-of-plane displacement
According to the corresponding displacement h of the light spot on the CCD obtained in the previous step1And h2Establishing a mathematical model for calculating the out-of-plane displacement, and calculating by a data processing module to obtain the out-of-plane displacement d of the reference plane assumed by the distance between the upper surface and the lower surface of the point1And d2(displacement square)To the direction shown, the upward arrow is positive and the downward is negative). Wherein d is adopted for the first optical path system1=k1h1Formula (I) is calculated, wherein k1Cos α/f (u-f). For the second optical path system <math> <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>l</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>l</mi> <mn>2</mn> </msub> <mo>=</mo> <msubsup> <mi>k</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>-</mo> <msubsup> <mi>k</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <msubsup> <mi>h</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>,</mo> </mrow></math> In the formula <math> <mrow> <msubsup> <mi>k</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>cos</mi> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> <mo>/</mo> <mi>f</mi> <mo>,</mo> </mrow></math> <math> <mrow> <msubsup> <mi>k</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <msqrt> <msubsup> <mi>n</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>n</mi> <mn>0</mn> <mn>2</mn> </msubsup> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> </msqrt> <mo>/</mo> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mn>0</mn> </msub> <mo>,</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow></math> u and a are object distances of two imaging systems, alpha is an included angle between the optical axis of the laser beam and the normal line of the measured surface, and theta1Is the incident angle of the laser beam incident on the air-sea interface, n0Is the refractive index of air, n1Is the refractive index of seawater, and f is the focal length of the imaging lens;
fifth, data processing and output
For two out-of-plane displacements d1And d2Subtracting the oil film thickness d ═ d of the point1-d2
The measurement principle of the method of the invention is as follows:
fig. 2 is a basic schematic diagram. In the oblique type shown in fig. 2(a), the laser beam is incident at the measurement point a at an angle α and forms an elliptical spot at that point. On the other side of the normal line at A, the optical axis of the imaging lens forms an angle beta with the normal line of the imaging lens, and the camera captures an image of the light spot and images the light spot at a on the CCD. When the incident light spot moves from the point A to the point M along the laser beam, the camera captures the laser light spot irradiated on the measured object surface M by the laser without any change and images the laser light spot on the M position of the CCD. The optical axis of the laser beam and the CCD are vertical to the optical axis of the imaging lens. According to the image shift h and the structural parameter k of the imaging system, the displacement of the measured object surface can be determined.
If the displacement of the light spot on the imaging plane is h1The displacement of the measured object surface is
<math> <mrow> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>u</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <msub> <mi>h</mi> <mn>1</mn> </msub> </mrow> <mi>f</mi> </mfrac> <mi>cos</mi> <mi>&alpha;</mi> <mo>=</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>h</mi> <mn>1</mn> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow></math>
Where u is the object distance of the imaging system and f is the focal length of the imaging lens. In the above formula, only h is measured1Then d is obtained1
For FIG. 2(b), the laser beam from the laser is at an angle of incidence θ1The air obliquely enters the point A on the air-sea water interface, and the reflection and the refraction occur at the air-water interface. The camera takes an image of the reflected light spot and images it at point D on the CCD. The refracted light is obliquely incident on the lower surface of the oil film, an elliptical laser spot is also formed at an incident point B on the lower surface of the oil film, and the camera captures an image of the spot and images the image at a point C on the CCD. If the displacement of the point C and the point D on the imaging surface is
Figure A200910068627D00052
The included angle between the CCD and the optical axis of the imaging lens is 90 degrees, then
<math> <mrow> <msub> <mi>l</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <msqrt> <msubsup> <mi>n</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>n</mi> <mn>0</mn> <mn>2</mn> </msubsup> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> </msqrt> </mrow> <mrow> <msub> <mi>n</mi> <mn>0</mn> </msub> <mi>f</mi> </mrow> </mfrac> <msubsup> <mi>h</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>=</mo> <msubsup> <mi>k</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <msubsup> <mi>h</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> </mrow></math>
In the formula n0Is the refractive index of air, n1Is the refractive index of seawater, a is the object distance of the imaging system, and f is the focal length of the imaging lens.
When the incident light spot moves from the point A to the point M, the imaging point correspondingly moves on the CCD. From the analysis of FIG. 2(a), the displacement of the object surface is shown
<math> <mrow> <msub> <mi>l</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> </mrow> <mi>f</mi> </mfrac> <mi>cos</mi> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> <mo>=</mo> <msubsup> <mi>k</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> </mrow></math>
Wherein
Figure A200910068627D00055
The image motion is the size of the image motion.
The displacement of the lower surface of the oil film relative to point M is then
<math> <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>l</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>l</mi> <mn>2</mn> </msub> <mo>=</mo> <msubsup> <mi>k</mi> <mn>1</mn> <mtext>&prime;</mtext> </msubsup> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>-</mo> <msubsup> <mi>k</mi> <mn>2</mn> <mtext>&prime;</mtext> </msubsup> <msubsup> <mi>h</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <mi>cos</mi> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> </mrow> <mi>f</mi> </mfrac> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>-</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <msqrt> <msubsup> <mi>n</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>n</mi> <mn>0</mn> <mn>2</mn> </msubsup> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> </msqrt> </mrow> <mrow> <msub> <mi>n</mi> <mn>0</mn> </msub> <mi>f</mi> </mrow> </mfrac> <msubsup> <mi>h</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow></math>
In the formula (2), only need to measureAnd
Figure A200910068627D00058
then d can be obtained2
The invention has the advantages and beneficial effects that:
the invention provides a non-contact measuring method, which has the characteristics of simple principle, convenient measurement, low cost, small volume, capability of eliminating errors introduced by an imaging system, liquid level change, oil film surface gloss and color and external environment, capability of effectively overcoming the influence of background light on measurement precision and the like, and can be used for dynamic online measurement of the thickness of an oil film on the surface of an aqueous solution.
The measuring method can quickly, accurately and quantitatively realize dynamic on-line automatic measurement of the thickness of the oil film of the spilled oil on the sea surface, provides a reliable measuring tool for ocean spilled oil detection, promotes the development of the existing spilled oil monitoring level in China, and provides important basis and technical support for spilled oil treatment.
[ description of drawings ]:
FIG. 1 is a schematic diagram of a symmetrical oblique-type dual-laser triangulation distance-measuring optical path system according to the present invention;
FIG. 2 is a schematic diagram of laser triangulation measurement according to the present invention, in which FIG. 2(a) is a first path of oblique single optical path system, and FIG. 2(b) is a second path of oblique single optical path system;
fig. 3 is an experimental light path diagram of a symmetrical oblique-incidence type double-laser triangulation measurement system.
Fig. 4 is an example of a light spot image obtained by a symmetric oblique-type dual-laser triangulation measurement system, fig. 4(1) is an image acquired by a first path system, and fig. 4(2) is an image acquired by a second path system.
Wherein, 1 is a laser, 2 is a convergent lens, 3 is an imaging lens, 4 is a CCD, 5 is a virtual reference plane, 6 is a data acquisition and processing module, 7 is a sample cell, 8 is a spectroscope, and 18, 19, 28 and 29 are all reflectors.
[ embodiments ] of the present invention:
example 1
The specific measurement process of the measurement method of the present invention is described below by taking the experimental optical path shown in fig. 3 as an example;
optical path system construction
A symmetrical oblique-incidence type double-laser triangular distance measuring system shown in figure 3 is built, and a light source part of the system comprises a laser 1, a converging lens 2, a spectroscope 3 and reflecting mirrors 18, 19, 28 and 29. The laser 1 is a He — Ne laser having a wavelength λ of 632.8nm, and the focusing lens 2 is a focusing lens having a focal length f of 300 mm. The receiving (imaging) system is composed of an imaging lens 3, a detector CCD 4 and a data acquisition and processing module 6. The imaging lens 3 is an imaging lens with the f being 250mm, the detector CCD 4 is an area array CCD, and the data acquisition and processing module 6 comprises an image acquisition card and a computer. The sample cell 7 is a transparent glass cell without a cover, water is filled in the cell, and oil floats on the water surface. In the optical path system, the laser beam, the optical axis of the imaging lens and the CCD are ensured to be in the same plane, the optical axis of the laser beam is strictly ensured to be vertical to the optical axis of the imaging lens, and the CCD is strictly ensured to be vertical to the optical axis of the imaging lens. Adjusting the light path to enable the first path of light beam and the second path of light beam to be intersected in the air, and enabling the upper and lower illuminating points to correspond to the upper and lower surfaces of any point of the object to be measured in the vertical direction;
② light spot image acquisition
Laser beams emitted by the He-Ne laser 1 are reflected by reflectors 18, 19, 28 and 29 after passing through a spectroscope 3, and are converged at a point S on the upper surface of an oil film to be measured by air obliquely entering through a first path of converging lens 2, so that scattering and refraction are generated at an air-oil interface. And adjusting the position of the converging lens 2 to enable the focal point to be positioned near the surface of the oil film to be detected, so that the light spot has strong brightness and a small radius to form a bright spot. The scattered light at point S forms an elliptical laser spot which is imaged by imaging lens 3 onto CCD 4. The reflected light is converged by the converging lens 2 of the second path, enters water through the glass (the bottom surface of the container) by air obliquely incident, is reflected and refracted at an interface point A of the air and the glass, is reflected and refracted at an interface point C of the glass and the water, reflected light at the point A, C forms an elliptical laser spot respectively, and is imaged on the CCD 4 by the imaging lens 3. The refracted light obliquely enters a point B on the lower surface of the oil film and is scattered at an oil-water interface. And adjusting the position of the converging lens 2 to enable the focus to be positioned near the lower surface of the oil film to be detected, forming an elliptical laser spot at an incident point B on the lower surface of the oil film, and imaging the spot on the CCD by an imaging system. Light spot image points of the two-path system on the CCD are collected, as shown in FIG. 4. The oil used in the experiment is petroleum, the image collected by the first path is shown in figure 4(1), and the image collected by the second path is shown in figure 4 (2). In the figure, a shows an image point corresponding to an air glass interface, b shows an image point corresponding to an oil-water interface, C shows an image point corresponding to a glass-water interface (not shown in fig. 1 and 2, if the refractive index of the selected sample cell material is close to that of water/seawater, the point a coincides with the point C), s shows an image point corresponding to an oil-gas interface, and a reference plane is selected as an M point in an oil layer, is a virtual reference plane, and is obtained by calibrating an experimental system. Meanwhile, in the measuring process, the relative spatial positions of the two distance measuring systems are fixed to be unchanged;
③ measuring image displacement
And carrying out image processing on the acquired light spot image by adopting filtering technologies such as morphological filtering and the like, and obtaining the centroid coordinate of the image point on the CCD by utilizing a centroid method. Fig. 4(1) acquires coordinates of the centroid of the light spot at the point S as (1896.9, 1186.2), and coordinates of the centroid of the light spot at the previously calibrated reference plane point M as (1581.1, 1127.6). The coordinates of the centroid of the spot at point B in fig. 4(2) (941.1, 1302.1), at point a (1669.3, 1352.1), at point C (1511.1, 1321.2), and at point M (205.0, 1159.3) are collected. The image displacement quantities h of the upper and lower light spot image points relative to the reference point can be obtained by utilizing a distance formula between the two points12.3768mm and <math> <mrow> <msub> <mi>h</mi> <mn>2</mn> </msub> <mo>=</mo> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>-</mo> <msubsup> <mi>h</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mn>10.9293</mn> <mi>mm</mi> <mo>-</mo> <mn>1.1928</mn> <mi>mm</mi> </mrow></math> <math> <mrow> <mo>-</mo> <mn>4.2204</mn> <mi>mm</mi> <mo>=</mo> <mn>5.5162</mn> <mi>mm</mi> <mrow> <mo>(</mo> <msubsup> <mi>h</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>=</mo> <msubsup> <mi>h</mi> <mi>bc</mi> <mo>&prime;</mo> </msubsup> <mo>+</mo> <msubsup> <mi>h</mi> <mi>ac</mi> <mo>&prime;</mo> </msubsup> <mo>)</mo> </mrow> <mo>.</mo> </mrow></math>
data processing
According to the corresponding displacement h of the light spot on the CCD obtained in the previous step, the data processing module calculates by using the formula d ═ kh, and the displacement d of the upper surface and the lower surface of the oil film relative to the point M can be obtained1And d2In the formula d1=(u-f)h1 cosα/f=k1h1 <math> <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>l</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>l</mi> <mn>2</mn> </msub> <mo>=</mo> <msubsup> <mi>k</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>-</mo> <msub> <mi>k</mi> <mi>bc</mi> </msub> <msubsup> <mi>h</mi> <mi>bc</mi> <mo>&prime;</mo> </msubsup> <mo>-</mo> <msub> <mi>k</mi> <mi>ac</mi> </msub> <msubsup> <mi>h</mi> <mi>ac</mi> <mo>&prime;</mo> </msubsup> <mo>,</mo> </mrow></math> Wherein, <math> <mrow> <msub> <mi>l</mi> <mn>1</mn> </msub> <mo>=</mo> <msubsup> <mi>k</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <msubsup> <mi>h</mi> <mn>1</mn> <mtext>&prime;</mtext> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <mi>cos</mi> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>/</mo> <mi>f</mi> <mo>,</mo> </mrow></math> <math> <mrow> <msub> <mi>l</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>bc</mi> </msub> <msubsup> <mi>h</mi> <mi>bc</mi> <mo>&prime;</mo> </msubsup> <mo>+</mo> <msub> <mi>k</mi> <mi>ac</mi> </msub> <msubsup> <mi>h</mi> <mi>ac</mi> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <msqrt> <msubsup> <mi>n</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>n</mi> <mn>0</mn> <mn>2</mn> </msubsup> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> </msqrt> <msubsup> <mi>h</mi> <mi>bc</mi> <mo>&prime;</mo> </msubsup> <mo>/</mo> <msub> <mi>n</mi> <mn>0</mn> </msub> <mi>f</mi> <mo>-</mo> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <msqrt> <msubsup> <mi>n</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>n</mi> <mn>0</mn> <mn>2</mn> </msubsup> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> </msqrt> <msubsup> <mi>h</mi> <mi>ac</mi> <mo>&prime;</mo> </msubsup> <mo>/</mo> <msub> <mi>n</mi> <mn>0</mn> </msub> <mi>f</mi> <mo>.</mo> </mrow></math> the object distance a, u, 508.38mm, a, 44.24 °, θ of the imaging system is obtained by calibration144.45 degrees, f 250mm, air refractive index n01, refractive index of water n11.33, glass refractive index n3Calculated as 1.5, d in fig. 41=1.7599mm,d21.4963mm, and 3.26mm oil film thickness d.

Claims (2)

1. A sea surface oil spill oil film thickness measuring method based on a differential laser triangulation method is characterized by comprising the following steps:
the method comprises the following steps that firstly, symmetrical oblique-type double-laser triangular distance measuring systems are selected, each laser triangular distance measuring system comprises a light source part consisting of a laser and a converging lens, and a receiving imaging system consisting of an imaging lens, a detector CCD or a position detector PSD (position sensitive detector) and a data acquisition and processing module;
secondly, designing an oil film thickness measuring optical path system, enabling a laser beam emitted by a first laser to be converged and obliquely incident on the upper surface of the oil film to be measured through a converging lens, forming an elliptical laser spot at an incident point on the upper surface of the oil film, and imaging the spot on a CCD (charge coupled device) through a receiving imaging system; the laser beam emitted by the second laser is converged by the convergent lens and obliquely enters the seawater through the air, reflection/scattering and refraction occur at an air-water interface, the refracted light obliquely enters the lower surface of the oil film, an elliptical laser spot is also formed at an incident point on the lower surface of the oil film, the imaging system images the spot on the CCD, and the upper incident point and the lower incident point correspond to the upper surface and the lower surface of any point of the oil film to be measured in the vertical direction; the reference plane is set as a virtual reference plane, and the selection of the reference plane can be arbitrary theoretically;
thirdly, determining the image displacement h of the upper and lower two light spot image points relative to the reference plane light spot image point according to the light spot image points of the two systems on the CCD1And <math> <mrow> <msub> <mi>h</mi> <mn>2</mn> </msub> <mo>=</mo> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>-</mo> <msubsup> <mi>h</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>;</mo> </mrow></math>
fourthly, according to the corresponding displacement h of the light spot on the CCD obtained in the previous step1And h2The data acquisition processing module calculates to obtain the out-of-plane displacement d of the assumed reference plane of the distance between the upper surface and the lower surface of the point1And d2(ii) a Wherein for the first optical path system d1=k1h1In the formula, k1The angle between the optical axis of the laser beam and the normal line of the measured surface is (u-f) cos alpha/f, alpha is the included angle, u is the object distance of the imaging system, and f is the focal length of the imaging lens; for the second optical path system <math> <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>l</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>l</mi> <mn>2</mn> </msub> <mo>=</mo> <msubsup> <mi>k</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <msubsup> <mi>h</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>-</mo> <msubsup> <mi>k</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <msubsup> <mi>h</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>,</mo> </mrow></math> In the formula <math> <mrow> <msubsup> <mi>k</mi> <mn>1</mn> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <msub> <mrow> <mi>cos</mi> <mi>&theta;</mi> </mrow> <mn>1</mn> </msub> <mo>/</mo> <mi>f</mi> <mo>,</mo> </mrow></math> <math> <mrow> <msubsup> <mi>k</mi> <mn>2</mn> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>f</mi> <mo>)</mo> </mrow> <msqrt> <msubsup> <mi>n</mi> <mn>1</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>n</mi> <mn>0</mn> <mn>2</mn> </msubsup> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mn>1</mn> </msub> </msqrt> <mo>/</mo> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mn>0</mn> </msub> <mi>f</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow></math> a is the object distance of the imaging system, θ1Is the incident angle of the laser beam incident on the air-sea interface, n0Is the refractive index of air, n1Is the refractive index of seawater;
fifthly, two out-of-plane displacements d are obtained according to the steps1And d2Subtracting the oil film thickness d ═ d of the point1-d2
2. The measuring method according to claim 1, wherein for the first optical path system, an included angle between an optical axis of the laser beam and an optical axis of the imaging lens is 90 degrees, and the CCD is perpendicular to the optical axis of the imaging lens; for the second optical path system, the optical axis of the laser beam is vertical to the optical axis of the imaging lens, and the angle formed by the CCD and the imaging optical axis is 90 degrees.
CN200910068627XA 2009-04-27 2009-04-27 Method for measuring thickness of sea surface spilled oil film based on differential laser triangulation method Expired - Fee Related CN101532825B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200910068627XA CN101532825B (en) 2009-04-27 2009-04-27 Method for measuring thickness of sea surface spilled oil film based on differential laser triangulation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200910068627XA CN101532825B (en) 2009-04-27 2009-04-27 Method for measuring thickness of sea surface spilled oil film based on differential laser triangulation method

Publications (2)

Publication Number Publication Date
CN101532825A true CN101532825A (en) 2009-09-16
CN101532825B CN101532825B (en) 2010-08-18

Family

ID=41103569

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910068627XA Expired - Fee Related CN101532825B (en) 2009-04-27 2009-04-27 Method for measuring thickness of sea surface spilled oil film based on differential laser triangulation method

Country Status (1)

Country Link
CN (1) CN101532825B (en)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101788272A (en) * 2010-04-06 2010-07-28 天津大学 Method for precisely measuring inner diameter of multi-direction shaft hole based on laser triangulation method
CN102221337A (en) * 2011-03-23 2011-10-19 大连海事大学 Oil film thickness inversion device of ship borne microwave scatterometer
CN102353684A (en) * 2011-06-23 2012-02-15 南京林业大学 Method for acquiring laser meat image by double-laser triangle method
CN102419198A (en) * 2011-09-04 2012-04-18 长春理工大学 High-precision real-time laser triangulation method and triangulation device for liquid level
CN102607434A (en) * 2012-03-08 2012-07-25 天津大学 Oil slick thickness measuring device
CN103217120A (en) * 2013-03-22 2013-07-24 浙江理工大学 Laser thickness-measuring method and device
CN103323466A (en) * 2012-03-20 2013-09-25 苏州中导光电设备有限公司 High-accuracy detection system for line trace of solar silicon wafer
CN103727887A (en) * 2013-12-16 2014-04-16 大连理工大学 Incoherent imaging glass thickness measuring method
CN104748690A (en) * 2013-12-25 2015-07-01 远景能源(江苏)有限公司 Raman spectrum method for measuring the thickness of non-crystal solid electrolyte interface film
CN105526874A (en) * 2015-12-03 2016-04-27 重庆三峡学院 Oil film thickness identification method based on spectral characteristic parameter
CN105716522A (en) * 2014-12-22 2016-06-29 财团法人工业技术研究院 Differential triangulation system and method thereof
CN109084693A (en) * 2018-07-20 2018-12-25 东莞富亚电子有限公司 Color image surveys film thickness method
CN110044280A (en) * 2019-04-23 2019-07-23 华中科技大学 A kind of side focal line method laser triangulation calibrator and method
CN110160452A (en) * 2019-06-19 2019-08-23 哈尔滨工业大学(威海) The measurement method of water-surface oil film thickness based on LR laser raman and laser fluorescence
CN110440700A (en) * 2018-05-02 2019-11-12 长沙青波光电科技有限公司 Target object core intersection detection device and skin oil layer thickness detection device
CN110470231A (en) * 2019-08-07 2019-11-19 上海交通大学 A kind of transparent substance thickness laser measurement method and system
CN110646807A (en) * 2019-09-19 2020-01-03 上海兰宝传感科技股份有限公司 Processing method for resisting ambient light interference based on laser triangulation system
CN111220588A (en) * 2020-03-24 2020-06-02 哈尔滨工业大学(威海) Flow field radiation convergence and dispersion measurement method based on oil film fluorescence brightness
CN113029012A (en) * 2021-04-09 2021-06-25 海南热带海洋学院 Combined type ocean oil spill oil film thickness measurement structure
CN113091665A (en) * 2021-05-10 2021-07-09 中国石油大学(北京) Device and method for measuring oil film thickness
CN114199525A (en) * 2021-12-10 2022-03-18 南京迈得特光学有限公司 Integrated lens measuring device and measuring method thereof
CN114526680A (en) * 2022-01-27 2022-05-24 太原理工大学 Thin ice thickness measuring device and method based on reflected light spot image recognition
CN117460125A (en) * 2023-10-26 2024-01-26 北华航天工业学院 Intelligent control system and method for UV light energy compensation

Cited By (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101788272A (en) * 2010-04-06 2010-07-28 天津大学 Method for precisely measuring inner diameter of multi-direction shaft hole based on laser triangulation method
CN102221337A (en) * 2011-03-23 2011-10-19 大连海事大学 Oil film thickness inversion device of ship borne microwave scatterometer
CN102221337B (en) * 2011-03-23 2012-07-04 大连海事大学 Oil film thickness inversion device of ship borne microwave scatterometer
CN102353684A (en) * 2011-06-23 2012-02-15 南京林业大学 Method for acquiring laser meat image by double-laser triangle method
CN102419198A (en) * 2011-09-04 2012-04-18 长春理工大学 High-precision real-time laser triangulation method and triangulation device for liquid level
CN102607434A (en) * 2012-03-08 2012-07-25 天津大学 Oil slick thickness measuring device
CN103323466B (en) * 2012-03-20 2016-04-13 苏州中导光电设备有限公司 Solar silicon wafers stria high-precision detection system
CN103323466A (en) * 2012-03-20 2013-09-25 苏州中导光电设备有限公司 High-accuracy detection system for line trace of solar silicon wafer
CN103217120B (en) * 2013-03-22 2016-03-09 浙江理工大学 A kind of thickness measurement with laser method and device
CN103217120A (en) * 2013-03-22 2013-07-24 浙江理工大学 Laser thickness-measuring method and device
CN103727887A (en) * 2013-12-16 2014-04-16 大连理工大学 Incoherent imaging glass thickness measuring method
CN103727887B (en) * 2013-12-16 2016-07-06 大连理工大学 Incoherent imaging glass thickness measuring method
CN104748690A (en) * 2013-12-25 2015-07-01 远景能源(江苏)有限公司 Raman spectrum method for measuring the thickness of non-crystal solid electrolyte interface film
CN104748690B (en) * 2013-12-25 2018-01-12 远景能源(江苏)有限公司 Determine the Raman optical spectrum method of the thickness of Amorphous solids electrolyte interface film
CN105716522A (en) * 2014-12-22 2016-06-29 财团法人工业技术研究院 Differential triangulation system and method thereof
CN105716522B (en) * 2014-12-22 2018-04-27 财团法人工业技术研究院 Differential triangulation system and method thereof
CN105526874A (en) * 2015-12-03 2016-04-27 重庆三峡学院 Oil film thickness identification method based on spectral characteristic parameter
CN105526874B (en) * 2015-12-03 2019-09-20 重庆三峡学院 A kind of oil film thickness recognition methods based on spectral signature parameter
CN110440700A (en) * 2018-05-02 2019-11-12 长沙青波光电科技有限公司 Target object core intersection detection device and skin oil layer thickness detection device
CN109084693A (en) * 2018-07-20 2018-12-25 东莞富亚电子有限公司 Color image surveys film thickness method
CN109084693B (en) * 2018-07-20 2020-03-31 东莞富亚电子有限公司 Color image film thickness measuring method
CN110044280A (en) * 2019-04-23 2019-07-23 华中科技大学 A kind of side focal line method laser triangulation calibrator and method
CN110160452A (en) * 2019-06-19 2019-08-23 哈尔滨工业大学(威海) The measurement method of water-surface oil film thickness based on LR laser raman and laser fluorescence
CN110160452B (en) * 2019-06-19 2020-09-01 哈尔滨工业大学(威海) Method for measuring thickness of water surface oil film based on laser Raman and laser fluorescence
CN110470231A (en) * 2019-08-07 2019-11-19 上海交通大学 A kind of transparent substance thickness laser measurement method and system
CN110646807A (en) * 2019-09-19 2020-01-03 上海兰宝传感科技股份有限公司 Processing method for resisting ambient light interference based on laser triangulation system
CN111220588B (en) * 2020-03-24 2023-05-16 哈尔滨工业大学(威海) Flow field spoke aggregation irradiance measurement method based on oil film fluorescence brightness
CN111220588A (en) * 2020-03-24 2020-06-02 哈尔滨工业大学(威海) Flow field radiation convergence and dispersion measurement method based on oil film fluorescence brightness
CN113029012A (en) * 2021-04-09 2021-06-25 海南热带海洋学院 Combined type ocean oil spill oil film thickness measurement structure
CN113091665A (en) * 2021-05-10 2021-07-09 中国石油大学(北京) Device and method for measuring oil film thickness
CN113091665B (en) * 2021-05-10 2022-05-10 中国石油大学(北京) Device and method for measuring oil film thickness
CN114199525A (en) * 2021-12-10 2022-03-18 南京迈得特光学有限公司 Integrated lens measuring device and measuring method thereof
CN114526680A (en) * 2022-01-27 2022-05-24 太原理工大学 Thin ice thickness measuring device and method based on reflected light spot image recognition
CN114526680B (en) * 2022-01-27 2023-07-14 太原理工大学 Thin ice thickness measuring device and measuring method based on reflection light spot image recognition
CN117460125A (en) * 2023-10-26 2024-01-26 北华航天工业学院 Intelligent control system and method for UV light energy compensation
CN117460125B (en) * 2023-10-26 2024-07-12 北华航天工业学院 Intelligent control system and method for UV light energy compensation

Also Published As

Publication number Publication date
CN101532825B (en) 2010-08-18

Similar Documents

Publication Publication Date Title
CN101532825A (en) Method for measuring thickness of sea surface spilled oil film based on differential laser triangulation method
CN108280824B (en) Laser shearing speckle interference defect detection system based on image registration and fusion
CN104808208B (en) Measurement system and method for detecting orientation and size of underwater target on basis of laser sound source
CN101387538B (en) Linear array CCD transmission type liquid level measuring device
CN100451678C (en) High spectrum full polarization three dimension imaging integrate detecting system
CN100504287C (en) Surface plasma resonance measurement mechanism and method thereof
US7570372B2 (en) Optical device for measuring the thickness of an at least partially transparent medium
CN102095469B (en) Device and method for measuring liquid level in storage tank by utilizing camera
CN109298410A (en) A kind of marine oil spill detecting laser radar
Qieni et al. A method for measuring the thickness of transparent oil film on water surface using laser trigonometry
CN204556500U (en) A kind of based on sharp light-struck device for measuring refractive index of transparent liquid
CN103884298B (en) Metal Surface Roughness based on guided mode measures system and method
CN101738243A (en) Method for detecting reflective liquid level displacement based on angular displacement tracking
CN102636457B (en) Trace liquid refractive index measuring system and measuring method
CN102840964A (en) Large-caliber long-focus collimator focal point real-time monitoring system
CN115266509B (en) Underwater vertical suspended matter concentration detection method and system based on laser radar
CN105352915A (en) Refractive index two-dimensional distribution dynamic measurement method
CN111735399B (en) Oil film thickness measuring device and method based on laser-induced fluorescence imaging
CN102621071A (en) Differential interference phase detecting method and device of vector beam excitation surface plasma
CN111189512A (en) Device and method for measuring liquid level and refractive index of transparent liquid based on laser ranging
CN102788678B (en) Oceanic turbulence parameter non-contact measuring method and system
CN100501319C (en) Photoelectric oblique angle measuring device by PSD
CN108776087B (en) Phase critical angle scattering measurement method and device for measuring bubble parameters
CN202599830U (en) Measurement system for refractive index of micro-liquid
CN108759689B (en) coating thickness meter based on double-light-path infrared reflection method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20100818

CF01 Termination of patent right due to non-payment of annual fee