多尺度分解与相位反褶积提高海上地震资料分辨率和信噪比

Offshore seismic exploration is full of high investment and risk. And there are many problems, such as multiple. The technology of high resolution and high S/N ratio on marine seismic data processing is becoming an important project. In this paper, the technology of multi-scale decomposition on both prestack and poststack seismic data based on wavelet and Hilbert-Huang transform and the theory of phase deconvolution is proposed by analysis of marine seismic exploration, investigation and study of literatures, and integration of current mainstream and emerging technology. Related algorithms are studied. The Pyramid algorithm of decomposition and reconstruction had been given by the Mallat algorithm of discrete wavelet transform In this paper, it is introduced into seismic data processing, the validity is shown by test with field data. The main idea of Hilbert-Huang transform is the empirical mode decomposition with which any complicated data set can be decomposed into a finite and often small number of intrinsic mode functions that admit well-behaved Hilbert transform. After the decomposition, a analytical signal is constructed by Hilbert transform, from which the instantaneous frequency and amplitude can be obtained. And then, Hilbert spectrum. This decomposition method is adaptive and highly efficient. Since the decomposition is based on the local characteristics of the time scale of data, it is applicable to nonlinear and non-stationary processes. The phenomenons of fitting overshoot and undershoot and end swings are analyzed in Hilbert-Huang transform. And these phenomenons are eliminated by effective method which is studied in the paper. The technology of multi-scale decomposition on both prestack and poststack seismic data can realize the amplitude preserved processing, enhance the seismic data resolution greatly, and overcome the problem that different frequency components can not restore amplitude properly uniformly in the conventional method. The method of phase deconvolution, which has overcome the minimum phase limitation in traditional deconvolution, approached the base fact well that the seismic wavelet is phase mixed in practical application. And a more reliable result will be given by this method. In the applied research, the high resolution relative amplitude preserved processing result has been obtained by careful analysis and research with the application of the methods mentioned above in seismic data processing in four different target areas of China Sea. Finally, a set of processing flow and method system was formed in the paper, which has been carried on in the application in the actual production process and has made the good progress and the huge economic benefit.

Simple ABCD matrix method for evaluating optical coupling system of laser diode to single-mode fiber with a lensed-tip

In this paper, we present a simple technique to determine the coupling efficiency between a laser diode and a lensed-tip based on the ABCD transformation matrix method. We have compared our analysis technique to that of previous work and have found that the presented method is reliable in predicting the coupling efficiency of lensed-tip and has the advantage of simplicity of coupling efficiency calculation even by a pocket calculator. The results can be useful for designing coupling optics. (c) 2005 Elsevier GmbH. All rights reserved.

Bonding and correlation analysis of various SiCO isomers

The structural properties for various SiCO isomers in the singlet and triplet states have been investigated using CASSCF methods with a 6-311 +G* basis set and also using three DFT and MP2 with same basis set for those systems except for the linear singlet state. The detailed bonding character is discussed, and the state-state correlations and the isomerization mechanism are also determined. Results indicate that there are four different isomers for each spin state, and for all isomers, the triplet state is more stable than the corresponding singlet state. The most stable is the linear SiCO ((3)Sigma(-)) species and may be refer-red to the ground state. At the CASSCF-MP2(full)/6-311+G* level, the state-state energy separations of the other triplet states relative to the ground state are 43.2 (cyclic), 45.2 (linear SiOC), and 75.6 kcal/mol (linear CSiO), respectively, whereas the triplet-singlet state excitation energies for each configuration are 17.3 (linear SiCO), 2.2 (cyclic SiCO), 10.2 (linear SiOC), and 18.5 kcal/mol (linear CSiO), respectively. SiCo ((3)Sigma(-)) may be classified as silene (carbonylsilene), and its COdelta- moiety possesses CO- property. The dissociation energy of the ground state is 42.5 kcal/mol at the CASSCF-MP2(full)/6-311+G* level and falls within a range of 36.5-41.5 kcal/mol at DFT level, and of 23.7-28.9 kcal/mol at the wave function-correlated level, whereas the vertical IP is 188.8 kcal/mol at the CASSCF-MP2(full)/6-311+G* level and is very close to the first IP of Si atom. Three linear isomers (SiCO, SiOC, and CSiO) have similar structural bonding character. SiOC may be referred to the iso-carbonyl Si instead of the aether compound, whereas the CSiO isomer may be considered as the combination of C (the analogue of Si) with SiO (the analogue of CO). The bonding is weak for all linear species, and the corresponding potential energy surfaces are flat, and thus these linear molecules are facile. Another important isomer is of cyclic structure, it may be considered as the combination of CO with Si by the side pi bond. This structure has the smallest triplet state-singlet state excitation energy (similar to2.2 kcal/mol); the C-O bonds are longer, and the corresponding vibrational frequencies are significantly smaller than those of the other linear species. This cyclic species is not classified as an epoxy compound. State-state correlation analysis and the isomerization pathway searches have indicated that there are no direct correlations among three linear structures for each spin state, but they may interchange by experiencing two transition states and one cyclic intermediate. The easiest pathway is to break the Si-O bond to go to the linear SiCO, but its inverse process is very difficult. The most difficult process is to break the C-O bond and to go to the linear CSiO.

Correlation analysis of the structures and gas-liquid chromatographic retention indices of amines

Quantitative structure-retention relationship(QSRR) was studied for amines to gas-liquid chromatography on three stationary phases of different polarities with the topological indices A(m) (A(m1), A(m2), A(m3)) and gravitational index GI. The algorithm of "Leaps and Bounds" was performed for selection of the variables. And the multi-regression and the quasi-Newton neural networks were employed for the calculation with better results.

CORRELATION-ANALYSIS IN STRUCTURE AND CHROMATOGRAPHIC DATA OF ORGANOPHOSPHORUS COMPOUNDS BY GAI

This approach is undertaken to examine the correlation ability of the general a(N)-index (GAI) to predict chromatographic behavior. The test is performed on various types of organophosphorus compounds. The results demonstrate that the GAI possesses a good correlation with chromatographic properties.

Analysis of the sea-level variability along the Chinese coast and estimation of the impact of a CO2-perturbed atmospheric circulation

The relationship between monthly sea-level data measured at stations located along the Chinese coast and concurrent large-scale atmospheric forcing in the period 1960-1990 is examined. It is found that sea-level varies quite coherently along the whole coast, despite the geographical extension of the station set. A canonical correlation analysis between sea-level and sea-level pressure (SLP) indicates that a great part of the sea-level variability can be explained by the action of the wind stress on the ocean surface. The relationship between sea-level and sea-level pressure is analyzed separately for the summer and winter half-years. In winter, one factor affecting sea-level variability at all stations is the SLP contrast between the continent and the Pacific Ocean, hence the intensity of the winter Monsoon circulation. Another factor that affects coherently all stations is the intensity of the zonal circulation at mid-latitudes. In the summer half year, on the other hand, the influence of SLP on sea-level is spatially less coherent: the stations in the Yellow Sea are affected by a more localized circulation anomaly pattern, whereas the rest of the stations is more directly connected to the intensity of the zonal circulation. Based on this analysis, statistical models (different for summer and winter) to hindcast coastal sealevel anomalies from the large-scale SLP field are formulated. These models have been tested by fitting their internal parameters in a test period and reproducing reasonably the sea-level evolution in an independent period. These statistical models are also used to estimate the contribution of the changes of the atmospheric circulation on sea-level along the Chinese coast in an altered climate. For this purpose the ouput of 150 year-long experiment with the coupled ocean-atmosphere model ECHAM1-LSG has been analyzed, in which the atmospheric concentration of greenhouse gases was continuously increased from 1940 until 2090, according to the Scenario A projection of the Intergovermental Panel on Climate Change. In this experiment the meridional (zonal) circulation relevant for sea-level tends to become weaker (stronger) in the winter half year and stronger (weaker) in summer. The estimated contribution of this atmospheric circulation changes to coastal sea-level is of the order of a few centimeters at the end of the integration, being in winter negative in the Yellow Sea and positive in the China Sea with opposite signs in the summer half-year.

中高层大气对低层大气声重波响应的分析

Acoustic Gravity waves (AGW) play an important role in balancing the atmospheric energy and momentum budget. Propagation of gravity wave in the atmosphere is one of the important factors of changing middle and upper atmosphere and ionosphere. The purpose of this dissertation is to study the propagation of gravity wave in a compression atmosphere whit means of numerical simulation and to analyze the response of middle and upper atmosphere to pulse disturbance from lower atmosphere. This work begins with the establishment of 2-D fully nonlinear compressible atmospheric dynamic model in polar coordinate, which is used ton numerically study gravity wave propagation. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. We also simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model and analyze the data we obtained by using Fourier Transform (FT), Short-time Fourier Transform (STFT) and Empirical Mode Decomposition (EMD) method which is an important part of Hilbert-Huang Transform (HHT). The research content is summarized in the following: 1. By using a two-dimensional full-implicit-continuous-Eulerian (FICE) scheme and taking the atmospheric basic motion equations as the governing equations, a numerical model for nonlinear propagation of acoustic gravity wave disturbance in two-dimensional polar coordinates is solved. 2. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. Results of numerical simulation show that the acoustic gravity wave packets propagate steadily upward and keep its shape well after several periods. 3. We simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model, and obtain the distribution of a certain physical quantity in time and space from earth’s surface to 300km above. The results reveal that the response of ionosphere occurs at a large horizontal distance from the source and the disturbance becomes greater with increasing of height. The situation when the direction of the background wind is opposite to or the same as the direction of disturbed velocity of gravity-wave is studied. The results show that gravity wave propagating against winds is easier than those propagating along winds and the background wind can accelerate gravity wave propagation. Just upon the source, an acoustic wave component with period of 6 min can be found. These images of simulation are similar to observations of the total electron content (TEC) disturbances caused by the great Sumatra-Andaman earthquake on December 26 in 2004. 4. Using the EMD method the disturbed velocity data of a certain physical quantity in time and space can be decomposed into a series of intrinsic mode function (IMF) and a trend mode respectively. The results of EMD reveal impact of the gravity wave frequency under the background winds.

HHT在电磁信号处理中的应用

In modem signal Processing，non－linear，non－Gaussian and non－stable signals are usually the analyzed and Processed objects，especially non－stable signals. The convention always to analyze and Process non－stable signals are： short time Fourier transform，Wigner－Ville distribution，wavelet Transform and so on. But the above three algorithms are all based on Fourier Transform，so they all have the shortcoming of Fourier Analysis and cannot get rid of the localization of it. Hilbert－Huang Transform is a new non－stable signal processing technology，proposed by N. E. Huang in 1998. It is composed of Empirical Mode Decomposition (referred to as EMD) and Hilbert Spectral Analysis (referred to as HSA). After EMD Processing，any non－stable signal will be decomposed to a series of data sequences with different scales. Each sequence is called an Intrinsic Mode Function (referred to as IMF). And then the energy distribution plots of the original non－stable signal can be found by summing all the Hilbert spectrums of each IMF. In essence，this algorithm makes the non－stable signals become stable and decomposes the fluctuations and tendencies of different scales by degrees and at last describes the frequency components with instantaneous frequency and energy instead of the total frequency and energy in Fourier Spectral Analysis. In this case，the shortcoming of using many fake harmonic waves to describe non－linear and non－stable signals in Fourier Transform can be avoided. This Paper researches in the following parts： Firstly，This paper introduce the history and development of HHT，subsequently the characters and main issues of HHT. This paper briefly introduced the basic realization principles and algorithms of Hilbert－Huang transformation and confirms its validity by simulations. Secondly, This paper discuss on some shortcoming of HHT. By using FFT interpolation, we solve the problem of IMF instability and instantaneous frequency undulate which are caused by the insufficiency of sampling rate. As to the bound effect caused by the limitation of envelop algorithm of HHT, we use the wave characteristic matching method, and have good result. Thirdly, This paper do some deeply research on the application of HHT in electromagnetism signals processing. Based on the analysis of actual data examples, we discussed its application in electromagnetism signals processing and noise suppression. Using empirical mode decomposition method and multi-scale filter characteristics can effectively analyze the noise distribution of electromagnetism signal and suppress interference processing and information interpretability. It has been founded that selecting electromagnetism signal sessions using Hilbert time-frequency energy spectrum is helpful to improve signal quality and enhance the quality of data.

A Q-swicthed all-solid-state single-longitudinal-mode laser with adjustable pulse-width and high repetition rate

A single-longitudinal-mode (SLM) laser-diode pumped Nd: YAG laser with adjustable pulse width is developed by using the techniques of pre-lasing and changing polarization of birefingent crystal. The Q-switching voltage is triggered by the peak of the pre-lasing pulse to achieve the higher stability of output pulse energy. The output energy of more than 1 mJ is obtained with output energy stability of 3% (rms) at 100 Hz. The pulse-width can be adjusted from 30 ns to 300 ns by changing the Q-switching voltage. The probability of putting out single-longitudinal-mode pulses is almost 100%. The laser can be run over four hours continually without mode hopping.

Environmental mechanism of change in cyanobacterial species composition in the northeastern part of Lake Dianchi (China)

Environmental mechanism of change in cyanobacterial species composition in the northeastern part of Lake Dianchi (also called Macun Bay and Haidong Bay) was studied using canonical correlation analysis (CCA), but also bottom-up control and top-down control were fully discussed. Results from CCA suggest: (1) the abundance and dominance of Microcystis aeruginosa in Macun Bay and Haidong Bay are influenced by total phosphorus (TP), nitrate (NO3--N), nitrite (NO2--N), dissolved oxygen (DO) and water temperature (WT); (2) water temperature has a positive correlation with the abundance of M. aeruginosa and it also has negative correlations with the abundances of Anabaena flos-aquae and Aphanizomenonon flos-aquae; and (3) abundances of both Anabaena flos-aquae and Aphanizomenon flos-aquae have positive correlations with ammonia-N (NH4+-N). Furthermore, cyanobacterial species composition has no significant correlations with light and size-fractioned iron in this study. Grazers, cyanophages and viruses were able to control cyanobacterial blooms and change the composition of cyanobacterial species. Though we studied physical and chemical factors intensely enough, we still are not able to predict the change in the composition of cyanobacterial blooms, because of plankton system in a chaotic behavior.

Simultaneous Q switching and mode locking with narrow envelope width in a microchip Nd : YVO4 laser with a composite semiconductor absorber

By using a composite semiconductor absorber and an output coupler, we demonstrated a Q-switched and mode-locked diode-pumped microchip Nd:YVO4 laser. With a 350-mu m-thick crystal, the width of the Q-switched envelope was as short as 12 ns; the repetition rate of the mode-locked pulses inside the Q-switched pulse was more than 10 GHz. The average output power was 335 mW at a maximum pump power of 1.6 W. Q-switched envelope widths of 21 and 31 ns were also achieved with crystals 0.7 and 1.0 mm thick, respectively.

Correlation analysis of the structures and stability constants of gadolinium(III) complexes

To simplify the abstraction of descriptors, for the correlation analysis of the stability constants of gadolinium(III) complexes and their ligand structures, aiming at gadolinium(III) complexes, we only considered the ligands and ignored the common parts of the structures, i.e., the metal ions. Quantum-chemical descriptors and topological indices were calculated to describe the structures of the ligands. Multiple regression analysis and neural networks were applied to construct the models between the ligands and the stability constants of gadolinium(III) complexes and satisfactory results were obtained.

Correlation analysis of the structures and pK(a) values of carboxylic acids

The quantum-chemical descriptors were used for QSPR study of the structures of carboxylic acids and their pK(a) values. The algorithm of "Leaps and Bounds" regression was performed for selection of the variables. The CoMFA method was carried out for 3D-QSPR. As the introduction of the charge of oxygen atom(Q(2)), the results obtained by CoMFA were improved greatly.

Simultaneous spectrophotometric determination of zinc, cadmium and mercury with hybrid linear analysis

A new method for simultaneous spectrophotometric determination of Zn, Cd and Hg using 2-(5-Br-2-pyridylazo)-5-diethylaminophenol as the color developing reagent was proposed. The absorption spectra of these three complexes have similar features with severe overlap in visible spectral range. For resolving these spectra, hybrid linear analysis was used, and the pure spectrum of each component was obtained from the calibration mixtures by least squares method. The effects of reaction condition, selection of wavelengths, determination of pure spectrum and additivity of absorbances etc. on the determination were discussed. The proposed method offers the advantages of simple, rapid, and accuracy. It has been successfully applied to the simultaneous determination of Zn, Cd and Hg in synthetic sample. A comparison was also made with the partial least squares method.