Phase-space analysis of wavefront coding imaging systems

We explore the use of the Radon-Wigner transform, which is associated with the fractional Fourier transform of the pupil function, for determining the point spread function (PSF) of an incoherant defocused optical system. Then we introduce these phase-space tools to analyse the wavefront coding imaging system. It is shown that the shape of the PSF for such a system is highly invarient to the defocous-related aberrations except for a lateral shift. The optical transfer function of this system is also investigated briefly from a new understanding of ambiguity function.

A new method for the estimation of oceanic mixed-layer depth using shipboard X-band radar images

Shipboard X-band radar images acquired on 24 June 2009 are used to study nonlinear internal wave characteristics in the northeastern South China Sea. The studied images show three nonlinear internal waves in a packet. A method based on the Radon Transform technique is introduced to calculate internal wave parameters such as the direction of propagation and internal wave velocity from backscatter images. Assuming that the ocean is a two-layer finite depth system, we can derive the mixed-layer depth by applying the internal wave velocity to the mixed-layer depth formula. Results show reasonably good agreement with in-situ thermistor chain and conductivity-temperature-depth data sets.

一种快速获取激光拼焊焊缝错配的新算法

本文讨论了基于视觉和结构光的焊缝形貌视觉检测系统的组成原理,针对拼板激光焊接工程中对焊缝检测的实时性要求高的特点,创新性地提出了一种快速获得焊缝错配缺陷的算法,主要步骤为:首先基于有关标准的检测指标阈值的设定,再次,针对在线提取到的图像进行预处理,主要是加窗和中值滤波;最后为Radon变换与错配检测。该算法能减少计算任务,从而快速获得错配检测指标。实验给出了不等厚板拼焊时的线性错配的数值和分布,同时验证了该算法的有效性。

叠前地震数据规则化、重建及噪音压制

I address of reconstruction of spatial irregular sampling seismic data to regular grids. Spatial irregular sampling data impairs results of prestack migration, multiple attenuations, spectra estimation. Prestack 5-D volumes are often divided into sub-sections for further processing. Shot gathers are easy to obtain from irregular sampling volumes. My strategy for reconstruction is as follows: I resort irregular sampling gathers into a form of easy to bin and perform bin regularization, then utilize F-K inversion to reconstruct seismic data. In consideration of poor ability of F-K regularization to fill in large gaps, I sort regular sampling gathers to CMP and proposed high-resolution parabolic Radon transform to interpolate data and extrapolate offsets. To strong interfering noise--multiples, I use hybrid-domain high-resolution parabolic Radon transform to attenuate it. F-K regularization demand ultimately for lower computing costs. I proposed several methods to further improve efficiency of F-K inversion: first I introduce 1D and 2D NFFT algorithm for a rapid calculation of DFT operators; then develop fast 1D and 2D CG method to solve least-square equations, and utilize preconditioner to accelerate convergence of CG iterations; what’s more, I use Delaunay triangulation for weight calculation and use bandlimit frequency and varying bandwidth technique for competitive computation. Numerical 2D and 3D examples are offered to verify reasonable results and more efficiency. F-K regularization has poor ability to fill in large gaps, so I rearrange data as CMP gathers and develop hybrid-domain high-resolution parabolic Radon transforms which be used ether to interpolate null traces and extrapolate near and far offsets or suppress a strong interfere noise: multiples. I use it to attenuate multiples to verify performances of our algorithm and proposed routines for industrial application. Numerical examples and field data examples show a nice performance of our method.

叠前时间偏移：各向异性、多次波衰减及应用

At present, in order to image complex structures more accurately, the seismic migration methods has been developed from isotropic media to the anisotropic media. This dissertation develops a prestack time migration algorithm and application aspects for complex structures systematically. In transversely isotropic media with a vertical symmetry axis (VTI media), the dissertation starts from the theory that the prestack time migration is an approximation of the prestack depth migration, based on the one way wave equation and VTI time migration dispersion relation, by combining the stationary-phase theory gives a wave equation based VTI prestack time migration algorithm. Based on this algorithm, we can analytically obtain the travel time and amplitude expression in VTI media, as while conclude how the anisotropic parameter influence the time migration, and by analyzing the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, we can update the velocity model and estimate the anisotropic parameter model through the time migration. When anisotropic parameter is zero, this algorithm degenerates to the isotropic time migration algorithm naturally, so we can propose an isotopic processing procedure for imaging. This procedure may keep the main character of time migration such as high computational efficiency and velocity estimation through the migration, and, additionally, partially compensate the geometric divergence by adopting the deconvolution imaging condition of wave equation migration. Application of this algorithm to the complicated synthetic dataset and field data demonstrates the effectiveness of the approach. In the dissertation we also present an approach for estimating the velocity model and anisotropic parameter model. After analyzing the velocity and anisotropic parameter impaction on the time migration, and based on the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, through migration we can update the velocity model and estimate the anisotropic parameter model by combining the advantages of velocity analysis in isotropic media and anisotropic parameter estimation in VTI media. Testing on the synthetic and field data, demonstrates the method is effective and very steady. Massive synthetic dataset、2D sea dataset and 3D field datasets are used for VTI prestack time migration and compared to the stacked section after NMO and prestack isotropic time migration stacked section to demonstrate that VTI prestack time migration method in this paper can obtain better focusing and less positioning errors of complicated dip reflectors. When subsurface is more complex, primaries and multiples could not be separated in the Radon domain because they can no longer be described with simple functions (parabolic). We propose an attenuating multiple method in the image domain to resolve this problem. For a given velocity model，since time migration takes the complex structures wavefield propagation in to account, primaries and multiples have different offset-domain moveout discrepancies, then can be separated using techniques similar to the prior migration with Radon transform. Since every individual offset-domain common-reflection point gather incorporates complex 3D propagation effects, our method has the advantage of working with 3D data and complicated geology. Testing on synthetic and real data, we demonstrate the power of the method in discriminating between primaries and multiples after prestack time migration, and multiples can be attenuated in the image space considerably.

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

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.

中国地下工程氡污染健康影响与防护研究-以内蒙古自治区为例

由于对人体健康的危害，氡越来越受到世界各国的广泛关注。经济的发展，使我国地下公共场所大量涌现，但其空气中氡浓度的大规模调查及系统研究还属空白。本论文采取瞬时法（美国SUN公司的1027型连续测氰仪）与累积法（CR一39SSNTD）相结合的测量方法，对我国23个城市的地下工程进行了为期一年的监测，获得了实测数据700多个。在监测的同时，收集了每个测点的各种相关情况资料。论文以内蒙古自治区的呼和浩特、包头、集宁、鄂尔多斯四个城市为例，从环境地球化学角度出发，在研究了各个测点的各种情况之后，探讨了氧在不同类型工程和不同城市的浓度规律。对长期有人在内工作的地下工程内氛所致工作人员的有效剂量当量进行了估算。最后对氨异常的测点进行了原因分析，并提出了地下工程内氰污染的治理和防范措施。本次研究主要得出以下结论：1、地下工程中氰浓度的季节变化规律为：从总平均值看，夏高冬低。与前人对坑道的调查结果一致，与地上氨浓度的变化规律相反。从两季浓度差来看，两季浓度差从大到小依次为：隧道类工程）单建式铜室类工程＞结建式嗣室类工程。2、总体来看内蒙三类工程的氰浓度值：隧道类工程）结建式铜室类工程）单建式铜室类工程。不管是单季还是年平均值，各城市三类工程的氧浓度都遵循隧道类．工程远大于结建式铜室工程和单建式铜室工程的规律。通过对结建式桐室工程和单建式铜室工程的对比分析，发现在地下工程内工程上部有无建筑物阻隔，并不是影响氨浓度高低的重要因素，即氰的纵向扩散作用小于横向扩散作用。3、呼和浩特、包头、集宁、鄂尔多斯四个城市的氛浓度年平均值分别为：190Bq／m3、700Bq／m3、190Bq／m3、270Bq／m3。去掉仅有一季结果的侧点值，包头、鄂尔多斯的氨浓度年平均值为：62oBq／m3、250Bq／m3。呼和浩特、包头、集宁、鄂尔多斯四个城市的氛浓度单季值的范围分别为：31-830Bq／m3、57-3500Bq／m3、26一65oBq／m，、64一69OBq／m'，；年平均值的范围分别为：39一580Bq／m：'、94一210oBq／m：，、35-460Bq／m3、110-390Bq／m3。呼和浩特、包头、集宁、鄂尔多斯四个城市分别有30.0%、44.4%、33.3%、50.0%的工程年平均氧浓度在200Bq／时以上。四个城市氦浓度的普遍水平：包头）鄂尔多斯）集宁）呼和浩特。4、对所有单建式铜室类地下工程和结建式铜室类工程内工作人员所受剂量当量进行了估算。单建式地下工程内工作人员所受剂量当量的范围是1.247mSv-6．SOgmsv，平均值是2.610Sv。结建式地下工程内工作人员所受剂量当量的范围是1．195Sv-8.531mSv，平均值是2.572mSv。两类地下工程内工作人员所受氨照射剂量当量都很高，有一半的工程剂量当量是地面上公众人员剂量当量0.924mSv的两倍以上。5、不管是从浓度角度还是从剂量估算角度分析，地下工程内的氛水平都比较偏高，需要采取一定防护措施，对现有状况进行改善。6、通风是地下工程防治氡污染最简便易行、见效最快的措施。使用防氡涂料是最经济的措施。