Three-dimensional dynamic imaging characterization device for micro-nano bubbles
Technical Field
The invention relates to a device and a method for measuring the individual size, the particle size distribution, the three-dimensional position and the three-dimensional motion track of micro-nano bubbles in a solution and on the surface.
Background
Compared with other particles such as solid particles, oil drops and the like, the micro-nano bubbles have the physical characteristics of small size, high possibility of disturbance, variable form, wide particle size distribution, small refractive index compared with water and the like. At present, there are several methods for characterizing micro-nano bubbles, such as atomic force microscopy for surface bubbles; as the solution bubble, there are a laser light scattering measurement method, a multi-wavelength extinction method, a two-dimensional image method, a low-field nuclear magnetic method, and the like. The characterization methods have respective advantages, but still have common limitations, namely the in-situ and real-time change process of the micro-nano bubble individual cannot be accurately and comprehensively observed under the three-dimensional scale, and bubbles and other solid or liquid impurities with similar sizes cannot be distinguished. Therefore, we developed a three-dimensional dynamic imaging technology for solution and surface micro-nano bubbles based on Digital Holographic Microscopy (DHM). The method records real-time holographic images of a plurality of micro-scale particles by utilizing light interference, carries out three-dimensional reconstruction on the images to obtain the motion tracks, the shapes and the phases of the micro-scale particles, and has the advantages of high flux, large depth of field, high precision, no damage, capability of giving material information (optical refractive index) and the like.
In the prior art, a system and a method for observing micro-nano bubbles in a water body disclose a function of imaging micro-nano bubbles with the diameter larger than 900 nanometers and obtaining particle size distribution and a motion track. The method can image the micro-nano bubbles with the diameter as small as 200 nanometers based on holography, and can obtain the particle size distribution and the three-dimensional dynamic track of the micro-nano bubbles. In addition, the device has the characteristic of high observation precision.
Disclosure of Invention
The invention provides a micro-nano bubble three-dimensional dynamic imaging and characterization device and method. The invention mainly aims to make up the general defects of the prior art, realize three-dimensional dynamic observation of micro-nano bubbles and have the advantages of simple and convenient optical path, large flux and high resolution. The method for dynamically observing the micro-nano bubbles in the solution and on the surface based on the three-dimensional holographic microscopic imaging can be used for representing the physicochemical characteristics of the micro-nano bubbles and researching the interaction mechanism of the micro-nano bubbles and the surrounding medium.
The invention is realized by the following technical scheme.
A three-dimensional dynamic imaging characterization device for micro-nano bubbles comprises: the device comprises a light source, a convex lens group, a sample cell, an objective lens and an image sensor; the light source, the convex lens group, the sample cell, the objective lens and the image sensor are sequentially positioned on the same horizontal line, the centers of the light source, the convex lens group, the objective lens and the image sensor are all positioned on a dot-dash line optical axis (dot-dash line) in the figure, and the surface of the sample cell is vertical to the optical axis; the light path is simple.
Light emitted by the light source passes through the convex lens group and is converged into converging parallel light with higher energy density. The incident light enters the sample cell, a part of the incident light is scattered after irradiating the micro-nano bubbles, the formed scattered light is used as object light, and the object light is interfered with the incident light at the defocused position on the image sensor through the objective lens to form a holographic interference pattern and is recorded by the image sensor.
The using method of the device comprises the following steps:
s1, turning on the power supply of the light source (1) and the image sensor (5);
s2, adjusting parameters of the image sensor (5) to make the brightness of the field uniform and appropriate;
s3, adding a solution rich in micro-nano bubbles into the sample cell (3), or after adding the solution, manufacturing the micro-nano bubbles by an injection or stirring method, standing for a proper time, and controlling the image sensor (5) to record a sample hologram in a high-speed continuous shooting manner;
s4, obtaining a background-subtracted hologram from the hologram by using a difference method;
s5, reconstructing a three-dimensional light field by the background-subtracted hologram by using a three-dimensional reconstruction algorithm;
s6, processing the reconstructed light field by using a track tracking algorithm to obtain motion track information of the micro-nano bubbles;
and S7, processing the reconstructed light field by using a particle size measurement algorithm to obtain the size information of each micro-nano bubble.
Further, the light source is coherent light or partially coherent light; the central wavelength of the light source is selected within the range of 200nm < lambda <1000 nm.
Further, the centers of the light source, the convex lens group, the objective lens and the image sensor are on the same vertical line. The distances of the units of the light source, the convex lens group, the objective lens and the image sensor in the optical axis direction satisfy the following imaging positions: and the focusing position is positioned on the inner wall of the sample pool close to the objective lens side, so that the micro-nano bubble imaging is properly defocused.
Furthermore, the convex lens group is replaced by a beam expander or a spatial filter is added in the light path, so that incident light spots can be further homogenized.
Further, the objective lens may be an air-medium objective lens or a liquid-medium objective lens, and when the liquid-medium objective lens is used, the front end of the objective lens may be directly contacted with the liquid in the sample cell by a special sample cell.
Further, the optical path length can be shortened by adding a prism group between the objective lens and the image sensor.
Further, the apparatus may add focusing and field of view shifting functionality by incorporating a shifting device in the optical path.
Further, the specific process of step S4 is: the hologram is averaged in light intensity to obtain a background image, which is calculated as follows:
wherein, Ib(x, y) is the gray scale value of the pixel at the (x, y) position in the background image, N is the total frame number of the hologram photographed at high speed, I is the specific frame number, I is the number of framesiAnd (x, y) is the gray value of the pixel at the (x, y) position in the original holographic image of the ith frame.
Further, in step S4, the specific process is as follows: and subtracting the background image from the hologram to obtain the background-subtracted hologram.
Further, the specific process of step S5 is: carrying out numerical reconstruction on the holographic image with the background removed to obtain the information of the three-dimensional reconstruction intensity distribution of the micro-nano bubbles, and calculating as follows:
U(r,-z)=FT-1(FT(Is(r,0)·H(q,-z)))
the initial coordinate of the micro-nano bubble in the horizontal direction is shown, and z is the initial axial coordinate of the micro-nano bubble; FT-1Is inverse Fourier transform; FT is Fourier transform; h (q, -z) is the Fourier transform of H (r, -z).
Further, in step S5, the value of the reconstruction axial interval is larger than the axial imaging range of the instrument, and the step is the size of the pixel divided by the magnification of the objective lens.
Further, in step S6, the three-dimensional positions of the signal points of the micro-nano bubbles are connected in series according to the correlation of the adjacent frames in the three-dimensional light field, so as to obtain a three-dimensional motion trajectory.
Further, in step S7, fitting a light intensity curve of the micro-nano bubbles in the light field with a function to obtain a half-peak width of the signal points;
and obtaining the size of each micro-nano bubble according to the corrected proportional coefficient of the half-peak width and the size of the micro-nano bubble and the half-peak width of each micro-nano bubble.
The image sensor collects interference light of object light and incident light of the micro-nano bubbles scattered by the incident light near a focal plane as a hologram.
And processing the hologram by adopting a difference method to obtain the background-subtracted hologram.
And reconstructing the interference image with the background subtracted into a three-dimensional light field by using a three-dimensional reconstruction algorithm comprising Fourier transform and inverse Fourier transform.
And processing the three-dimensional light field by using a three-dimensional track tracking algorithm to obtain a three-dimensional track of the micro-nano bubbles.
And processing the three-dimensional light field by using a particle size measurement algorithm to obtain the individual size and the particle size distribution of the micro-nano bubbles.
Further, the Gaussian fitting position of the micro-nano bubble obtained by the particle size measurement algorithm can be used for correction and verification of the trajectory tracking algorithm.
In conclusion, by utilizing the hologram of the micro-nano bubbles, the individual size distribution, the three-dimensional density distribution and the three-dimensional motion track of the micro-nano bubbles can be obtained, and the three-dimensional observation of a large number of micro-nano bubbles with micron and submicron sizes by an image method can be realized.
Compared with the prior art, the invention has the advantages that: the micro-nano bubbles with the diameter as small as 200 nanometers in the liquid surface interface and the bulk phase can be imaged by holography, and the particle size distribution and the three-dimensional dynamic track of the micro-nano bubbles can be obtained. And the device has the characteristics of simplicity, convenience, feasibility and high observation precision in the design aspect.
Drawings
FIG. 1 is a schematic diagram of a three-dimensional digital holographic imaging characterization device of micro-nano bubbles;
fig. 2 is an original hologram of the micro-nano bubble collected by the three-dimensional digital holographic imaging characterization device for the micro-nano bubble in embodiment 1 of the present invention, where the number represents the time sequence, and the hologram of a part of the bubble is framed;
fig. 3 is a hologram of the micro-nano bubbles with background subtraction obtained by the three-dimensional digital holographic imaging characterization device for the micro-nano bubbles in embodiment 1 through data processing, wherein the number represents a time sequence, and diffraction rings of part of the bubbles are framed;
fig. 4 is a three-dimensional trajectory diagram of the micro-nano bubbles collected by the three-dimensional digital holographic imaging characterization device of the micro-nano bubbles in embodiment 1 of the invention, wherein the size of a sphere is the size of the micro-nano bubbles;
fig. 5 is a size distribution of micro-nano bubbles collected by the three-dimensional digital holographic imaging characterization device for micro-nano bubbles in embodiment 1 of the present invention.
Fig. 6 is a hologram of a background subtraction of micro-nano bubbles and a three-dimensional trajectory diagram of a certain nano bubble in the hologram in embodiment 2 of the present invention, where the size of a micro-sphere in the figure is the size of the micro-nano bubble;
fig. 7 is a hologram of a background subtraction of micro-nano bubbles and a three-dimensional trajectory diagram of a certain nano bubble in the hologram in embodiment 3 of the present invention, where the size of a micro-sphere in the figure is the size of the micro-nano bubble;
fig. 8 is a hologram of a background subtraction of micro-nano bubbles and a three-dimensional trajectory diagram of a certain nano bubble in the hologram in embodiment 4 of the present invention, where the size of a micro-sphere in the figure is the size of the micro-nano bubble;
fig. 9 is a size distribution of the micro-nano bubbles acquired by 6 times of recording in the three-dimensional digital holographic imaging characterization device of the micro-nano bubbles in embodiment 4 of the invention.
Detailed Description
The following is a detailed description of a micro-nano bubble three-dimensional imaging characterization device and an image acquisition mode according to the present invention with reference to the accompanying drawings.
A three-dimensional dynamic imaging characterization device for micro-nano bubbles comprises: the device comprises a light source 1, a convex lens group 2, a sample cell 3, an objective lens 4 and an image sensor 5; the light source 1, the convex lens group 2, the sample cell 3, the objective lens 4 and the image sensor 5 are sequentially positioned on the same horizontal line, the centers of the light source (1), the convex lens group 2, the objective lens 4 and the image sensor 5 are all positioned on the optical axis of the dot-dash line in the figure, and the surface of the sample cell 3 is vertical to the optical axis; the distances of the units of the light source 1, the convex lens group 2, the objective lens 4 and the image sensor 5 in the optical axis direction satisfy the following imaging positions: and the focusing position is positioned on the inner wall of the sample cell close to the objective lens side, so that the micro-nano bubble imaging is properly defocused. The light source is coherent light or partially coherent light. The central wavelength selection range of the light source is as follows: 200nm < lambda <1000 nm.
The using method of the device comprises the following steps:
s1, turning on the power supply of the light source 1 and the image sensor 5;
s2, adjusting parameters of the image sensor 5 to make the brightness of the field uniform and appropriate;
s3, adding a solution rich in micro-nano bubbles into the sample cell 3, or after adding the solution, manufacturing the micro-nano bubbles by an injection or stirring method, standing, and controlling the image sensor 5 to record a sample hologram in a high-speed continuous shooting manner;
s4, obtaining a background-subtracted hologram from the original hologram;
s5, reconstructing a three-dimensional light field by the background-subtracted hologram by using a three-dimensional reconstruction algorithm;
s6, processing the three-dimensional reconstruction light field by using a track tracking algorithm to obtain motion track information of the micro-nano bubbles;
and S7, processing the reconstructed light field by using a particle size measurement algorithm to obtain the size information of each micro-nano bubble.
The specific process of step S4 is: and (3) carrying out light intensity averaging on the hologram to obtain a background image, and subtracting the background image from the hologram to obtain a background-subtracted hologram.
The specific process of step S5 is: carrying out numerical reconstruction on the holographic image with the background removed to obtain the information of the three-dimensional reconstruction intensity distribution of the micro-nano bubbles, and calculating as follows:
U(r,-z)=FT-1(FT(Is(r,0)·H(q,-z)))
the initial coordinate of the micro-nano bubble in the horizontal direction is shown, and z is the initial axial coordinate of the micro-nano bubble; FT-1Is inverse Fourier transform; FT is Fourier transform; h (q, -z) is the Fourier transform of H (r, -z).
In step S5, the value of the reconstruction axial interval is larger than the axial imaging range of the instrument, and the step is the size of the pixel divided by the magnification of the objective lens.
In step S6, connecting the three-dimensional positions of adjacent micro-nano bubble signal points according to the correlation of adjacent frames in the three-dimensional light field, and obtaining a three-dimensional motion trajectory.
In step S7, for the signal points of the micro-nano bubbles in the optical field, the light intensity is fitted by using a function, and the half-peak width of the signal points is obtained.
And obtaining the size of each micro-nano bubble according to the corrected proportional coefficient of the half-peak width and the size of the micro-nano bubble and the half-peak width of each micro-nano bubble.
Example 1
FIG. 1 is a schematic diagram of a micro-nano bubble three-dimensional imaging characterization device. The device comprises a light source 1, a convex lens group 2, a sample cell 3, an objective lens 4 and an image sensor 5.
The light source 1 uses a partially coherent LED light source with a center wavelength selected in the range of 200nm < λ <1000nm, where a blue LED light source is used. Light rays emitted by the light source pass through the beam expander on the optical axis to obtain parallel light beams.
The parallel light beam becomes the incident light irradiating the sample cell 3 through the convex lens group 2 on the optical axis, the convex lens group 2 is composed of 2 or 3 convex lenses arranged in front and back on the optical axis, the convex lens group 2 needs to collect the light intensity of the homogenized parallel light beam as much as possible. And the light beam is converged and adjusted again to be parallel light, and the obtained converged parallel light needs to have excellent parallelism in the space where the sample cell 3 is located.
The converged parallel light is used as incident light to irradiate the sample cell 3, two surfaces of the sample cell 3 in the optical axis direction need to be perpendicular to the optical axis, and the two surfaces are made of high-light-transmittance quartz glass with the flatness superior to that of the optical surface with the central wavelength 1/4 of the light source 1.
Pure water containing micro-nano bubbles generated by a micro-nano bubble generator is injected into the sample cell 3. Incident light irradiating the micro-nano bubbles is scattered to become object light, namely the propagation direction of the dotted object light in the figure is related to the micro-nano bubbles, the three-dimensional position and size information of the micro-nano bubbles are contained, and the wavelength of the object light is the same as that of the original incident light.
The object light is converged by the objective lens 4 on the optical axis, and a 40-fold objective lens is used in the present embodiment to obtain high resolution by taking both the magnification and the field definition in high-speed shooting.
The object light converged by the objective lens 4 is recorded as an original hologram by the high-speed image sensor 5 on the optical axis, resulting in a hologram set, as shown in fig. 2. The image sensor 5 may be a CMOS, CCD, or the like image sensor.
The computer controls the image acquisition of the image sensor 5 as a control device. The light source 1 supporting digital control can also be controlled.
The computer acts as a data analysis device and processes the acquired set of original holograms using an algorithm in accordance with the principles of holography. Let the gray value of the pixel at the (x, y) position in the ith frame image in the group image of the original hologram be Ii(x, y) obtaining a background map by calculating:
Ib(x, y) is the gray scale value of the pixel at the (x, y) position in the background image, and N is the total number of frames of the hologram photographed at high speed. And subtracting the background image from the hologram group image by using a difference value method to obtain the background-subtracted hologram.
The obtained background-subtracted holographic interference pattern is a plurality of interference rings with alternate light and shade. The holographic image with the background removed is subjected to numerical reconstruction by using a three-dimensional reconstruction algorithm based on Fourier transform, and the calculation is as follows:
U(r,-z)=FT-1(FT(Is(r,0)·H(q,-z)))
r is an initial coordinate of the micro-nano bubble in the horizontal direction for a propagation operator,z is an initial axial coordinate of the micro-nano bubble; i issIs the light intensity of the sample; FT-1Is inverse Fourier transform; FT is Fourier transform; h (q, -z) is the Fourier transform of the light field propagation factor H (r, -z). In the aspect of parameter setting, the value of the reconstruction axial interval is larger than the imaging range of the instrument in the axial direction, and the stepping distance is the same as the pixel size of the image.
And running a three-dimensional reconstruction algorithm to obtain and derive the information of the three-dimensional reconstruction light intensity distribution of each frame of hologram in the micro-nano bubble group diagram.
And (3) searching a local maximum value of light intensity through threshold filtering on the information of the three-dimensional reconstruction intensity distribution of the micro-nano bubbles obtained by reconstructing the light field, and filtering noise points with low reconstruction light intensity to obtain the three-dimensional position of the micro-nano bubble point.
And (3) according to the correlation of the signal points of the micro-nano bubbles of adjacent frames in the background-subtracted hologram group diagram in the three-dimensional light field, executing a track tracking algorithm on the information of the three-dimensional reconstruction light intensity distribution, setting a tracking step length according to the approximate movement speed of the micro-nano bubbles, and connecting the three-dimensional positions of the signal points of the micro-nano bubbles in series to obtain a three-dimensional movement track.
And (4) running a size measurement algorithm on the signal points of the micro-nano bubbles in the optical field. Firstly, a function is used for fitting a light intensity curve of the signal point to obtain the half-peak width of the signal point. And obtaining the size of each micro-nano bubble according to the corrected proportional coefficient of the half-peak width and the size of the micro-nano bubble and the half-peak width of each micro-nano bubble.
And finally obtaining a three-dimensional locus diagram of the micro-nano bubbles containing the size information, such as fig. 4, wherein the size of the micro-sphere is the size of the micro-nano bubbles, and the size distribution diagram of the micro-nano bubbles is shown in fig. 5.
Example 2
The light source 1 in example 1 is a green LED light source, and an ethanol aqueous solution containing micro-nano bubbles generated by an alcohol-water substitution method is injected into the sample cell 3, and the rest is unchanged, so as to obtain a background-reduced micro-nano bubble hologram and a three-dimensional trajectory diagram of one of the bubbles, as shown in fig. 6.
Example 3
The light source 1 in example 1 is an ultraviolet LED light source, ethanol aqueous solution containing micro-nano bubbles generated by an alcohol-water replacement method is injected into the sample cell 3, and the rest is unchanged, so as to obtain a background-reduced micro-nano bubble hologram and a three-dimensional track diagram of one of the bubbles, as shown in fig. 7.
Example 4
The objective lens 4 in example 3 was changed to a 100 x oil medium objective lens to pursue as high a magnification as possible. Then, an image is collected under a 100-fold mirror, a background-subtracted micro-nano bubble hologram and a three-dimensional track diagram of a bubble are obtained as shown in fig. 8, and 649 micro-nano bubble tracks obtained by 6-time total 6-min shooting are statistically distributed in size as shown in fig. 9.
The above embodiments are preferred embodiments of the present invention, but the present invention is not limited to the above embodiments, and any other changes, modifications, substitutions, combinations, and simplifications which do not depart from the spirit and principle of the present invention should be construed as equivalents thereof, and all such changes, modifications, substitutions, combinations, and simplifications are intended to be included in the scope of the present invention.