改善三维波动方程叠前深度偏移处理效果的方法研究

3D wave equation prestack depth migration is the effective tool for obtaining the exact imaging result of complex geology structures. It's a part of the 3D seismic data processing. 3D seismic data processing belongs to high dimension signal processing, and there are some difficult problems to do with. They are: How to process high dimension operators? How to improve the focusing? and how to construct the deconvolution operator? The realization of 3D wave equation prestack depth migration, not only realized the leap from poststack to prestack, but also provided the important means to solve the difficult problems in high dimension signal processing. In this thesis, I do a series research especially for the solve of the difficult problems around the 3D wave equation prestack depth migration and using it as a mean. So this thesis service for the realization of 3D wave equation prestack depth migration for one side and improve the migration effect for another side. This thesis expatiates in five departs. Summarizes the main contents as the follows: In the first part, I have completed the projection from 3D data point area to low dimension are using de big matrix transfer and trace rearrangement, and realized the liner processing of high dimension signal. Firstly, I present the mathematics expression of 3D seismic data and the mean according to physics, present the basic ideal of big matrix transfer and describe the realization of five transfer models for example. Secondly, I present the basic ideal and rules for the rearrange and parallel calculate of 3D traces, and give a example. In the conventional DMO focusing method, I recall the history of DM0 process firstly, give the fundamental of DMO process and derive the equation of DMO process and it's impulse response. I also prove the equivalence between DMO and prestack time migration, from the kinematic character of DMO. And derive the relationship between DMO base on wave equation and prestack time migration. Finally, I give the example of DMO process flow and synthetic data of theoretical models. In the wave equation prestak depth migration, I firstly recall the history of migration from time to depth, from poststack to prestack and from 2D to 3D. And conclude the main migration methods, point out their merit and shortcoming. Finally, I obtain the common image point sets using the decomposed migration program code.In the residual moveout, I firstly describe the Viterbi algorithm based on Markov process and compound decision theory and how to solve the shortest path problem using Viterbi algorithm. And based on this ideal, I realized the residual moveout of post 3D wave equation prestack depth migration. Finally, I give the example of residual moveout of real 3D seismic data. In the migration Green function, I firstly give the concept of migration Green function and the 2D Green function migration equation for the approximate of far field. Secondly, I prove the equivalence of wave equation depth extrapolation algorithms. And then I derive the equation of Green function migration. Finally, I present the response and migration result of Green function for point resource, analyze the effect of migration aperture to prestack migration result. This research is benefit for people to realize clearly the effect of migration aperture to migration result, and study on the Green function deconvolution to improve the focusing effect of migration.

基于信息重建技术的地震信号处理方法研究

The dissertation addressed the problems of signals reconstruction and data restoration in seismic data processing, which takes the representation methods of signal as the main clue, and take the seismic information reconstruction (signals separation and trace interpolation) as the core. On the natural bases signal representation, I present the ICA fundamentals, algorithms and its original applications to nature earth quake signals separation and survey seismic signals separation. On determinative bases signal representation, the paper proposed seismic dada reconstruction least square inversion regularization methods, sparseness constraints, pre-conditioned conjugate gradient methods, and their applications to seismic de-convolution, Radon transformation, et. al. The core contents are about de-alias uneven seismic data reconstruction algorithm and its application to seismic interpolation. Although the dissertation discussed two cases of signal representation, they can be integrated into one frame, because they both deal with the signals or information restoration, the former reconstructing original signals from mixed signals, the later reconstructing whole data from sparse or irregular data. The goal of them is same to provide pre-processing methods and post-processing method for seismic pre-stack depth migration. ICA can separate the original signals from mixed signals by them, or abstract the basic structure from analyzed data. I surveyed the fundamental, algorithms and applications of ICA. Compared with KL transformation, I proposed the independent components transformation concept (ICT). On basis of the ne-entropy measurement of independence, I implemented the FastICA and improved it by covariance matrix. By analyzing the characteristics of the seismic signals, I introduced ICA into seismic signal processing firstly in Geophysical community, and implemented the noise separation from seismic signal. Synthetic and real data examples show the usability of ICA to seismic signal processing and initial effects are achieved. The application of ICA to separation quake conversion wave from multiple in sedimentary area is made, which demonstrates good effects, so more reasonable interpretation of underground un-continuity is got. The results show the perspective of application of ICA to Geophysical signal processing. By virtue of the relationship between ICA and Blind Deconvolution , I surveyed the seismic blind deconvolution, and discussed the perspective of applying ICA to seismic blind deconvolution with two possible solutions. The relationship of PC A, ICA and wavelet transform is claimed. It is proved that reconstruction of wavelet prototype functions is Lie group representation. By the way, over-sampled wavelet transform is proposed to enhance the seismic data resolution, which is validated by numerical examples. The key of pre-stack depth migration is the regularization of pre-stack seismic data. As a main procedure, seismic interpolation and missing data reconstruction are necessary. Firstly, I review the seismic imaging methods in order to argue the critical effect of regularization. By review of the seismic interpolation algorithms, I acclaim that de-alias uneven data reconstruction is still a challenge. The fundamental of seismic reconstruction is discussed firstly. Then sparseness constraint on least square inversion and preconditioned conjugate gradient solver are studied and implemented. Choosing constraint item with Cauchy distribution, I programmed PCG algorithm and implement sparse seismic deconvolution, high resolution Radon Transformation by PCG, which is prepared for seismic data reconstruction. About seismic interpolation, dealias even data interpolation and uneven data reconstruction are very good respectively, however they can not be combined each other. In this paper, a novel Fourier transform based method and a algorithm have been proposed, which could reconstruct both uneven and alias seismic data. I formulated band-limited data reconstruction as minimum norm least squares inversion problem where an adaptive DFT-weighted norm regularization term is used. The inverse problem is solved by pre-conditional conjugate gradient method, which makes the solutions stable and convergent quickly. Based on the assumption that seismic data are consisted of finite linear events, from sampling theorem, alias events can be attenuated via LS weight predicted linearly from low frequency. Three application issues are discussed on even gap trace interpolation, uneven gap filling, high frequency trace reconstruction from low frequency data trace constrained by few high frequency traces. Both synthetic and real data numerical examples show the proposed method is valid, efficient and applicable. The research is valuable to seismic data regularization and cross well seismic. To meet 3D shot profile depth migration request for data, schemes must be taken to make the data even and fitting the velocity dataset. The methods of this paper are used to interpolate and extrapolate the shot gathers instead of simply embedding zero traces. So, the aperture of migration is enlarged and the migration effect is improved. The results show the effectiveness and the practicability.

ACOUSTO-OPTIC IMAGING IN DIFFUSE MEDIA USING PULSED ULTRASOUND AND THE PHOTOREFRACTIVE EFFECT

Acousto-optic imaging (AOI) in optically diffuse media is a hybrid imaging modality in which a focused ultrasound beam is used to locally phase modulate light inside of turbid media. The modulated optical field carries with it information about the optical properties in the region where the light and sound interact. The motivation for the development of AOI systems is to measure optical properties at large depths within biological tissue with high spatial resolution. A photorefractive crystal (PRC) based interferometry system is developed for the detection of phase modulated light in AOI applications. Two-wave mixing in the PRC creates a reference beam that is wavefront matched to the modulated optical field collected from the specimen. The phase modulation is converted to an intensity modulation at the optical detector when these two fields interfere. The interferometer has a high optical etendue, making it well suited for AOI where the scattered light levels are typically low. A theoretical model for the detection of acoustically induced phase modulation in turbid media using PRC based interferometry is detailed. An AOI system, using a single element focused ultrasound transducer to pump the AO interaction and the PRC based detection system, is fabricated and tested on tissue mimicking phantoms. It is found that the system has sufficient sensitivity to detect broadband AO signals generated using pulsed ultrasound, allowing for AOI at low time averaged ultrasound output levels. The spatial resolution of the AO imaging system is studied as a function of the ultrasound pulse parameters. A theoretical model of light propagation in turbid media is used to explore the dependence of the AO response on the experimental geometry, light collection aperture, and target optical properties. Finally, a multimodal imaging system combining pulsed AOI and conventional B- mode ultrasound imaging is developed. B-mode ultrasound and AO images of targets embedded in both highly diffuse phantoms and biological tissue ex vivo are obtained, and millimeter resolution is demonstrated in three dimensions. The AO images are intrinsically co-registered with the B-mode ultrasound images. The results suggest that AOI can be used to supplement conventional B-mode ultrasound imaging with optical information.

Identification of fluorescent beads using a coded aperture snapshot spectral imager.

We apply a coded aperture snapshot spectral imager (CASSI) to fluorescence microscopy. CASSI records a two-dimensional (2D) spectrally filtered projection of a three-dimensional (3D) spectral data cube. We minimize a convex quadratic function with total variation (TV) constraints for data cube estimation from the 2D snapshot. We adapt the TV minimization algorithm for direct fluorescent bead identification from CASSI measurements by combining a priori knowledge of the spectra associated with each bead type. Our proposed method creates a 2D bead identity image. Simulated fluorescence CASSI measurements are used to evaluate the behavior of the algorithm. We also record real CASSI measurements of a ten bead type fluorescence scene and create a 2D bead identity map. A baseline image from filtered-array imaging system verifies CASSI's 2D bead identity map.

Imaging through turbulence using compressive coherence sensing.

Previous studies have shown that the isoplanatic distortion due to turbulence and the image of a remote object may be jointly estimated from the 4D mutual intensity across an aperture. This Letter shows that decompressive inference on a 2D slice of the 4D mutual intensity, as measured by a rotational shear interferometer, is sufficient for estimation of sparse objects imaged through turbulence. The 2D slice is processed using an iterative algorithm that alternates between estimating the sparse objects and estimating the turbulence-induced phase screen. This approach may enable new systems that infer object properties through turbulence without exhaustive sampling of coherence functions.

Acoustic radiation force impulse imaging (ARFI) on an IVUS circular array.

Our long-term goal is the detection and characterization of vulnerable plaque in the coronary arteries of the heart using intravascular ultrasound (IVUS) catheters. Vulnerable plaque, characterized by a thin fibrous cap and a soft, lipid-rich necrotic core is a precursor to heart attack and stroke. Early detection of such plaques may potentially alter the course of treatment of the patient to prevent ischemic events. We have previously described the characterization of carotid plaques using external linear arrays operating at 9 MHz. In addition, we previously modified circular array IVUS catheters by short-circuiting several neighboring elements to produce fixed beamwidths for intravascular hyperthermia applications. In this paper, we modified Volcano Visions 8.2 French, 9 MHz catheters and Volcano Platinum 3.5 French, 20 MHz catheters by short-circuiting portions of the array for acoustic radiation force impulse imaging (ARFI) applications. The catheters had an effective transmit aperture size of 2 mm and 1.5 mm, respectively. The catheters were connected to a Verasonics scanner and driven with pushing pulses of 180 V p-p to acquire ARFI data from a soft gel phantom with a Young's modulus of 2.9 kPa. The dynamic response of the tissue-mimicking material demonstrates a typical ARFI motion of 1 to 2 microns as the gel phantom displaces away and recovers back to its normal position. The hardware modifications applied to our IVUS catheters mimic potential beamforming modifications that could be implemented on IVUS scanners. Our results demonstrate that the generation of radiation force from IVUS catheters and the development of intravascular ARFI may be feasible.

Near-field spatial imaging of a Ni-like Ag 140-angstrom x-ray laser

We measure the two-dimensional, near-field spatial distribution of a 140-Angstrom nickel-like silver x-ray laser at the output aperture with high magnification using a curved multilayer x-ray mirror to image the output onto an x-ray charge-coupled device camera. Lasing is created by illuminating silver slab targets with a pair of 75 ps laser pulses separated by 2.2 nsec from the Vulcan laser. The two-dimensional, high-resolution, spatial image shows the x-ray laser source size and its position relative to the target surface. A dramatic change in both the position and source size are observed for the refraction compensating curved target as compared with the flat targets.

Near Field Enhancement and Subwavelength Imaging Using Resonantly Loaded Apertures

It is demonstrated that the electromagnetic (EM) transmission through a subwavelength or non-resonant aperture in a conductive screen can be dramatically enhanced by loading it with folded metallic strips exhibiting resonant properties. When illuminated by an EM plane wave these loaded apertures enable very tight, subwavelength, collimation of the EM power in the near field zone. We propose planar and quasi-planar resonant insertion geometries that should allow, for the first time, two-dimensional dual-polarization subwavelength field confinement along with ability to focus both electric and magnetic fields. The proposed technique for resonance transmission enhancement and near field confinement forms a basis for a new class of microwave near field imaging probe with subwavelength resolution capable of operating over a wide range of imaging distances (0.05–$0.25lambda$). Measurement results demonstrate the possibility of high contrast (more than 3 dB in amplitude and 40 degrees in phase) near field subwavelength imaging of 2D and 3D resonant and non-resonant metallic and dielectric targets in free space and in moderately lossy layered media.

Resonantly loaded apertures for high-resolution near field surface imaging

A novel type of microwave probes based on the loaded aperture geometry has been proposed and experimentally evaluated for dielectrics characterisation and high-resolution near-field imaging. Experimental results demonstrate the possibility of very accurate microwave spectroscopic characterisation of thin lossy dielectric samples and biological materials containing water. High-resolution images of the subwavelength lossy dielectric strips and wet and dry leaves have been obtained with amplitude contrast around 10-20 dB and spatial resolution better than one-tenth of a wavelength in the near-field zone. A microwave imaging scenario for the early-stage skin cancer identification based on the artificial dielectric model has also been explored. This model study shows that the typical resolution of an artificial malignant tumour with a characteristic size of one-tenth of a wavelength can be discriminated with at least 6 dB amplitude and 50° phase contrast from the artificial healthy skin and with more than 3 dB contrast from a benign lesion of the same size. It has also been demonstrated that the proposed device can efficiently deliver microwave energy to very small, subwavelength, focal areas which is highly sought in the microwave hyperthermia applications.

High-Resolution Microwave Near-Field Imaging Using Resonance Probes

A novel microwave high-resolution near-field imaging technique is proposed and experimentally evaluated in reflectometry imaging scenarios involving planar metal-dielectric structures. Two types of resonance near field probes-a small helix antenna and a loaded subwavelength slot aperture are studied in this paper. These probes enable very tight spatial field localization with the full width at half maximum around one tenth of a wavelength, λ, at λ/100-λ/10 standoff distance. Importantly, the proposed probes permit resonance electromagnetic coupling to dielectric or printed conductive patterns, which leads to the possibility of very high raw image resolution with imaged feature-to-background contrast greater than 10-dB amplitude and 50° phase. In addition, high-resolution characterization of target geometries based on the cross correlation image processing technique is proposed and assessed using experimental data. It is shown that printed elements features with subwavelength size ~λ/15 or smaller can be characterized with at least 10-dB resolution contrast.

MPGDs based radiation imaging devices and applications

Este trabalho descreve o desenvolvimento e aplicação de sistemas baseados em detetores gasosos microestruturados, para imagiologia de fluorescência de raios-X por dispersão em energia (EDXRF). A técnica de imagiologia por fluorescência de raios-X assume-se como uma técnica poderosa, não-destrutiva, em análises da distribuição espacial de elementos em materiais. Os sistemas para imagiologia de EDXRF desenvolvidos são constituídos por: um tubo de raios-X, usado para excitar os elementos da amostra; um detetor gasoso microestruturado; e uma lente pinhole que foca a radiação de fluorescência no plano do detetor formando assim a imagem e permitindo a sua ampliação. Por outro lado é estudada a influência do diâmetro da abertura do pinhole bem como do fator de ampliação obtido para a imagem, na resolução em posição do sistema. Foram usados dois conceitos diferentes de detetores gasosos microestruturados. O primeiro, baseado na microestrutura designada por 2D-Micro-Hole & Strip Plate (2D-MHSP) com uma área ativa de 3 3 cm2, enquanto que o segundo, baseado na estrutura 2D-Thick-COBRA (2D-THCOBRA) apresenta uma área ativa de deteção de 10 10 cm2. Estes detetores de raios-X de baixo custo têm a particularidade de funcionar em regime de fotão único permitindo a determinação da energia e posição de interação de cada fotão que chega ao detetor. Deste modo permitem detetar a energia dos fotões X de fluorescência, bem como obter imagens 2D da distribuição desses fotões X para o intervalo de energias desejado. São por isso adequados a aplicações de imagiologia de EDXRF. Os detetores desenvolvidos mostraram resoluções em energia de 17% e 22% para fotões incidentes com uma energia de 5.9 keV, respectivamente para o detetor 2D-MHSP e 2D-THCOBRA e resoluções em posição adequadas para um vasto número de aplicações. Ao longo deste trabalho é detalhado o desenvolvimento, o estudo das características e do desempenho de cada um dos detetores, e sua influência na performance final de cada sistema proposto. Numa fase mais avançada apresentam-se os resultados correspondentes à aplicação dos dois sistemas a diversas amostras, incluindo algumas do nosso património cultural e também uma amostra biológica.

Non-invasive imaging of roots with high resolution X-ray micro-tomography

ray micro-tomography is a well-established technique for non-invasive imaging and evaluation of heterogeneous materials. An inexpensive X-ray micro-tomography system has been designed and built for the specific purposes of examining root growth and root/soil interactions. The system uses a silver target X-ray source with a focal spot diameter of 80 mum, an X-ray image intensifier with a sampling aperture of about 100 mum, and a sample with a diameter of 25 mm. Pre-germinated wheat and rape seeds were grown for up to 8-10 days in plastic containers in a sandy loam soil sieved to < 250 μm, and imaged with the X-ray system at regular intervals. The quality of 3 D image obtained was good allowing the development and growth of both root axes and some first-order laterals to be observed. The satisfactory discrimination between soil and roots enabled measurements of root diameter (wheat values were 0.48-1.22 mm) in individual tomographic slices and, by tracking from slice to slice, root lengths were also measured. The measurements obtained were generally within 10% of those obtained from destructive samples measured manually and with a flat-bed scanner. Further developments of the system will allow more detailed examination of the root: soil interface.

Forward modelling to determine the observational signatures of white-light imaging and interplanetary scintillation for the propagation of an interplanetary shock in the ecliptic plane

Recent coordinated observations of interplanetary scintillation (IPS) from the EISCAT, MERLIN, and STELab, and stereoscopic white-light imaging from the two heliospheric imagers (HIs) onboard the twin STEREO spacecraft are significant to continuously track the propagation and evolution of solar eruptions throughout interplanetary space. In order to obtain a better understanding of the observational signatures in these two remote-sensing techniques, the magnetohydrodynamics of the macro-scale interplanetary disturbance and the radio-wave scattering of the micro-scale electron-density fluctuation are coupled and investigated using a newly constructed multi-scale numerical model. This model is then applied to a case of an interplanetary shock propagation within the ecliptic plane. The shock could be nearly invisible to an HI, once entering the Thomson-scattering sphere of the HI. The asymmetry in the optical images between the western and eastern HIs suggests the shock propagation off the Sun–Earth line. Meanwhile, an IPS signal, strongly dependent on the local electron density, is insensitive to the density cavity far downstream of the shock front. When this cavity (or the shock nose) is cut through by an IPS ray-path, a single speed component at the flank (or the nose) of the shock can be recorded; when an IPS ray-path penetrates the sheath between the shock nose and this cavity, two speed components at the sheath and flank can be detected. Moreover, once a shock front touches an IPS ray-path, the derived position and speed at the irregularity source of this IPS signal, together with an assumption of a radial and constant propagation of the shock, can be used to estimate the later appearance of the shock front in the elongation of the HI field of view. The results of synthetic measurements from forward modelling are helpful in inferring the in-situ properties of coronal mass ejection from real observational data via an inverse approach.

Spectroscopic imaging of arteries and atherosclerotic plaques

Fourier transform infrared (FTIR) spectroscopic imaging using a focal plane array detector has been used to study atherosclerotic arteries with a spatial resolution of 3-4 mum, i.e., at a level that is comparable with cellular dimensions. Such high spatial resolution is made possible using a micro-attenuated total reflection (ATR) germanium objective with a high refractive index and therefore high numerical aperture. This micro-ATR approach has enabled small structures within the vessel wall to be imaged for the first time by FTIR. Structures observed include the elastic lamellae of the tunica media and a heterogeneous distribution of small clusters of cholesterol esters within an atherosclerotic lesion, which may correspond to foam cells. A macro-ATR imaging method was also applied, which involves the use of a diamond macro-ATR accessory. This study of atherosclerosis is presented as an illustrative example of the wider potential of these A TR imaging approaches for cardiovascular medicine and biomedical applications. (C) 2004 Wiley Periodicals, Inc.

Using GPUs for beamforming acceleration on SAFT imaging

SAFT techniques are based on the sequential activation, in emission and reception, of the array elements and the post-processing of all the received signals to compose the image. Thus, the image generation can be divided into two stages: (1) the excitation and acquisition stage, where the signals received by each element or group of elements are stored; and (2) the beamforming stage, where the signals are combined together to obtain the image pixels. The use of Graphics Processing Units (GPUs), which are programmable devices with a high level of parallelism, can accelerate the computations of the beamforming process, that usually includes different functions such as dynamic focusing, band-pass filtering, spatial filtering or envelope detection. This work shows that using GPU technology can accelerate, in more than one order of magnitude with respect to CPU implementations, the beamforming and post-processing algorithms in SAFT imaging. ©2009 IEEE.