Plasmonic surface-wave splitter

The authors present an analysis of a plasmonic surface-wave splitter, simulated using a two-dimensional finite-difference time-domain technique. A single subwavelength slit is employed as a high-intensity nanoscale excitation source for plasmonic surface waves, resulting in a miniaturized light-surface plasmon coupler. With different surface structures located on the two sides of the slit, the device is able to confine and guide light waves of different wavelengths in opposite directions. Within the 15 mu m simulation region, it is found that the intensity of the guided light at the interface is roughly two to eight times the peak intensity of the incident light, and the propagation length can reach approximately 42 and 16 mu m and at the wavelengths of 0.63 and 1.33 mu m, respectively. (c) 2007 American Institute of Physics.

Effect of disorder on self-collimated beam in photonic crystal

We report a numerical analysis of various types of disorder effects on self-collimated beam in two-dimensional photonic crystal. Finite-difference time-domain (FDTD) method is used to simulate the process by using a pulse propagation technique. The position disorders along the directions parallel and perpendicular to the incidence are considered. We show that random disorder along the perpendicular direction will have a lesser effect on the performance of the dispersion waveguides than those along the parallel direction. Furthermore, the self-collimation waveguide (SCW) has new characteristics when compared with the photonic crystal line defect waveguide. (c) 2006 Elsevier B.V. All rights reserved.

A differential recursive scheme for suppression of Peak to average power ratio for OFDM

A differential recursive scheme for suppression of Peak to average power ratio (PAPR) for Orthogonal frequency division multiplexing (OFDM) signal is proposed in this thesis. The pseudo-randomized modulating vector for the subcarrier series is differentially phase-encoded between successive components in frequency domain first, and recursion manipulates several samples of Inverse fast Fourier transformation (IFFT) output in time domain. Theoretical analysis and experimental result exhibit advantage of differential recursive scheme over direct output scheme in PAPR suppression. And the overall block diagram of the scheme is also given.

Resonant frequencies and quality factors for optical equilateral triangle resonators calculated by FDTD technique and the Pade approximation

The mode wavelength and quality factor (Q-factor) for resonant modes in optical equilateral triangle resonators (ETR's) are calculated by the finite-difference time-domain (FDTD) technique and the Pade approximation, For an ETR with the side length of 3 mu m and the refractive index of 3.2, we get the mode wavelength interval of about 70 nm and the Q-factor of the fundamental mode over 10(3), The results show that the ETR is suitable to realize single-mode operation, and that the radiation loss in the corner regions of ETR is rather low, In addition, the numerical results of the mode wavelength agree very well with our analytical formula.

Experimental investigation of slow light phenomenon in photonic crystal waveguide line defect laser

The guide mode whose frequency locates in the band edge in photonic crystal single line defect waveguide has very low group velocity. So the confinement and gain of electromagnetic field in the band edge are strongly enhanced. Photonic crystal waveguide laser is fabricated and the slow light phenomenon is investigated. The laser is pumped by pulsed pumping light at 980nm whose duty ratio is 0.05%. The active layer in photonic crystal slab is InGaAsP multiple quantum well. Light is transimited by a photonic crystal chirp waveguide in one facet of the laser. Then the output light is coupled to a fiber and the character of laser is analysis by an optical spectrometer. It is found that single mode and multimode happens with different power of pumping light. Meanwhile the plane wave expansion and finite-difference time-domain methods are used to simulate the phenomenon of slow light. And the result of the experiment is compared with the theory which proves the slow light results in lasing oscillation.

A serpin from Chinese shrimp Fenneropenaeus chinensis is responsive to bacteria and WSSV challenge

Arthropod defence responses (e.g. prophenoloxidase (proPO) activation and Toll pathway initiation) are mediated by serine proteinase cascades and regulated by serpins in haemolymph. A serpin (Fc-serpin) cDNA was cloned from the haemocytes of Fenneropenaeus chinensis by rapid amplification of cDNA ends (RACE) PCR and haemocyte cDNA library screening. The full-length cDNA consists of 1734 bp, encoding 411 amino acids with a calculated molecular mass of 46.55 kDa and a theoretical isoelectric point of 7.70. Fc-serpin contains a typical serpin-like homologue (serine proteinase inhibitors domain). The deduced protein contains a putative signal peptide of 19 amino acids and the serpin's signature sequence ((FHCNRPFLFLI389)-F-379). Fc-serpin showed some identity with Pacifastacus leniusculus serpin (42%) and Manduca sexta serpin-6 (34%). The reactive centre loop (RCL) sequences of Fc-serpin, P leniusculus serpin, M. sexta serpin-6 and Bombyx mori serpin-2 are highly similar. An Arg at the PI position of the reactive site indicates that Fc-serpin may have inhibitory activity against prophenoloxidase activating proteinase (PAP) and clotting enzyme. Transcripts of Fc-serpin mRNA were mainly detected in haemocytes and the lymphoid organ by RT-PCR. The variation of the mRNA transcription level in haemocytes followed by artificial infection with bacteria OF white spot syndrome virus (WSSV) was quantified by SYBR Green real-time PCR analysis. Expression profiles of Fc-serpin greatly fluctuated after challenge. This work represents the first report Of a serpin in penaeid shrimp. The data provide clues that Fc-serpin might play potential roles in the innate immunity of shrimp. (C) 2008 Elsevier Ltd. All rights reserved.

Characterization of an ETS transcription factor in the sea scallop Chlamys farreri

We have cloned and characterized a cDNA encoding a putative ETS transcription factor, designated Cf-ets. The Cf-ets encodes a 406 amino acid protein containing a conserved ETS domain and a Pointed domain. Phylogenetic analysis revealed that Cf-ets belongs to the ESE group of ETS transcription factor family. Real-time PCR analysis of Cf-ets expression in adult sea scallop tissues revealed that Cf-ets was expressed mainly in gill and hemocytes, in a constitutive manner. Cf-ets mRNA level in hemocytes increased drastically after microbial challenge indicated its indispensable role in the anti-infection process. Simultaneously, the circulating hemocyte number decreased. In mammals, most ETS transcription factors play indispensable roles in blood cell differentiation and linage commitment during hematopoisis. Cf-ets is therefore likely to be a potential biomarker for hematopoiesis studies in scallops. (C) 2009 Elsevier Ltd. All rights reserved.

Molecular cloning, expression and biochemical property analysis of AtKP1, a kinesin gene from Arabidopsis thaliana

Kinesins are common in a variety of eukaryotic cells with diverse functions. A cDNA encoding a member of the Kinesin-14B subfamily is obtained using X-RACE technology and named AtKP1 (for Arabidopsis kinesin protein 1). This cDNA has a maximum open reading frame of 3.3 kb encoding a polypeptide of 1087 aa. Protein domain analysis shows that AtKP1 contains the motor domain and the calponin homology domain in the central and amino-terminal regions, respectively. The carboxyl-terminal region with 202 aa residues is diverse from other known kinesins. Northern blot analysis shows that AtKP1 is widely expressed at a higher level in seedlings than in mature plants. 2808 bp of the AtKP1 promoter region is cloned and fused to GUS. GUS expression driven by the AtKP1 promoter region shows that AtKP1 is mainly expressed in vasculature of young organs and young leaf trichomes, indicating that AtKP1 may participate in the differentiation or development of Arabidopsis thaliana vascular bundles and trichomes. A truncated AtKP1 protein containing the putative motor domain is expressed in E. coli and affinity-purified. In vitro characterizations indicate that the polypeptide has nucleotide-dependent microtubule-binding ability and microtubule-stimulated ATPase activity.

基于匹配追踪方法的地震信号时频分析研究

Seismic signal is a typical non-stationary signal, whose frequency is continuously changing with time and is determined by the bandwidth of seismic source and the absorption characteristic of the media underground. The most interesting target of seismic signal’s processing and explaining is to know about the local frequency’s abrupt changing with the time, since this kind of abrupt changing is indicating the changing of the physical attributes of the media underground. As to the seismic signal’s instantaneous attributes taken from time-frequency domain, the key target is to search a effective, non-negative and fast algorithm time-frequency distribution, and transform the seismic signal into this time-frequency domain to get its instantaneous power spectrum density, and then use the process of weighted adding and average etc. to get the instantaneous attributes of seismic signal. Time-frequency analysis as a powerful tool to deal with time variant non-stationary signal is becoming a hot researching spot of modern signal processing, and also is an important method to make seismic signal’s attributes analysis. This kind of method provides joint distribution message about time domain and frequency domain, and it clearly plots the correlation of signal’s frequency changing with the time. The spectrum decomposition technique makes seismic signal’s resolving rate reach its theoretical level, and by the method of all frequency scanning and imaging the three dimensional seismic data in frequency domain, it improves and promotes the resolving abilities of seismic signal vs. geological abnormal objects. Matching pursuits method is an important way to realize signal’s self-adaptive decomposition. Its main thought is that any signal can be expressed by a series of time-frequency atoms’ linear composition. By decomposition the signal within an over completed library, the time-frequency atoms which stand for the signal itself are selected neatly and self-adaptively according to the signal’s characteristics. This method has excellent sparse decomposition characteristics, and is widely used in signal de-noising, signal coding and pattern recognizing processing and is also adaptive to seismic signal’s decomposition and attributes analysis. This paper takes matching pursuits method as the key research object. As introducing the principle and implementation techniques of matching pursuits method systematically, it researches deeply the pivotal problems of atom type’s selection, the atom dictionary’s discrete, and the most matching atom’s searching algorithm, and at the same time, applying this matching pursuits method into seismic signal’s processing by picking-up correlative instantaneous messages from time-frequency analysis and spectrum decomposition to the seismic signal. Based on the research of the theory and its correlative model examination of the adaptively signal decomposition with matching pursuit method, this paper proposes a fast optimal matching time-frequency atom’s searching algorithm aimed at seismic signal’s decomposition by frequency-dominated pursuit method and this makes the MP method pertinence to seismic signal’s processing. Upon the research of optimal Gabor atom’s fast searching and matching algorithm, this paper proposes global optimal searching method using Simulated Annealing Algorithm, Genetic Algorithm and composed Simulated Annealing and Genetic Algorithm, so as to provide another way to implement fast matching pursuit method. At the same time, aimed at the characteristics of seismic signal, this paper proposes a fast matching atom’s searching algorithm by means of designating the max energy points of complex seismic signal, searching for the most optimal atom in the neighbor area of these points according to its instantaneous frequency and instantaneous phase, and this promotes the calculating efficiency of seismic signal’s matching pursuit algorithm. According to these methods proposed above, this paper implements them by programmed calculation, compares them with some open algorithm and proves this paper’s conclusions. It also testifies the active results of various methods by the processing of actual signals. The problems need to be solved further and the aftertime researching targets are as follows: continuously seeking for more efficient fast matching pursuit algorithm and expanding its application range, and also study the actual usage of matching pursuit method.

积分偏移的计算几何与并行实现

The Second Round of Oil & Gas Exploration needs more precision imaging method, velocity vs. depth model and geometry description on Complicated Geological Mass. Prestack time migration on inhomogeneous media was the technical basic of velocity analysis, prestack time migration on Rugged surface, angle gather and multi-domain noise suppression. In order to realize this technique, several critical technical problems need to be solved, such as parallel computation, velocity algorithm on ununiform grid and visualization. The key problem is organic combination theories of migration and computational geometry. Based on technical problems of 3-D prestack time migration existing in inhomogeneous media and requirements from nonuniform grid, parallel process and visualization, the thesis was studied systematically on three aspects: Infrastructure of velocity varies laterally Green function traveltime computation on ununiform grid, parallel computational of kirchhoff integral migration and 3D visualization, by combining integral migration theory and Computational Geometry. The results will provide powerful technical support to the implement of prestack time migration and convenient compute infrastructure of wave number domain simulation in inhomogeneous media. The main results were obtained as follows: 1. Symbol of one way wave Lie algebra integral, phase and green function traveltime expressions were analyzed, and simple 2-D expression of Lie algebra integral symbol phase and green function traveltime in time domain were given in inhomogeneous media by using pseudo-differential operators’ exponential map and Lie group algorithm preserving geometry structure. Infrastructure calculation of five parts, including derivative, commutating operator, Lie algebra root tree, exponential map root tree and traveltime coefficients , was brought forward when calculating asymmetry traveltime equation containing lateral differential in 3-D by this method. 2. By studying the infrastructure calculation of asymmetry traveltime in 3-D based on lateral velocity differential and combining computational geometry, a method to build velocity library and interpolate on velocity library using triangulate was obtained, which fit traveltime calculate requirements of parallel time migration and velocity estimate. 3. Combining velocity library triangulate and computational geometry, a structure which was convenient to calculate differential in horizontal, commutating operator and integral in vertical was built. Furthermore, recursive algorithm, for calculating architecture on lie algebra integral and exponential map root tree (Magnus in Math), was build and asymmetry traveltime based on lateral differential algorithm was also realized. 4. Based on graph theory and computational geometry, a minimum cycle method to decompose area into polygon blocks, which can be used as topological representation of migration result was proposed, which provided a practical method to block representation and research to migration interpretation results. 5. Based on MPI library, a process of bringing parallel migration algorithm at arbitrary sequence traces into practical was realized by using asymmetry traveltime based on lateral differential calculation and Kirchhoff integral method. 6. Visualization of geological data and seismic data were studied by the tools of OpenGL and Open Inventor, based on computational geometry theory, and a 3D visualize system on seismic imaging data was designed.

粘弹性介质中的地震波传播与波形反演研究

The real earth is far away from an ideal elastic ball. The movement of structures or fluid and scattering of thin-layer would inevitably affect seismic wave propagation, which is demonstrated mainly as energy nongeometrical attenuation. Today, most of theoretical researches and applications take the assumption that all media studied are fully elastic. Ignoring the viscoelastic property would, in some circumstances, lead to amplitude and phase distortion, which will indirectly affect extraction of traveltime and waveform we use in imaging and inversion. In order to investigate the response of seismic wave propagation and improve the imaging and inversion quality in complex media, we need not only consider into attenuation of the real media but also implement it by means of efficient numerical methods and imaging techniques. As for numerical modeling, most widely used methods, such as finite difference, finite element and pseudospectral algorithms, have difficulty in dealing with problem of simultaneously improving accuracy and efficiency in computation. To partially overcome this difficulty, this paper devises a matrix differentiator method and an optimal convolutional differentiator method based on staggered-grid Fourier pseudospectral differentiation, and a staggered-grid optimal Shannon singular kernel convolutional differentiator by function distribution theory, which then are used to study seismic wave propagation in viscoelastic media. Results through comparisons and accuracy analysis demonstrate that optimal convolutional differentiator methods can solve well the incompatibility between accuracy and efficiency, and are almost twice more accurate than the same-length finite difference. They can efficiently reduce dispersion and provide high-precision waveform data. On the basis of frequency-domain wavefield modeling, we discuss how to directly solve linear equations and point out that when compared to the time-domain methods, frequency-domain methods would be more convenient to handle the multi-source problem and be much easier to incorporate medium attenuation. We also prove the equivalence of the time- and frequency-domain methods by using numerical tests when assumptions with non-relaxation modulus and quality factor are made, and analyze the reason that causes waveform difference. In frequency-domain waveform inversion, experiments have been conducted with transmission, crosshole and reflection data. By using the relation between media scales and characteristic frequencies, we analyze the capacity of the frequency-domain sequential inversion method in anti-noising and dealing with non-uniqueness of nonlinear optimization. In crosshole experiments, we find the main sources of inversion error and figure out how incorrect quality factor would affect inverted results. When dealing with surface reflection data, several frequencies have been chosen with optimal frequency selection strategy, with which we use to carry out sequential and simultaneous inversions to verify how important low frequency data are to the inverted results and the functionality of simultaneous inversion in anti-noising. Finally, I come with some conclusions about the whole work I have done in this dissertation and discuss detailly the existing and would-be problems in it. I also point out the possible directions and theories we should go and deepen, which, to some extent, would provide a helpful reference to researchers who are interested in seismic wave propagation and imaging in complex media.

小波域地震资料的处理方法（WSSP）研究

The modeling formula based on seismic wavelet can well simulate zero - phase wavelet and hybrid-phase wavelet, and approximate maximal - phase and minimal - phase wavelet in a certain sense. The modeling wavelet can be used as wavelet function after suitable modification item added to meet some conditions. On the basis of the modified Morlet wavelet, the derivative wavelet function has been derived. As a basic wavelet, it can be sued for high resolution frequency - division processing and instantaneous feature extraction, in acoordance with the signal expanding characters in time and scale domains by each wavelet structured. Finally, an application example proves the effectiveness and reasonability of the method. Based on the analysis of SVD (Singular Value Decomposition) filter, by taking wavelet as basic wavelet and combining SVD filter and wavelet transform, a new de - noising method, which is Based on multi - dimension and multi-space de - noising method, is proposed. The implementation of this method is discussed the detail. Theoretical analysis and modeling show that the method has strong capacity of de - noising and keeping attributes of effective wave. It is a good tool for de - noising when the S/N ratio is poor. To give prominence to high frequency information of reflection event of important layer and to take account of other frequency information under processing seismic data, it is difficult for deconvolution filter to realize this goal. A filter from Fourier Transform has some problems for realizing the goal. In this paper, a new method is put forward, that is a method of processing seismic data in frequency division from wavelet transform and reconstruction. In ordinary seismic processing methods for resolution improvement, deconvolution operator has poor part characteristics, thus influencing the operator frequency. In wavelet transform, wavelet function has very good part characteristics. Frequency - division data processing in wavelet transform also brings quite good high resolution data, but it needs more time than deconvolution method does. On the basis of frequency - division processing method in wavelet domain, a new technique is put forward, which involves 1) designing filter operators equivalent to deconvolution operator in time and frequency domains in wavelet transform, 2) obtaining derivative wavelet function that is suitable to high - resolution seismic data processing, and 3) processing high resolution seismic data by deconvolution method in time domain. In the method of producing some instantaneous characteristic signals by using Hilbert transform, Hilbert transform is very sensitive to high - frequency random noise. As a result, even though there exist weak high - frequency noises in seismic signals, the obtained instantaneous characteristics of seismic signals may be still submerged by the noises. One method for having instantaneous characteristics of seismic signals in wavelet domain is put forward, which obtains directly the instantaneous characteristics of seismic signals by taking the characteristics of both the real part (real signals, namely seismic signals) and the imaginary part (the Hilbert transfom of real signals) of wavelet transform. The method has the functions of frequency division and noise removal. What is more, the weak wave whose frequency is lower than that of high - frequency random noise is retained in the obtained instantaneous characteristics of seismic signals, and the weak wave may be seen in instantaneous characteristic sections (such as instantaneous frequency, instantaneous phase and instantaneous amplitude). Impedance inversion is one of tools in the description of oil reservoir. one of methods in impedance inversion is Generalized Linear Inversion. This method has higher precision of inversion. But, this method is sensitive to noise of seismic data, so that error results are got. The description of oil reservoir in researching important geological layer, in order to give prominence to geological characteristics of the important layer, not only high frequency impedance to research thin sand layer, but other frequency impedance are needed. It is difficult for some impedance inversion method to realize the goal. Wavelet transform is very good in denoising and processing in frequency division. Therefore, in the paper, a method of impedance inversion is put forward based on wavelet transform, that is impedance inversion in frequency division from wavelet transform and reconstruction. in this paper, based on wavelet transform, methods of time - frequency analysis is given. Fanally, methods above are in application on real oil field - Sansan oil field.

弱各向异性介质中传播的地震波及其速度重建

We begin our studies to make the best of information of seismic data and carry out the description of cracks parameters by extracting anisotropic information. The researching contents are: (1) velocity and polarization anomaly of seismic wave (qP and qSV wave) in weak anisotropic media; (2) reflection seismic synthetic record in anisotropic media; (3) multiple scattering induced by cracks; (4) anisotropic structure inversion and velocity reconstruction with VSP (Vertical Seismic Profile) data; (5) multi-parameters analysis of anisotropy in time-domain and depth-domain. Then we obtain results as follows: (1) We achieve approximate relation of qP and qSV wave's velocity and polarization property in weak anisotropic media. At the same time, we calculate anisotropic velocity factors and polarization anomaly of several typical sedimentary rocks. The results show there are different anisotropic velocity factors and polarization anomaly in different rocks. It is one of the primary theoretical foundation which is expected to identify lithology; (2) We calculate reflection seismic synthetic record with theoretical model; (3) We simulate scattering induced by cracks with Boundary Element Method. Numerical studies show that in the presence of cracks; spatial and scale-length distributions are important and cannot be ignored in modeling cracked solids; (4) From traveltimes information of VSP data, we study the velocity parameter inversion of seismic wave under isotropic and anisotropic models, and its result indicate that the inversion imaging under anisotropic model will not destroy the original features of isotropic model, but it will bring on some bigger error if we adopt the method of isotropic model for anisotropic model data. Further more, basing on the study we develop the CDP mapping technology of reflecting structure under isotropic and anisotropic models, and we process real data as a trial of the methods; (5) We study the problem of initial model reconstruction of anisotropic parameters structure represented by Anderson parameter in depth domain for surface data.

利用分频与反演提高地震资料信噪比与分辨率方法研究

This thesis mainly talks about the wavelet transfrom and the frequency division method. It describes the frequency division processing on prestack or post-stack seismic data and application of inversion noise attenuation, frequency division residual static correction and high resolution data in reservoir inversion. This thesis not only describes the frequency division and inversion in theory, but also proves it by model calculation. All the methods are integrated together. The actual data processing demonstrates the applying results. This thesis analyzes the differences and limitation between t-x prediction filter and f-x prediction filter noise attenuation from wavelet transform theory. It considers that we can do the frequency division attenuation process of noise and signal by wavelet frequency division theory according to the differences of noise and signal in phase, amplitude and frequency. By comparison with the f-x coherence noise, removal method, it approves the effects and practicability of frequency division in coherence and random noise isolation. In order to solve the side effects in non-noise area, we: take the area constraint method and only apply the frequency division processing in the noise area. So it can solve the problem of low frequency loss in non-noise area. The residual moveout differences in seismic data processing have a great effect on stack image and resolutions. Different frequency components have different residual moveout differences. The frequency division residual static correction realizes the frequency division and the calculation of residual correction magnitude. It also solves the problems of different residual correction magnitude in different frequency and protects the high frequency information in data. By actual data processing, we can get good results in phase residual moveout differences elimination of pre-stack data, stack image quality and improvement of data resolution. This thesis analyses the characters of the random noises and its descriptions in time domain and frequency domain. Furthermore it gives the inversion prediction solution methods and realizes the frequency division inversion attenuation of the random noise. By the analysis of results of the actual data processing, we show that the noise removed by inversion has its own advantages. By analyzing parameter's about resolution and technology of high resolution data processing, this thesis describes the relations between frequency domain and resolution, parameters about resolution and methods to increase resolution. It also gives the processing flows of the high resolution data; the effect and influence of reservoir inversion caused by high resolution data. Finally it proves the accuracy and precision of the reservoir inversion results. The research results of this thesis reveal that frequency division noise attenuation, frequency residual correction and inversion noise attenuation are effective methods to increase the SNR and resolution of seismic data.

认知电位时空模式研究

This report mainly focused on methodology of spatiotemporal patterns (STP) of cognitive potentials or event-related potentials (ERP). The representation of STP of brain wave is an important issue in the research of neural assemblies. This paper described methods of parametric 3D head or brain modeling and its corresponding interpolation for functional imaging based on brain waves. The 3D interpolation method is an extension of cortical imaging technique. It can be used with transformed domain features of brain wave on realistic head or brain models. The simulating results suggests that it is a better method in comparison with the global nearest neighbor technique. A stable and definite STP of brainwave referred as microstate may become basic element for comprehending sophisticated cognitive processes. Fuzzy c-mean algorithm was applied to segmentation STPs of ERP into microstates and corresponding membership functions. The optimal microstate number was estimated with both the trends of objective function against increasing clustering number and the decorrelation technique base don microstate shape similarity. Comparable spatial patterns may occur at different moments in time with fuzzy indices and thus the serial processing limit generated from behavioral methods has been break through. High-resolution frequency domain analysis was carried out with multivariate autoregressive model. Bases on a 3D interpolation mentioned above, visualization of dynamical coordination of cerebral network was realized with magnitude-squared partial coherence. Those technique illustrated with multichannel ERP of 9 subjects when they undertook Strop task. Stroop effects involves several regions during post-perception stage with technique of statistical parameter mapping based F-test [SPM(F)]. As SPM(F) suggested task effects occurred within 100 ms after stimuli presentation involved several sensory regions, it may reflect the top-down processing effect.