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

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.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

Hybrid Finite Compact-WENO Schemes for Shock Calculation

Hybrid finite compact (FC)-WENO schemes are proposed for shock calculations. The two sub-schemes (finite compact difference scheme and WENO scheme) are hybridized by means of the similar treatment as in ENO schemes. The hybrid schemes have the advantages of FC and WENO schemes. One is that they possess the merit of the finite compact difference scheme, which requires only bi-diagonal matrix inversion and can apply the known high-resolution flux to obtain high-performance numerical flux function; another is that they have the high-resolution property of WENO scheme for shock capturing. The numerical results show that FC-WENO schemes have better resolution properties than both FC-ENO schemes and WENO schemes. In addition, some comparisons of FC-ENO and artificial compression method (ACM) filter scheme of Yee et al. are also given.

Non-Linear Wave-Induced Transient Response of Soil Around a Trenched Pipeline

Submarine pipelines are always trenched within a seabed for reducing wave loads and thereby enhancing their stability. Based on Biot’s poroelastic theory, a two-dimensional finite element model is developed to investigate non-linear wave-induced responses of soil around a trenched pipeline, which is verified with the flume test results by Sudhan et al. [Sudhan, C.M., Sundar, V., Rao, S.N., 2002. Wave induced forces around buried pipeline. Ocean Engineering, 29, 533–544] and Turcotte et al. [Turcotte, B.R., Liu, P.L.F., Kulhawy, F.H., 1984. Laboratory evaluation of wave tank parameters for wave-sediment interaction. Joseph H. Defree Hydraulic Laboratory Report 84-1, School of Civil and Environmental Engineering, Cornell University]. Non-linear wave-induced transient pore pressure around pipeline at various phases of wave loading is examined firstly. Unlike most previous investigations, in which only a single sediment layer and linear wave loading were concerned, in this study, the influences of the non-linearity of wave loading, the physical properties of backfill materials and the geometry profile of trenches on the excess pore pressures within the soil around pipeline, respectively, were explored, taking into account the in situ conditions of buried pipeline in the shallow ocean zones. Based on the parametric study, it is concluded that the shear modulus and permeability of backfill soils significantly affect the wave-induced excess pore pressures around trenched pipeline, and that the effect of wave non-linearity becomes more pronounced and comparable with that of trench depth, especially at high wave steepness in shallow water.

Governing Equations and Numerical Solutions of Tension Leg Platform with Finite Amplitude Motion

It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP. The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theoretical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displacement, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.

Analysis of instability of a cylindrical soap film of finite dimensions

By means of experiments of instability of a uniform cylindrical soap film, Boys had showed that the bubble molded by the film is unstable when its length is greater than its circumference. Recently that is generally called the Rayleigh Criterion. In this paper, a linear theory in hydrodynamics is applied to analyze the stability of the cylindrical soap film supported by two equal size disks; all conditions of the stationary wave on the end plates of two disks are given. From here we get that the Rayleigh Criterion on the stability of the cylindrical soap film is proved.

Relationships Between Initial Unloading Slope, Contact Depth, and Mechanical Properties for Conical Indentation in Linear Viscoelastic Solids

Using analytical and finite element modeling, we examine the relationships between initial unloading slope, contact depth, and mechanical properties for spherical indentation in viscoelastic solids with either displacement or load as the independent variable. We then investigate whether the Oliver-Pharr method for determining the contact depth and contact radius, originally proposed for indentation in elastic and elastic-plastic solids, is applicable to spherical indentation in viscoelastic solids. Finally, the analytical and numerical results are used to answer questions raised in recent literature about measuring viscoelastic properties from instrumented spherical indentation experiments.

Three-dimensional finite element analysis of composites with coated spherical inclusions

A three-dimensional finite element analysis has been used to determine the internal stresses in a three-phase composite. The stresses have been determined for a variety of interphase properties, the thicknesses of the interphase and the volume fractions of particles. Young's modulus has been calculated from a knowledge of these stresses and the applied deformation. The calculations show that stress distributions in the matrix and the mechanical properties are sensitive to the interphase property in the three-phase composites. The interfacial stresses in the three-dimensional analysis are in agreement with results obtained by an axisymmetric analysis. The predicted bulk modulus in three-dimensional analysis agrees well with the theoretical solution obtained by Qui and Weng, but it presents a great divergence from that in axisymmetric analyses. An investigation indicates that this divergence may be caused by the difference in the unit cell structure between two models. A comparison of the numerically predicted bulk and shear modulus for two-phase composites with the theoretical results indicates that the three-dimensional analysis gives quite satisfactory results.

A Finite Element Analysis of a Crack Penetrating or Deflecting into an Interface in a Thin Laminate

The crack tip driving force of a crack growing from a pre-crack that is perpendicular to and terminating at an interface between two materials is investigated using a linear fracture mechanics theory. The analysis is performed both for a crack penetrating the interface, growing straight ahead, and for a crack deflecting into the interface. The results from finite element calculations are compared with asymptotic solutions for infinitesimally small crack extensions. The solution is found to be accurate even for fairly large amounts of crack growth. Further, by comparing the crack tip driving force of the deflected crack with that of the penetrating crack, it is shown how to control the path of the crack by choosing the adhesion of the interface relative to the material toughness.

Relationships between Initial Unloading Slope, Contact Depth, and Mechanical Properties for Spherical Indentation in Linear Viscoelastic Solids

Using analytical and finite element modeling, we examine the relationships between initial unloading slope, contact depth, and mechanical properties for spherical indentation in viscoelastic solids with either displacement or load as the independent variable. We then investigate whether the Oliver-Pharr method for determining the contact depth and contact radius, originally proposed for indentation in elastic and elastic-plastic solids, is applicable to spherical indentation in viscoelastic solids. Finally, the analytical and numerical results are used to answer questions raised in recent literature about measuring viscoelastic properties from instrumented spherical indentation experiments.

A surface energy model and application to mechanical behavior analysis of single crystals at sub-micron scale

Size effects of mechanical behaviors of materials are referred to the variation of the mechanical behavior due to the sample sizes changing from macroscale to micro-/nanoscales. At the micro-/nanoscale, since sample has a relatively high specific surface area (SSA) (ratio of surface area to volume), the surface although it is often neglected at the macroscale, becomes prominent in governing the energy effect, although it is often neglected at the macroscale, becomes prominent in governing the mechanical behavior. In the present research, a continuum model considering the surface energy effect is developed through introducing the surface energy to total potential energy. Simultaneously, a corresponding finite element method is developed. The model is used to analyze the axial equilibrium strain problem for a Cu nanowire at the external loading-free state. As another application of the model, from dimensional analysis, the size effects of uniform compression tests on the microscale cylinder specimens for Ni and Au single crystals are analyzed and compared with experiments in literatures. (C) 2009 Elsevier B.V. All rights reserved.

激光二极管端面抽运Nd∶YVO4固体激光器热致光束畸变非对称的平衡

由于a轴切割Nd∶YVO4晶体的非对称性,使得激光二极管(LD)端面抽运Nd∶YVO4固体激光器不同于Nd∶YAG激光器,输出的激光经常产生非对称结果。用有限元法分析激光二极管端面抽运a轴切割Nd∶YVO4固体激光器的晶体热效应,包括温度分布、内部应力和产生的形变。分析结果表明端面抽运a轴切割Nd∶YVO4晶体产生了椭球热透镜效应。从结构方面和抽运方面提出了热透镜非对称性的平衡方法,实验验证了方法的可行性。

Rigorous transmittance analysis of a sub-wavelength Sb thin film lens

We quantitatively analysed the factors contributing to the optical transmission enhancement of a sub-wavelength Sb thin film lens, using the finite-difference time-domain (FDTD) method. The results show that the transmission enhancement of the dielectric with a Gaussian distributed refractive index loaded in a sub-wavelength circular hole is not only due to the high refractive index dielectric, but also due to the specific distributions of refractive index. It is the first study about the effects of the refractive index distribution on the transmission of a sub-wavelength aperture. This kind of lens has practical applications in the very small aperture lasers and for near-field optical storage and lithography.

Numerical study of sub-wavelength plasmonic waveguide

The authors present an analysis of a plasmonic waveguide, simulated using a two-dimensional finite-difference time-domain technique. With the surface structures located on the surface of the metal, the device is able to confine and guide light waves in a sub-wavelength scale. And two waveguides can be placed within 150 nm (similar to 6% of the incident wavelength) that will helpful for the optoelectronic integration. Within the 20 mu m simulation region, it is found that the intensity of the guided light at the interface is roughly two to four times the peak intensity of the incident light, and the propagation length can reach approximately 40 Pm at the wavelength of 2.44 mu m. (c) 2007 Elsevier B.V. All rights reserved.

Mode-coupling analysis of three-dimensional microdisk resonators by the finite-difference time-domain technique

Quality factor enhancement due to mode coupling is observed in a three-dimensional microdisk resonator. The microdisk, which is vertically sandwiched between air and a substrate, with a radius of 1 mu m, a thickness of 0.2 mu m, and a refractive index of 3.4, is considered in a finite-difference time-domain (FDTD) numerical simulation. The mode quality factor of the fundamental mode HE71 decreases with an increase of the refractive index of the substrate, n(sub), from 2.0 to 3.17. However, the mode quality factor of the first-order mode HE72 reaches a peak value at n(sub) = 2.7 because of the mode coupling between the fundamental and the first-order modes. The variation of mode field distributions due to the mode coupling is also observed. This mechanism may be used to realize high-quality-factor modes in microdisks with high-refractive-index substrates. (c) 2006 Optical Society of America.