Laser Beam Processing Apparatuses and Correspondent Method Using Multi-beam Interference

Abstract

The invention relates to apparatuses and correspondent method of laser beam processing for various materials with strong light absorption and scattering. The apparatuses can be used for medical no incision laser surgery, long distance underwater or atmosphere light communication, less attenuation light energy delivery in optical turbid media, and so on. This invention is a new use of the imaging method using multiple beam interference to create destructive interference in the beam propagation path to reduce the illumination light intensity and so to reduce absorption and scattering of the materials, and to create constructive interference to produce high composite light intensity which forms an inner light layer to illuminate, process the inner object in the materials. The apparatuses are practicable and have great performance. Compared with the traditional apparatuses, for example, the created laser scalpel can treat tissue at depths of more than 5 cm in human body without incision and with high 3D precision of about 1 μm. The effective light energy delivery distance is more than 1000 m in clear seawater.

Claims

1. A method of laser beam processing for light absorption and scattering materials, the method comprising: using a negative dispersion generation device to broaden the full width of half maximum of a short light pulse; then making the broadened light pulse enter the material containing the object for processing; also utilizing positive dispersion of the material containing the object to compress the broadened light pulse in the propagation path, and to create a short light pulse again which forms an inner light layer to illuminate the object in the material; if needed, making the short signal light pulse reflected from the object return along the incident path reversely; during the return path, the full width of half maximum of the short signal light pulse is broadened by positive dispersion of the material again; then, the broadened signal light pulse is compressed by the negative dispersion generation device; and the broadened signal light pulse becomes short signal light pulse again and is received by imaging receiver placed at the observing position for image-guided processing; and at last increasing the power of the illumination short light pulse to a required value, the object can be processed.

2. The method of laser beam processing for light absorption and scattering materials of claim 1, wherein the frequency range of said short light pulse is in visible region, or/and in infrared, or/and in ultraviolet, or/and in X-ray region(s).

3. The method of laser beam processing for light absorption and scattering materials of claim 1, wherein said image receiving position is designed by making the absolute value of the dispersion generated by the negative dispersion generation device be equal to the absolute value of the dispersion generated by the material within the path reaching the object, but with the opposite sign, and making two optical path distances of the output port of the said short light pulse source and the image receiving position to the object in the material be equal.

4. The method of laser beam processing for light absorption and scattering materials of claim 1, wherein said inner light layer is plane, or cylindrical, or spherical layer located in the material, the thickness of the layer is much thinner than the processing or imaging distance in the material, or this layer is focused to be point or line located in the material.

5. The method of laser beam processing for light absorption and scattering materials of claim 1, wherein said short light pulse may be used to process or/and image the object repeatedly in the material.

6. The method of laser beam processing for light absorption and scattering materials of claim 1, wherein said increasing the power of the illumination short light pulse to a required value is to make the peak intensity of the formed illumination short light pulse in the material be high enough for processing the object, and to make the light intensity of the illumination light pulse in the propagation path be low enough for not damaging the material by light pulse broadening.

7. A apparatus of laser beam processing for light absorption and scattering materials designed based on the method of claim 1, the apparatus comprising: laser generating short light pulse which contains N polarized light beams with different frequencies, the same or approximately the same polarization states, and zero initial phases at a certain moment; the angular frequency intervals Δω of these N beams are equal or not equal but are equal usually; an optical adjusting means to make the amplitudes of the N polarized light beams become the same or approximately the same; a mirrored negative dispersion generation device; an processing or imaging distance adjuster to adjust the processing or imaging distance in the material containing the object; a means to move the formed inner light layer, or light point, or light line in the material in three dimensions.

8. The apparatus of laser beam processing for light absorption and scattering materials of claim 7, wherein said light absorption and scattering materials include human body, animal body, seawater, river water, lake water, pond water, fog, smog, snow, ice, cloud, atmosphere, and any gaseous, liquid or solid materials which have light absorption or/and scattering, especially have strong light absorption or/and scattering.

9. The apparatus of laser beam processing for light absorption and scattering materials of claim 7, wherein said laser beam processing includes medical laser beam treatments, medical laser beam surgery, light communications in atmosphere or water, various light energy delivery in bulk gas, bulk liquid and bulk solid materials for heating, denaturing, ablating, etching, welding, drilling, vaporizing, hitting, cutting, destroying, and so on.

10. The apparatus of laser beam processing for light absorption and scattering materials of claim 7, wherein said laser is mode-locked laser.

11. The apparatus of laser beam processing for light absorption and scattering materials of claim 7, in the wherein said apparatus, the optical elements including the prisms, triangular components, lenses, beam splitters, and so on, all are made of the same material as the material for processing or/and imaging, or all are made of the material which has the same or very approximate same dispersion property as that of the material for processing or/and imaging.

12. The apparatus of laser beam processing for light absorption and scattering materials of claim 7, wherein the number N of said N polarized light beams is from 3 to 10.sup.12 or more, the angular frequency intervals Δω of any two frequency adjacent beams of these N beams are equal or not equal but are equal usually.

13. The apparatus of laser beam processing for light absorption and scattering materials of claim 7, wherein said N polarized light beams are plane polarized, or elliptically polarized, or circularly polarized light beams, wherein said polarization states include polarization directions of the plane polarized light beams, ellipticities of the elliptically polarized light beams.

14. The apparatus of laser beam processing for light absorption and scattering materials of claim 7, wherein said N light beams are plane, or cylindrical, or spherical light beams.

15. The apparatus of laser beam processing for light absorption and scattering materials of claim 7, wherein said mirrored negative dispersion generation device consists chiefly of the prisms and lenses, the output surface of the prism generating the negative dispersion is shaped by computer-controlled high precision grounding and polishing to satisfy the requirements of optical path difference compensations for all pairs of two frequency adjacent beams of the said N beams, and the retroreflective micro-mirror layer is used to make the output surface of the prism generating the negative dispersion become retroreflective surface for reflecting the said N beams with different incident angels reversely.

16. The apparatus of laser beam processing for light absorption and scattering materials of claim 7, wherein said processing or/and imaging distance adjuster consists of two triangular components, which move in opposite directions to adjust the processing or/and imaging distance in the material by changing an additional distance outside the material and offset the extra dispersions caused by component triangular shapes.

17. The apparatus of laser beam processing for light absorption and scattering materials of claim 7, wherein said means to move the said inner light layer, or said point, or said line in the material in three dimensions is by moving reflective mirror(s) to make N polarized light beams scan in two dimensional plane, and adjusting the imaging and processing distance in third dimension.

18. The apparatus of laser beam processing for light absorption and scattering materials of claim 7, wherein said an optical adjusting means to make the amplitudes of N polarized light beams become the same or approximately the same is using dye to make dispersion compensation to laser cavity gain.

19. A apparatus of acoustic processing with or without image guide for sound wave absorption or/and scattering materials designed based on the method of laser beam processing for light absorption and scattering materials of claim 1, the apparatus comprising: a sound wave generator generating N sound waves with different frequencies and the same or not same frequency intervals, the same or approximately the same amplitudes, and the zero initial phases at a certain moment; a mirrored negative dispersion generation device for sound waves; an processing or/and imaging distance adjuster to adjust the expected processing or/and imaging distance in the material, a means to move the processing or/and imaging area in the material in three dimensions.

20. The apparatus of acoustic processing with or without image guide for sound wave absorption or/and scattering materials of claim 19, wherein said mirrored negative dispersion generation device for sound waves generates the acoustic path difference compensations for all pairs of two frequency adjacent sound waves of the said N sound waves for the acoustic path differences produced in the material contains the object for all pairs of two frequency adjacent sound waves of the said N sound waves.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0031] The preferred embodiments of the invented apparatuses will be described underneath. Obviously, these embodiments are not the all apparatuses which can be designed based on the method of this invention. Basing on the method of this invention and using existing technical knowledge, the said apparatus embodiments may be modified and alternated. Therefore, the applicant of this invention reserves the rights of all modifications, alternatives, and equivalent arrangements of the invented apparatus embodiments described underneath.

[0032] The aforementioned aspects and advantages of the invention will be appreciated from the following descriptions of preferred embodiments and accompanying drawings wherein:

[0033] FIG. 1 illustrates a short light pulse is broadened by positive dispersion medium only (in the top), and is broadened by the negative dispersion generation device first and then is compressed by the positive dispersion medium next (in the bottom).

[0034] FIG. 2 illustrates the schematic diagram of optical structure of the medical image-guided no incision laser surgery apparatus.

[0035] FIG. 3 illustrates the schematic diagram of optical structure of the underwater image-guided laser energy delivery apparatus.

[0036] FIG. 4 illustrates the change of composite light intensity of multi-beam interference in dispersive medium. The light intensity I changes with the parameter K. Not that just taking N=20.

DETAILED DESCRIPTION OF INVENTED APPARATUSES AND CORRESPONDENT METHOD

[0037] The invented apparatuses of laser beam processing for light absorption and scattering materials have been designed based on the method described above. The apparatus comprising: laser generating short light pulse which contains N polarized light beams with different frequencies, the same or approximately the same polarization states, and zero initial phases at a certain moment; the angular frequency intervals Δω of these N beams are equal or not equal but are equal usually; an optical adjusting means to make the amplitudes of the N polarized light beams become the same or approximately the same; a mirrored negative dispersion generation device; an processing or imaging distance adjuster to adjust the processing or imaging distance in the material containing the object; a means to move the formed inner light layer, or light point, or light line in the material in three dimensions.

[0038] The short light pulse comes from the mode-locked laser. The same or approximately the same amplitudes of the N polarized light beams is obtained by using dye to make dispersion compensation to laser cavity gain.

[0039] Then, the short pulse enters the input surface of a dispersive medium at t=0. Since the pulse contains N frequency components (the number N is from 3 to 10.sup.12 or more), that is, N light beams, the phase difference of any pair of two beams corresponding to the angular frequencies ω.sub.j and ω.sub.j-1 (j=0,1,2,3, . . . , N−1) in these N beams is zero when the pulse enters the input surface of the medium. Also supposing the input surface of the medium is located at the position of x=0, thus, the initial phases φ.sub.j are zero at x=0 and t=0. The input surface is perpendicular to the x direction.

[0040] In the most situations, the optical media are positive dispersive media including human body and seawater. When the light pulse enters the positive dispersive medium, the different beams constituting the light pulse travel at different speeds. The higher the beam frequency is, the lower the beam travels. Thus, the pulse broadens and become a beam group with weaker and weaker composite light intensity since destructive interference of the multiple beams as shown in the top of FIG. 1.

[0041] In FIG. 1, the light beams are from a mode-locked laser 2. The solid line, dashed line and dotted line represent three planar wavefronts of the beams corresponding to the angular frequencies ω.sub.0, ω.sub.j and ω.sub.N-1. If these three wavefronts travel in the positive dispersive turbid medium 4, since ω.sub.0<ω.sub.j<ω.sub.N-1, the wavefront represented by the solid line travels fastest.

[0042] When these three wavefronts travel in a negative dispersion generation device 8 (see the bottom of FIG. 1), their traveling speeds are reversed, that is, the higher the beam frequency is, the faster the beam travels. The light beams are from mode-locked laser 6. If the light pulse enters the device 8 at t=0, three wavefronts are overlapped at the position x=0. In the device 8, since ω.sub.0<ω.sub.j<ω.sub.N-1, the wavefront presented by the dotted line travels fastest.

[0043] Thus, the phase difference Δϕ.sub.j=ϕ.sub.j−ϕ.sub.j-1 will change with x from zero to negative value although ω.sub.j>ω.sub.j-1. The destructive interference will occur and grow in the device 8 with change of Δϕ.sub.j from zero to negative value.

[0044] Note that the shorter the light pulse duration is, the faster the pulse broadens, and the quicker the pulse peak light intensity decreases. It is because the shorter pulse has wider frequency range and contains more frequency components. If defining the decrease time length of the pulse peak intensity from 100% to a significantly small percentage, such as 0.1%, as the initial broadening period T.sub.ib, outside the initial broadening distance D.sub.ib=V.sub.aT.sub.ib, the light absorption and scattering will become significantly small because the light peak intensity has dropped significantly, where V.sub.a is the average speed of the N beams in the dispersive medium. The required initial broadening period T.sub.ib or the initial broadening distance D.sub.ib depends on the absorption and scattering coefficients of the medium. For the turbid medium with larger absorption or/and scattering, the required T.sub.ib or D.sub.ib should be shorter. For example, the D.sub.ib value should be of the scale of millimeters for medical processing and imaging and of the scale of the meters for underwater processing and imaging. In the same way, the last shortening period T.sub.ls is defined, which is the increase time length of the pulse peak intensity from a very small percentage, such 0.1%, of its maximum value to the 100% of its maximum value. Because the pulse shortening is the reverse process of the pulse broadening completely, T.sub.ib should be equal to T.sub.ls in the same optical medium.

[0045] After leaving the negative dispersion device 8, the broadened light pulse enters the positive dispersive turbid medium 9. The optical path difference of any pair of two frequency adjacent beams of the pulse decreases gradually in the turbid medium 9. The optical path difference between any two of the three wavefronts shown in FIG. 1 also decreases gradually. Thus, in other words, the light pulse broadening terminates, and the light pulse shortening begins.

[0046] Let the phase difference between two beams corresponding to angular frequencies ω.sub.j and ω.sub.j-1 in the turbid medium be Δϕ′.sub.j, the propagation distance of the beam corresponding to angular frequency ω.sub.j be x′.sub.j in the turbid medium, and the refractive indexes of the turbid medium corresponding to the angular frequencies ω.sub.j and ω.sub.j-1 be n′.sub.j and n′.sub.j-1. Then, if the medium used in the device for generating the negative dispersion is the same medium as the turbid medium, or it has the same or very approximate dispersion property as the turbid medium, then, n′.sub.j=n.sub.j and n′.sub.j-1=n.sub.j-1. Under this condition, if taking x′.sub.j=x.sub.j, because Δϕ.sub.j is produced by negative dispersion generation device, we get Δϕ′.sub.j=−Δϕ.sub.j. Thus, because Δϕ′.sub.j−Δϕ.sub.j=0 for every pair of two frequency adjacent beams, the broadened light pulse will be compressed completely. A thin inner light layer will be formed in the turbid medium 9. Therefore, by making mirrored negative dispersion compensation, the expected thin inner light layer can be formed in the turbid medium (please see more detailed descriptions in the attached appendix).

[0047] The said light absorption and scattering materials include human body, animal body, seawater, river water, lake water, pond water, fog, smog, snow, ice, cloud, atmosphere and any gaseous, liquid or solid materials which have light absorption or/and scattering, especially have strong light absorption or/and scattering.

[0048] The said laser beam processing includes medical laser beam treatments, medical laser beam surgery, light communications in atmosphere or water, various light energy delivery in bulk gas, bulk liquid and bulk solid materials for heating, denaturing, ablating, etching, welding, drilling, vaporizing, hitting, cutting, destroying, and so on

[0049] The medical image-guided no incision laser surgery apparatus is described below as the first preferred embodiment of the apparatuses according to the invention. FIG. 2 is the schematic diagram of this apparatus optical structure. Because the sizes of various parts and components differ largely, in order to show necessary details, the shown structure is not drawn in actual proportion.

[0050] The short light pulse comes from a mode-locked fiber laser 100, which is pumped by a light-emitting diode 102. The pump light enters a doped fiber 106 through a coupling element 108. When the total spectral bandwidth of the laser output beams needs to be wide, several light-emitting diodes with different emitting frequencies may be used jointly to pump the fiber 106. An optical isolator 110 and a polarization controller 112 are used to ensure unidirectional beam oscillation. The mode-locked fiber laser with multi-wavelength output has been developed maturely [see the references: N. Li, et al, “Cavity-length optimization for high energy pulse generation in a long cavity passively mode-locked all-fiber ring laser,” Applied Optics, 51, 17, 2012, pp.3726- 3730].

[0051] The laser beams go out through an optical coupler 114, and then the beam diameters are enlarged by lens 116 and 122. After passing through the beam splitter 124, 10% of the light energy enters the second beam splitter 126. The transmittance of the beam splitter 126 is 50%. Then, 5% of the total incident light energy enters two lenses 128 and 130 for beam shrinking and focusing. To use a beam splitter 124 with low transmittance is for less signal light energy loss when the signal light is reflected by beam splitter 124 later. Then, the shrank and slight focused beams enter a right-angle prism 132 normally.

[0052] The vertex angle of the prism 132 is β, therefore, the incident angles of the central lines of the N thin and slight focused beams to the output surface of the prism 132 are β too. When these beams are refracted by the prism 132, the refractive angle θ.sub.j of the central line of the beam corresponding to angular frequencies ω.sub.j is [see the reference: D. S. Goodman, “General principles of Geometric Optics,” in Handbook of Optics, McGRAW-Hill, 1995, Vol. I].

n.sub.j sin β=sin θ.sub.j,   (1)

where n.sub.j is the refractive index of the prism 132 corresponding to the angular frequency ω.sub.j, and the refractive index of atmosphere is approximately 1. After being refracted by the prism 132, all beams enter a thin lens 134. f.sub.3 is the focal distance of lens 134. In order to simplify the analyses and calculations, the central line of the beam corresponding to the angular frequency ω.sub.N-1 is arranged along the optical axis of the lens 134 and through the center of the lens 134. U is the distance from the beam centers on the output surface of the prism 132 to the center of the lens 134. If U<f.sub.3, according to Newton equation for thin lens [see the reference: D. S. Goodman, “General principles of Geometric Optics,” in Handbook of Optics, McGRAW-Hill, 1995, Vol. I], the angle of θ.sub.N-1−θ.sub.j, which is the angle between two central lines of the beams corresponding to the angular frequencies ω.sub.N-1 and ω.sub.j, is magnified by M times to θ′.sub.j, that is

θ′.sub.j=M(θ.sub.N-1=θ.sub.j).   (2)

[0053] Where M=f.sub.3/(U−f.sub.3). When M is negative, the image is a virtual image.

[0054] After the refractive angles are amplified, the N beams enter thin lens 136 all along the directions parallel to the optical axis of the lens 136. f.sub.4 is the focal distance of the lens 136. Then, the N parallel beams all enter prism 138. The input surface of the prism 138 is planar and perpendicular to the optical axis of the lens 136.

[0055] In the prism 138, the height of the central line of the beam corresponding to the angular frequency ω.sub.j is H.sub.j (from the optical axis of the lens 136), the travel distance of the beam corresponding to the angular frequency ω.sub.j is D.sub.j. From FIG. 2 (see more descriptions in the attached appendix),

[00001]$\begin{array}{cc}\frac{H_j}{f_4}=\mathrm{tg}{\theta }_j^{\prime }.& \left(3\right)\end{array}$

[0056] When the beams travel a distance D in the prism 138, two beams corresponding to the angular frequencies ω.sub.j and ω.sub.j-1 will produce an optical path difference ΔP.sub.j as

ΔP.sub.j=DΔn.sub.j, j=0,1,2,3, . . . N−1,   (4)

where

Δn.sub.j=n.sub.j-1−n.sub.j.   (5)

[0057] If the positive dispersion generated in the turbid medium needs to be compensated ideally, all optical path differences produced by all pairs of two frequency adjacent beams in the turbid medium must be generated in the prism 138 equally but with the opposite signs, that is, a mirrored negative dispersion must be generated.

[0058] In order to generate mirrored negative dispersion, the prism 138 must be made of the same medium as the turbid medium, or the prism 138 must have the same or very approximate dispersion property as the turbid medium. To satisfy such a requirement has become relatively easy in recent years. For example, to find a material whose optical property is approximate to the human tissues is not difficult because of the development of tissue simulating phantoms [see the reference: B. W. Pogue, and M. S. Patterson, “Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry,” J. Biomed. Opt., 11(4), 02.1-02.16, 2006]. In the optical spectroscopy, imaging, and therapy research fields, such simulating materials have been widely used. The dispersion, absorption and scattering properties of these materials are characteristic of human tissues. Of course, to choose the material to make the prism 138, its light absorption and scattering should be small for saving light energy. If the chosen material is soft, the prism 138 may be made to be a transparent container and filled with the chosen material.

[0059] Suppose the optical path difference produced by two beams corresponding to the angular frequencies ω.sub.j and ω.sub.j-1 in the turbid medium is ΔP′.sub.j=D′Δn′.sub.j, where Δn′.sub.j=n′.sub.j−n′.sub.j-1, and n′.sub.j and n′.sub.j-1 are the refractive indexes of the turbid medium corresponding to the angular frequencies ω.sub.j and ω.sub.j-1. D′ is the traveling distance of the beam corresponding to the angular frequency ω.sub.j in the turbid medium (D′ may be regarded as the processing or imagine distance in the turbid medium). After satisfying the requirement for medium dispersive property, the negative dispersion generation device needs to generate the following optical path differences

ΔP.sub.j=−ΔP′.sub.j=−D′Δn.sub.j, j=0,1,2,3, . . . , N−1,   (6)

[0060] In addition, except for the prism 138, other optical elements used in the negative dispersion generation device also produce positive dispersions. If the total value of the traveling path lengths of the beams in these elements is much less than D′, these additional positive dispersions can be ignored. Otherwise, they need to be compensated too. Usually, all of the optical elements used in the negative dispersion generation device should be made of the same medium or have the same or very approximate dispersive property as the turbid medium.

[0061] Because the prism 138 can only be made from positive dispersive medium, to produce negative optical path difference for any pair of two beams corresponding to the angular frequencies ω.sub.j and ω.sub.j-1 when ω.sub.j>ω.sub.j-1, the only way is to change the propagation distance difference of any pair of two beams in the prism 138.

[0062] If the propagation distances of the two beams corresponding to the angular frequencies ω.sub.j and ω.sub.j-1 are D.sub.j and D.sub.j-1, respectively, D.sub.j and D.sub.j-1 must satisfy the condition

D.sub.jn.sub.j−D.sub.j-1n.sub.j-1≈ΔD.sub.jn.sub.j=ΔP.sub.j=−ΔP′.sub.j=−D′Δn.sub.j, j=0,1,2,3, . . . , N−1.   (7)

[0063] Where ΔD.sub.j=D.sub.j−D.sub.j-1.

[0064] Rearranging Eq.(7), we have

[00002]$\begin{array}{cc}\Delta D_j=-D^{\prime }\frac{\Delta n_j}{n_j}=D^{\prime }\frac{n_{j-1}-n_j}{n_j},j=0,1,2,3,.Math.,N-1.& \left(8\right)\end{array}$

[0065] Eq. (8) gives the negative optical path difference required for two beams corresponding to the angular frequencies ω.sub.j and ω.sub.j-1 in the prism 138. Thus, the total required propagation distance D.sub.j of the beam corresponding to the angular frequency ω.sub.j in the prism 138 is

D.sub.j=ΔD.sub.N-1+ΔD.sub.N-2+ΔD.sub.N-3+ . . . +ΔD.sub.j.   (9)

[0066] In order to produce mirrored negative dispersion compensations, a method has been crated. It is based on using computer-controlled high precision optical machining and retroreflective micro-mirrors.

[0067] The method is accomplished by measuring the turbid medium refractive indexes corresponding to different frequencies within the required range first. Because the number of the frequencies is large, only partial and discrete data need to be measured. Then, one can use a computer to fit refractive index change curve with frequency from the obtained data. There are several dispersion equations for fitting the refractive index change curves, such as Cauchy, Hartmann, Conrady and Kettler-Drude equations, etc. [see the reference: W. J. Smith, “Optical Materials and Interference Coatings,” in Modern Optical Engineering, McGRAW-Hill, 2000, Chapter 7, p. 176].

[0068] Then ΔD.sub.j can be calculated from the fitted refractive index curve and the required D′ value according to Eq. (8). The D′ is the expected processing or imaging distance in the turbid medium. At last, the computer is used to obtain the total propagation distance D.sub.j for each beam corresponding to the angular frequency ω.sub.j by Eq.(9).

[0069] From Eq.(2) and Eq.(3), we have

H.sub.j=f.sub.4tg[M(θ.sub.N-1=θ.sub.j)].   (10)

[0070] From Eq.(1), we have θ.sub.j=n.sub.j arcsin β and θ.sub.N-1=n.sub.N-1 arcsin β. Thus,

H.sub.j=f.sub.4tg[M(n.sub.N-1 arcsin β−n.sub.j arcsin β)].   (11)

[0071] Using H.sub.j value as the position for a point along the direction perpendicular to the optical axis of the lens 136, and using D.sub.j value corresponding to that H.sub.j as the position for that point along the direction of the optical axis of the lens 136, thus, a data group that contains N data pairs for N points can be produced in the same way. Then, the computer is used to fit a smooth curve which connects the N point positions of (D.sub.j, H.sub.j) from the produced data group. At last, the computer-controlled high precision grinding and polishing are used to shape the output surface of the prism 138 according to the fitted smooth curve.

[0072] In recent years, the computer-controlled high-precision optical grinding and polishing have been developed significantly, which can fabricate the optical elements with very high accuracy [see the reference: D. W. Kim, H. M. Martin, and J. H. Burgea, “Calibration and optimization of computer-controlled optical surfacing for large optics,” Proc. SPIE, 8126, 15.1-15.10, 2011].

[0073] When the N beams are incident on the output surface of the prism 138, because the output surface shape is generally non-planar and non-spherical, the incident angles of the N beams are different. Especially, since the incident beams are focused beams, even the light rays of each beam have different incident angles. Therefore, these beams can't return back along their incident paths by simply making the prism output surface become a reflective surface.

[0074] This problem can be solved by using retroreflective micro-prism mirrors. In recent years, the micro-prisms are used as the tiny tetrahedrons and are placed in arrays on thin hard or soft sheet surface. These retroreflective sheets can reflect light beams within wide spectral range. When the beam incident angle is less than 30°, the reflectivity can be >90%. The average diameter of these micro-prisms is less than 45 μm [see the reference: A. Lundvall, F. Nikolajeff, and T. Lindstrom, “High performing micromachined retroreflector,” Opt. Express, 11(20), 2459-2473, 2003; A. Poscik, J. Szkudlarek, and G. Owczarek, “Photometric properties of retroreflective materials in dependence on their structure and angle of illumination,” Fibres Text. East. Eur. 3(135), 58-64, 2019].

[0075] When a soft micro-prism retroreflective mirror layer is pasted on the smooth output surface of the prism 138 by optical glue, the focused N beams with different incident angles can be returned to travel along their incident paths completely. Note that the returned beams will travel along their previous path once more, so D′ in Eq.(8) should be reduced to 0.5 D′.

[0076] The returned N beams from the prism 138 recombine in the prism 132. After being reflected by beam splitter 126 and mirror 140, these beams become parallel beams and enter an imaging distance adjuster consisting of two triangular components 142 and 146. Because the length and shape of the prism 138 have determined a fixed processing or imaging distance D′ in the turbid medium by the Eq.(8) and Eq.(9), every apparatus has a fixed imaging distance D′. Therefore, if the expected processing or imaging depth in the human body is D.sub.b, the distance changes by adjusting two components 142 and 146 are D.sub.3 and D.sub.4, then to make

D′=D.sub.3+D.sub.4+D.sub.b.   (12)

[0077] The expected processing or imaging depth D.sub.b can be adjusted by changing D.sub.3+D.sub.4. Note that the two components 142 and 146 should be made from the same medium as the turbid medium or have the same or very approximate dispersive property as the turbid medium too. Because two triangular components have symmetrical shapes, no unwanted dispersions will be produced by the distance adjuster.

[0078] Then, the N parallel beams enter the human body 156 by reflection of the mirror 148. If there is no lens 150, an inner light layer will be created at the depth of D.sub.b. The layer thickness is determined by two factors. One is the number N of the beams and the frequency interval Δν(2πΔν=Δω) of the N beams according to the relation of NΔνΔτ=1 [see the reference: W. H. Carters, “Coherence theory,” in Handbook of Optics, McGRAW-Hill, 1995, Vol. I, p.4.3], where Δτ is the duration of the light pulse when Δϕ.sub.j=0, which determines the inner light layer thickness δH by δH=V.sub.hΔτ, where V.sub.h is the average speed of the N beams in the human body. If the total spectral width of the N beams is wide enough, the layer thickness can be very thin, such as less than 1 μm. Another factor is the initial broadening period T.sub.ib. The pulse broadening due to the chromatic dispersion can be estimated as [see the reference: C.-A. Bunge, M. Beckers, and B. Lustermann, Polymer Optical Fibres, Fibre Types, Materials, Fabrication, Characterization and Applications, Elsevier Ltd, Woodhead Publishing, 2017, pp.47-118]

ΔT′=L′Δλd.sub.c,   (13)

where Δλ is the pulse spectral width in wavelength, d.sub.c is chromatic dispersion coefficient, and L′ is the propagation distance of the pulse in the dispersive medium, ΔT′ is full width of half maximum (FWHM) of the pulse. For example, the typical value of d.sub.c is 20 ps/nm.Math.km at 1550 nm for telecom fibers. Thus, if Δλ=1000 nm, which corresponds to 2 fs light pulse, when L′=1 mm, ΔT′=20 fs. Because 2 fs pulse broadens to 20 fs, the peak light intensity of the pulse should drop to below 10% of its maximum value. For seawater, typical d.sub.c values are from 60 ps/nm.Math.km to 300 ps/nm.Math.km [see the reference: “Seawater intrusion and mixing in estuaries,” Marine Species Introduced Traits Wiki, 2020, marinespecies.org/introduced/wiki/Seawater_intrusion_and_mixing_in_estuaries#Experimental_determination_of_the_longitudinal_dispersion_coefficient]. Considering that about 60% of human body is water by weight, thus roughly speaking, the ultra-short light pulse of fs level can broaden fast enough in the human body too (the dispersion coefficients of the human tissues have not been found temporarily). Therefore, if a pulse of fs level broadens by negative dispersion first, then it will shorten fast enough in the human body during the last shortening period T.sub.ls. Thus, the light energy loss due to light absorption and scattering during the last shortening period T.sub.ls is small. Fortunately, obtaining ultrafast, high power fs lasers is not difficult nowadays.

[0079] In FIG. 2, the N beams are focused by lens 150 when they enter the human body 156. It is for confocal imaging to improve the imaging longitudinal resolution, which will be explained later.

[0080] The signal light beams reflected by the targeted tissue return along the incident path reversely. Generally the N beams constituting the incident pulse will all be reflected by the target tissue. The reflections occur on the interface on the targeted tissue surface and between the two areas with different refractive indexes. The reflectivities of the interface for N beams don't differ much usually. During the return path, the signal light pulse will broaden by positive dispersive tissues again as the signal light pulse still contains multiple frequency components, that is, the multiple beams, which results in decrease of the composite intensity of signal light beams, and so results in decrease of the light absorption and scattering in the body again. Then, the signal light beams exit the body. The optical path differences of the signal light beams are further enlarged by positive dispersive imaging distance adjuster. After reflected by beam splitter 126, the signal light beams enter the prism 138 again. This time, the broadened signal pulse will be compressed by negative dispersion by the prism 138. Because the return process of the signal pulse is a completely reverse process of the laser pulse illumination process, the detailed mathematic description does not needed.

[0081] When the signal beams reach the beam splitter 124 again, the signal light beams travel a distance which equals D′ exactly. Thus, the expected signal pulse appears by constructive interference of N signal beams. After reflected by beam splitter 124, as beam splitter 124 has high reflectivity, most of the energy of the signal light pulse is focused on the image plane 166 by lens 162.

[0082] If the lens 150 is not used, the designed apparatus has the most popular imaging structure, which can make one point on the object plane become one point on the image plane directly. This structure can easily combine existing ultra-resolution technologies [see the reference: G. Huszka, and M. A. M. Gijs, “Super-resolution optical imaging: A comparison,” Micro and Nano Eng. 2, 7-28, 2019], such as to place a phase filter before the focusing lens 162. In this way, the imaging resolution along the object plane can exceed the theoretical diffraction limit, which is significantly less than the beam wavelengths.

[0083] Using the lens 150 is for improving the longitudinal resolution. Because during the last shortening period T.sub.ls in the human body, the light pulse intensity will be significantly large. For example, as described above, within the range of 1 mm, the intensity of a 2 fs pulse is about 10% of its maximum value in the telecom fiber. As the initial broadening period and the last shortening period has equal length in the same medium, the effective thickness of a 2 fs pulse will be much larger than its theoretical thickness of approximate 0.2 μm in the human body, which will reduce the longitudinal resolution of the imaging. The confocal imaging can solve this problem [see the reference: S. Inoue, and R. Oldenbourg, “Microscopes,” in Handbook of Optics, McGRAW-Hill, 1995, Vol. II]. By using lens 150 to focus the illumination light beams to scan the targeted tissue, and using a spatial pinhole 168 placed before the image plane 166 to block out-of-focus light in image formation, the imaging longitudinal resolution can be increased to wavelength level, that is, about 1 μm, and with better contrast.

[0084] When using the lens 150, the focal point of the lens 150 is at the position with the depth of D.sub.F in the human body. It should be indicated here, the formed inner light layer and the focal point 160 locate at different positions, which gives a special benefit of more conveniently changing the processing area size and controlling the processing power density, because the area size, and so the power density of the formed inner light layer can be also changed by adjusting the distance difference of D.sub.F−D.sub.b.

[0085] This apparatus creates a thin inner light layer in the human body 156 in the approximate 2D XY plane. The layer area size may also be changed by the distance difference of D.sub.F−D.sub.b. In contrast to the focal point scanning illumination, which is used by many existing 3D imaging or medical treating technologies, this 2D illumination simplifies the apparatus optical structure and improves the imaging and treating speed.

[0086] The mirror 148 can move in the X direction. The apparatus or the human body 156 can move in the Y direction. Thus, by adjusting the processing or/and imaging depth D.sub.b, the 3D imaging or treating can be achieved in the human body 156. The change of the distance difference of D.sub.F−D.sub.b is by moving the lens 150 in the Z direction.

[0087] Many medical treatments need imaging the targeted tissue for guidance. This apparatus has excellent imaging ability. The signal light reflected by the target tissue returns along the incident path reversely. The returned signal light is almost entirely consisted of the previous illumination N beams with reduced amplitudes and somewhat changed polarizations. In the returning path, because the returned beams have different frequencies and different phase velocities, thus, the destructive interference makes the composite light intensity small.

[0088] When the targeted tissue is observed clearly by observer, it also means that the targeted tissue is aimed by the laser beam exactly. Then by raising the output power of the mode-locked laser to the required level for treating the targeted tissue, the targeted tissue can be vaporized, or ablated, or incised by formed inner light layer, which may be called as laser scalpel. When the illumination light power is raised for treating the tissue, the light intensity in the beam propagation path can be still under the safe threshold value since multiple beam interference (see further descriptions below). Afterwards, by reducing the output power of the mode-locked laser to previous level, the result of such no incision laser surgery can be checked by imaging the targeted tissue again.

[0089] The method of said invention may further combine a variety of existing technologies to produce a variety of laser processing with or without image guidance. As these are existing technologies and knowledge, no further explanation is needed too.

[0090] The underwater image-guided laser energy delivery apparatus is described below as the second preferred embodiment of the apparatuses according to the invention. FIG. 3 is the schematic diagram of this apparatus optical structure. Because the sizes of various parts and components differ largely, in order to show necessary details, the shown structure is not drawn in actual proportion too.

[0091] In FIG. 3, the mode-locked fiber laser 200 consists of light-emitting diode 202, doped fiber 204, coupling element 206, optical isolator 208, polarization controller 210 and optical coupler 212. The laser beam diameters are enlarged by lens 214 and 216. After passing through beam splitters 218 and 220, the laser beams enter negative dispersion generation device consisting of lenses 222, 226, 230 and 236, prisms 228 and 238. In the prism 238, the negative dispersion is produced as that in the above described medical laser surgery apparatus.

[0092] The processing or/and imaging distance adjuster is composed of two paralleled mirrors 246 and 248, and two triangular components 250 and 252. Two mirror planes are inclined at an angle θ.sub.M to Z axis.

[0093] Since the desired processing or/and imaging distance underwater is long, after N parallel beams entering the distance adjuster, each beam will be reflected multiple times in the adjuster. Making the diameter of each beam be small, thus each beam can obtain a larger number of reflections between two mirrors. If the expected processing or/and imaging distance in the seawater 266 is D.sub.b , the distance change by adjusting the components 250 and 252 is ΠD.sub.w, and the designed processing or/and imaging distance of the apparatus is D′. By making D′=D.sub.b+ΠD.sub.w, the D.sub.b can be adjusted by changing ΠD.sub.w. Here, Π is the number of the reflection times of a beam between two mirrors, D.sub.w is the travelling distance of a beam in two components between two reflections. Because the diameters of the beams are small, such as a diameter of 5 mm, the required thickness of the imaging distance adjuster is thin, such as less than 10 mm, so the distance adjuster can have moderate volume and light weight. Furthermore, if required, the distance adjuster can be a composite distance adjuster composed of multiple distance adjusters.

[0094] The change of the value of ΠD.sub.w is by moving the components 250 and 252 simultaneously along the mirror planes in the opposite directions. Moving the components 250 and 252 simultaneously and in the opposite directions is for offset extra dispersions caused by the triangular shapes of two components.

[0095] After the diameters of the N beams are expanded by lens 254 and 256, the N parallel beams enter the seawater 266 and form an inner light layer on the object plane 260 at the position with the distance of D.sub.b.

[0096] This apparatus forms a 2D illumination, which will simplify the imaging optical structure and improve imaging speed, since a 2D inner light layer is formed in the YZ plane in the seawater 266. The layer area is determined by the cross-sectional area of the N beams group. The Group of the N beams can scan up and down, from the right to the left, in order to achieve large range 3D illumination and energy delivery in the seawater 266.

[0097] Compared with the medical imaging and treating apparatus, the depth resolution requirement for the underwater imaging and delivery is generally much lower than that for the medical imaging and treating. In the seawater, the expected imaging and delivery distance is from several meters to even kilometers, thus the depth resolution of 0.1 mm to 1 mm is very enough generally, which is 2 to 3 orders of magnitude lower than the requirement for medical apparatus. Therefore, the total spectral bandwidth of the N beams is 2 to 3 orders of magnitude narrower than that of the medical apparatus too.

[0098] The signal light produced by reflection from the object in the formed light layer returns back along the incident path reversely, and going through a process similar to that of the medical imaging described above. At last, the signal light beams are reflected by the beam splitter 218 and create the expected signal light pulse, which is focused on the image plane 270 by lens 268.

[0099] In the same way, in order to improve the longitudinal resolution of imaging, a spatial pinhole 276 is placed before the image plane 270 to form the confocal imaging.

[0100] The lens 268 and the 2D image plane 270 also form the most common camera structure, which makes one object point become one image point, and so it is easy to get high imaging resolution and fast imaging speed.

[0101] The returned signal light rays from the points of the object plane are not drawn in FIG. 3. In addition, the parallel illumination light beams may be focused on the object plane, then by point scanning way to produce the image on the image plane 270 or to focus the light energy to a tiny point on the object plane.

[0102] Also when the object in the seawater is observed clearly by observer, it means that the light energy can be delivered to that object by raising or not raising the output power of the mode-locked laser. For example, the underwater light communication may not need to raise the laser beam power. Some other applications, such as to hit shark or other dangerous creatures needs to raise the laser beam power. If the laser beam power is raised to the level high enough, more objects including non-life ones can be processed or destroyed too. Afterwards, by reducing the output power of the mode-locked laser to the imaging level, the result of light energy delivery can be checked by imaging the object.

[0103] Similarly, based on the method of the invention, a variety of existing technologies can combined to create a variety of new functions for underwater apparatus. Since these works may be done by using the existing knowledge, no more descriptions are given here.

Brief Description of Performances of the Invented Apparatuses

[0104] In the following description, the imaging and laser beam processing performances and signal enhancement by using the said apparatuses are given.

[0105] Following the way of the analyses and derivations described in the existing invention titled” Method of inner light layer illumination by multi-beam interference and apparatuses for imaging in turbid media” (see the attached appendix), we have the composite light intensity of the multi-beam interference of the N beams in dispersive medium is

[00003]$\begin{array}{cc}I=A^2\left(x\right)\frac{{\cos }^2\left[\frac{\left(N-1\right)\left(K\Delta \omega \right)}{2}\right]{\sin }^2\left(\frac{\mathrm{NK}\Delta \omega }{2}\right)}{{\sin }^2\left(\frac{K\Delta \omega }{2}\right)}.& \left(14\right)\end{array}$

where K=t−(2xn.sub.0/C), A is the amplitude of the N beams (supposing the amplitudes of the N beams have the same value), t is the time, x is the beam traveling distance in the turbid medium, n.sub.0 is the refractive index of the turbid medium corresponding to the angular frequency ω.sub.0, C is light speed in vacuum.

[0106] When KΔω becomes zero, the value of I goes to the maximum. The results of numerical calculations by Eq.(14) are shown in FIG. 4 and Table 1. In FIG. 4, I changes with parameter K. 300 and 302 are two composite light intensity maximums. Using K as the unit of transverse coordinate is for avoiding complicated theoretical derivations, and the essential characteristics of the multi-beam interference in optical dispersive medium can still be shown. The I curve in FIG. 4 may be regarded as the change of the composite light intensity of N beams with x at the moment when the inner light layer is recreated. In the calculations, A=1. In order to show the width of the composite light intensity maximum obviously, just taking N=20.

[0107] In the Table 1, N is the number of the beams participating in the interference. γ is the enhancement factor of the composite light intensity maximum γI.sub.0. ε is the attenuation factor of the remaining composite light intensity εI.sub.0 between two composite light intensity maximums (see FIG. 4). I.sub.0 is the incident intensity of each beam of the N beams. We see that γ increases with N increase rapidly (see attached appendix for the detailed descriptions about the calculations).

TABLE-US-00001 TABLE 1 Signal Intensity Enhancement at Different Imaging Depth in Human Body by Multi-beam Interference. N γ ε Factor D′ = 2 cm D′ = 5 cm D′ = 10 cm D′ = 15 cm D′ = 20 cm 1.00E+01 1.00E+02 5.05E−14 ξ 1.10E−01 1.10E−01 1.10E−01 1.10E−01 1.10E−01 1.00E+02 1.00E+04 6.11E−10 ξ 1.10E+01 1.10E+01 1.10E+01 1.10E+01 1.10E+01 1.00E+03 9.97E+05 6.21E−06 ξ 1.10E+03 1.10E+03 1.09E+03 1.09E+03 1.09E+03 1.00E+04 7.08E+07 2.80E−06 ξ 7.79E+04 7.79E+04 7.78E+04 7.78E+04 7.78E+04 1.00E+05 7.08E+09 1.13E−05 ξ 7.78E+06 7.78E+06 7.76E+06 7.75E+06 7.74E+06 1.00E+06 7.08E+11 9.35E−06 ξ 7.79E+08 7.78E+08 7.77E+08 7.76E+08 7.75E+08 1.00E+07 7.08E+13 5.31E−06 ξ 7.79E+10 7.78E+10 7.78E+10 7.77E+10 7.77E+10 1.00E+01 1.00E+02 5.05E−14 α (dB) 2.63E+02 6.27E+02 1.23E+03 1.84E+03 2.45E+03 1.00E+02 1.00E+04 6.11E−10 α (dB) 2.83E+02 6.47E+02 1.25E+03 1.86E+03 2.47E+03 1.00E+03 9.97E+05 6.21E−06 α (dB) 3.03E+02 6.67E+02 1.27E+03 1.88E+03 2.49E+03 1.00E+04 7.08E+07 2.80E−06 α (dB) 3.21E+02 6.85E+02 1.29E+03 1.90E+03 2.51E+03 1.00E+05 7.08E+09 1.13E−05 α (dB) 3.41E+02 7.05E+02 1.31E+03 1.92E+03 2.53E+03 1.00E+06 7.08E+11 9.35E−06 α (dB) 3.61E+02 7.25E+02 1.33E+03 1.94E+03 2.55E+03 1.00E+07 7.08E+13 5.31E−06 α (dB) 3.81E+02 7.45E+02 1.35E+03 1.96E+03 2.57E+03

[0108] The calculation results shown in FIG. 4 and Table 1 indicate that the larger the N value, the higher the light peak intensity γI.sub.0, the narrower the full-width of half maximum of the light peak intensity γI.sub.0. At the same time, the larger the N value, the smaller the remaining light intensity εI.sub.0 between two peak light intensities.

[0109] The numerical calculations show that when the N changes from 10.sup.1 to 10.sup.7, the enhancement factor γ of the peak intensity γI.sub.0 changes from 10.sup.2 to 10.sup.14, and the attenuation factor ε of the remain light intensity εI.sub.0 changes from 10.sup.−14 to 10.sup.−6. The calculation results are shown in Table 1.

[0110] In the Table 1, the difference between the composite peak light intensity γI.sub.0 and the remaining composite light intensity εI.sub.0 may be over 18 orders of magnitude. Such large intensity difference can certainly give plentiful room to avoid light injury for human tissues or light damage for materials in the laser beam propagation path. In the medical surgery applications, the light power density of less than 10 mW/cm2 is safe for human tissues including skin [see reference: T. J. Dougherty, et al, “Photodynamic Therapy (Review),” Journal of the National Cancer Institute, Vol. 90, No. 12, 1998, pp.889-905]. And the power density of higher than 10 W/cm2 can ablate most targeted tissues in human body without problem. The difference between 10 mW/cm2 and 10 w/cm2 is just 3 orders of magnitude. Therefore, no incision laser surgery can be completed by the invented apparatus. For the underwater light energy delivery, hitting a shark perhaps needs a light power density of about 100 W/mm2 because that a light power density of 500 W/mm2 can cut a steel plate [see reference: Miyamoto and H. Maruo, “Mechanism of laser cutting,” Welding in the World, Le Soudage Dans Le Monde, Vol. 29, No. 9/10, 1991, pp.283-294]. Therefore, if the light power density of the laser beam is 50 mW/mm2 in the propagation path, which should not heat the water obviously. The difference between 50 mW/cm2 and 500 W/cm2 is 4 orders of magnitude. Therefore, good underwater light energy delivery can also be completed by the invented apparatus.

[0111] Now, based on the actual absorption and scattering situations of the seawater and human body, the signal intensity change and imaging sensitivity enhancement of the said apparatuses can be calculated. In the human body, the absorption coefficient and scattering coefficient are different for different tissues. Here, the average absorption coefficient μ.sub.ba=0.397 mm.sup.−1 and scattering coefficient μ.sub.bs=1 mm.sup.−1 of the human blood are taken for the whole body temporarily [see above Dr. M. C. Hillman's doctoral thesis]. Although taking the coefficients of blood for whole human body is differ from the actual situation, as mentioned above, since there is a large amount of blood in the human body and the coefficients of the blood have approximate order of magnitude as those of the most human tissues, such treatments can give approximate results. In addition, there is a practical method for determining the expected processing or imaging distance in the media consisting of the compositions with different refractive indices. it will be given at the last of this invention. The light reflectance R is assumed for the case that light is reflected from the interface between the blood and adipose. The refractive indices of the blood and adipose are taken as 1.3771 and 1.4714, respectively (correspondent to the wavelength of 680 nm). Then, according to the Fresnel formula, at the boundary of two media with different refractive indices of n.sub.1 and n.sub.2, the amplitude reflectance r is

[00004]$\begin{array}{cc}r=\frac{n_1-n_2}{n_1+n_2},& \left(15\right)\end{array}$

[see the reference: J. M. Bennett, “Polarization,” Chapter 5, Handbook of Optics, Vol. I, 2ed, Edited by M. Bass, and et al, McGRAW-Hill, New York, 1995, p.5.7], and the light intensity reflectance R=r.sup.2. We get the light intensity reflectance R=0.0011 for this interface, Also according to the intensity enhancement factor γ and the intensity attenuation factor ε shown in the Table 1, we get the signal intensity change factor ξ and the imaging sensitivity enhancement factor α values with different processing or imaging distances of 2 cm, 5 cm, 10 cm, 15 cm and 20 cm in the human body as the follows in Table 1.

[0112] The absorption coefficient and scattering coefficient of the seawater vary a lot according to the different situations, here taking the clear seawater as the example. The absorption coefficient μ.sub.wa of clear seawater is 0.0196 m.sup.−1, and its scattering coefficient μ.sub.ws is 0.0212 m.sup.−1 [see above reference written by C. D. Mobley]. The refractive index of water at the wavelength of 550 nm is 1.336. Assuming the refractive index of the object is 1.6 (note that the optical glass refractive index range is 1.5 to 2.0), then an assumed light intensity reflectance of R=0.00809 for object in clear seawater is obtained. According to the intensity enhancement factor γ of the peak light intensity and the intensity attenuation factor ε of the remain light intensity in the Table 2, we get the signal intensity change factor ξ and the imaging sensitivity enhancement factor α values with different processing or imaging distances of 200 m, 500 m, 1000 m, 1500 m and 2000 m in the clear seawater shown in Table 2 (also see attached appendix for the detailed descriptions about the calculations).

[0113] From Table 1 and 2, we can see that the signal composite light intensity maximum is much higher than the signal intensity of the normal imaging. For example, when D′=5 cm for medical processing or imaging, or D′=500 m for underwater processing or imaging, the intensity enhancement factor is more than 1.1×10.sup.3 when N>10.sup.3. It means that when N>10.sup.3, and the total spectra width of N beams is wide enough, the peak intensity of the signal light pulse can be higher than NI.sub.0.Note that NI.sub.0 is the average value of the total intensity of the N beams (if the N beams are incoherent light beams). Of course, the extreme high pulse peak intensity is always with the extreme narrow pulse duration usually, and so the energy of each pulse may be very low.

[0114] However, as long as the signal to noise ratio is high, the required signal energy can be got by receiving repeated signal pulses. It can be seen that the value ξ is still high even when D′=20 cm for medical processing or imaging, or D′=2000 m for underwater processing or imaging. Therefore, there is good potential to get the signal intensity enhancement factor of near 2000 dB, whose corresponding processing or imaging depth and distance are 20 cm in the human body and 2000 m in the clear seawater. Considering the approximations are made in the calculations, the signal intensity enhancement factor of 600 dB is taken for representing the apparatus performances, whose corresponding processing or imaging depth is 5 cm in the human body, and 1000 m in the clear seawater.

TABLE-US-00002 TABLE 2 Signal Intensity Enhancement at Different Imaging Depth in Seawater by Multi-beam Interference N Γ ε Factor D′ = 200 m D′ = 500 m D′ = 1000 m D′ = 1500 m D′ = 2000 m 1.00E+01 1.00E+02 5.05E−14 ξ 1.10E−01 1.10E−01 1.10E−01 1.10E−01 1.10E−01 1.00E+02 1.00E+04 6.11E−10 ξ 1.10E+01 1.10E+01 1.10E+01 1.10E+01 1.10E+01 1.00E+03 9.97E+05 6.21E−06 ξ 1.10E+03 1.10E+03 1.09E+03 1.09E+03 1.09E+03 1.00E+04 7.08E+07 2.80E−06 ξ 7.79E+04 7.79E+04 7.78E+04 7.78E+04 7.78E+04 1.00E+05 7.08E+09 1.13E−05 ξ 7.78E+06 7.78E+06 7.76E+06 7.75E+06 7.74E+06 1.00E+06 7.08E+11 9.35E−06 ξ 7.79E+08 7.78E+08 7.77E+08 7.76E+08 7.75E+08 1.00E+07 7.08E+13 5.31E−06 ξ 7.79E+10 7.78E+10 7.78E+10 7.77E+10 7.77E+10 1.00E+01 1.00E+02 5.05E−14 α (dB) 2.63E+02 6.27E+02 1.23E+03 1.84E+03 2.45E+03 1.00E+02 1.00E+04 6.11E−10 α (dB) 2.83E+02 6.47E+02 1.25E+03 1.86E+03 2.47E+03 1.00E+03 9.97E+05 6.21E−06 α (dB) 3.03E+02 6.67E+02 1.27E+03 1.88E+03 2.49E+03 1.00E+04 7.08E+07 2.80E−06 α (dB) 3.21E+02 6.85E+02 1.29E+03 1.90E+03 2.51E+03 1.00E+05 7.08E+09 1.13E−05 α (dB) 3.41E+02 7.05E+02 1.31E+03 1.92E+03 2.53E+03 1.00E+06 7.08E+11 9.35E−06 α (dB) 3.61E+02 7.25E+02 1.33E+03 1.94E+03 2.55E+03 1.00E+07 7.08E+13 5.31E−06 α (dB) 3.81E+02 7.45E+02 1.35E+03 1.96E+03 2.57E+03

[0115] As shown in Table 1 and Table 2, the enhancement of the imaging sensitivity of the said apparatuses is very large. The said apparatuses have such great performance is not strange, because it is created by multiple beam interference. In the past, multiple beam interference has demonstrated its astonishing abilities, such as to create extremely short light pulse of dozens of attoseconds (1 attosecond=10.sup.−18 s) and extremely strong light power of several terawatts (1 terawatt=10.sup.12 watts). They are the fastest-ever and strongest-ever man-made events until now. In the future, the multiple beam interference will surely make more technical contributions.

[0116] At last, we give the practical method for adjusting the additional distance to get accurate expected imaging and laser processing distance in the turbid media or human body. In a turbid medium or human body, there are various compositions with different refractive indices, so it is difficult to find the accurate value of the additional distance determined by all refractive indices of the compositions in the turbid medium or human body. However, there is a practical way to overcome this difficulty. Because when the angular phase differences between any frequency adjacent beams become zero, the composite light intensity maximum emerges certainly. Therefore, somewhat like to search for a music station by tuning the frequency of a radio, no matter what the accurate distance of the object or the targeted tissue position is, one just needs observing into the turbid medium or human body and adjusting the distance adjuster at the same time. When the searched object or targeted tissue appears in the visual field (by camera) and becomes clear, the accurate expected processing or imaging distance is achieved.

[0117] The acoustic wave has good ability to permeate dense liquid and solid materials when the wave frequency is not high. It is used widely for detecting the information in dense materials, such as the injury in bulk solid materials. Therefore, to increase acoustic imaging and processing distance in dense materials is very useful.

[0118] Because optical waves and acoustic waves all are oscillation waves, such waves all cosinusoidally vary, that is, these waves all have cosine or sine forms. Therefore, multiple beam (wave) interference can certainly happen to them, and result in the similar physical results, that is, using multiple wave interference to reduce the composite wave intensity in the propagation path, and increase composite wave intensity at expected positions. Therefore, the apparatus of acoustic processing or imaging/detection in dense materials, that is, in turbid media can be designed based on the principle of this invented method.

[0119] The apparatus of acoustic processing with or without image-guild for sound wave absorption or/and scattering materials designed based on the method of laser beam processing for light absorption and scattering materials comprising: a sound wave generator generating N sound waves with different frequencies and the same or not same frequency intervals, the same or approximately the same amplitudes, and the zero initial phases at a certain moment; a mirrored negative dispersion generation device for sound waves; an processing or/and imaging distance adjuster to adjust the expected processing or/and imaging distance in the material, a means to move the processing or/and imaging area in the material in three dimensions.

[0120] The said mirrored negative dispersion generation device for sound waves generates the acoustic path difference compensations for all pairs of two frequency adjacent sound waves of the said N sound waves for the acoustic path differences produced in the material contains the object for all pairs of two frequency adjacent sound waves of the said N sound waves.

[0121] The said sound wave generator may be a mode-locked acoustic laser. The said imaging or detection distance adjustor consists of the components similar to those in the above described medical and underwater imaging apparatuses, but can make acoustic wave traveling distance change. The said means to move imaging or detecting area in the turbid media in three dimensions is using the reflector(s) for acoustic waves. Because the acoustic imaging or detection apparatus can be designed based on the same principle of this invented method and existing knowledge, no more description is given here.

[0122] Similarly, based on the principle of this invented method, a variety of existing technologies can combined to create a variety of new functions. Since these works may be done by using the existing knowledge, no more descriptions are given here.