Multiple quasi-phase-matching for enhanced generation of selected high harmonics in aperiodic modulated fibers

A technique for enhanced generation of selected high harmonics in a gas medium, in a high ionization limit, is proposed in this paper. An aperiodically corrugated hollow-core fiber is employed to modulate the intensity of the fundamental laser pulse along the direction of propagation, resulting in multiple quasi-phase-matched high harmonic emissions at the cutoff region. Simulated annealing (SA) algorithm is applied for optimizing the aperiodic hollow-core fiber. Our simulation shows that the yield of selected harmonics is increased equally by up to 2 orders of magnitude compared with no modulation and this permits flexible control of the quasi-phase-matched emission of selected harmonics by appropriate corrugation. (c) 2007 Optical Society of America.

Far-field spot compression without energy loss in main lob in wireless laser communication

Far-field spot compression without energy loss in main lob is of great significance to wireless laser communication. In this letter, we propose two schemes to obtain far-field spot compression without energy loss in main lob. One scheme is based on the simulated annealing (SA) algorithm. Using SA algorithm, we design the phase profile of the diffractive phase element (DPE). Using the designed DPE, far-field spot compression without energy loss in main lob is achieved. The other scheme is based on YG algorithm. By means of YG algorithm, we appropriately designed the DPE in the emitting plane. Using the DPE, far-field spot compression without energy loss in main lob is obtained. (c) 2007 Elsevier GmbH. All rights reserved.

A metal-mirror-based reflecting polarizing beam splitter

We theoretically investigated the design of a metal-mirror-based reflecting polarizing beam splitter (RPBS). The metal mirror is a silver slab, which is embedded in the substrate of a rectangular silica transmission grating. By using a modal analysis and rigorous coupled-wave analysis, an RPBS grating is designed for operation at 1550 nm. When it is illuminated in Littrow mounting, the transverse electric (TE) and transverse magnetic (TM) waves will be mainly reflected in the minus-first and zeroth orders, respectively. Moreover, a wideband RPBS grating is obtained by adopting the simulated annealing algorithm. The RPBS gratings exhibit high diffraction efficiencies (similar to 95%) and high extinction ratios over a certain angle and wavelength range, especially for the minus-first-order reflection. This kind of RPBS should be useful in practical optical applications.

Nanoelectronic Circuit Architectures Based on Single-Electron Turnstiles

Single-electron devices (SEDs) have ultra-low power dissipation and high integration density, which make them promising candidates as basic circuit elements of the next generation VLSI circuits. In this paper, we propose two novel circuit single-electron architectures: the single-electron simulated annealing algorithm (SAA) circuit and the single-electron cellular neural network (CNN). We used the MOSFET-based single-electron turnstile [1] as the basic circuit element. The SAA circuit consists of the voltage-controlled single-electron random number generator [2] and the single-electron multiple-valued memories (SEMVs) [3]. The random-number generation and variable variations in SAA are easily achieved by transferring electrons using the single-electron turnstile. The CNN circuit used the floating-gate single-electron turnstile as the neural synapses, and the number of electrons is used to represent the cells states. These novel circuits are promising in future nanoscale integrated circuits.

地形可视性分析中观察点设置问题研究

观察点设置问题是地形可视性分析中的一类重要问题，对该问题的研究可以在空间信息辅助决策、通信、旅游、野生动物保护等领域发挥重大作用。本文在对地形可视性分析中观察点设置问题现有研究成果总结和分析基础上对该问题展开深入研究。 首先，针对现有解决方法只从智能算法或地形数据表示方法单一角度进行分析和研究的局限性，提出了一种问题相关的智能算法和数据表示方法相结合的解决问题新框架。该框架考虑了解决观察点设置问题时智能算法的优点和数据表示方式的特点相互配合问题，目的是充分发挥二者各自的优势以提高观察点设置问题解决的准确度与效率。 其次，在深入分析观察点设置问题本身特点的基础上，结合隶属云理论的基本理论和方法，对经典模拟退火算法从退温函数设计、温度产生过程、状态生成过程三方面进行了问题相关的改进，提出了一种适于观察点设置问题的改进模拟退火算法(Improved Simulated Annealing algorithm, ISA)。该算法一方面保持了经典模拟退火算法的稳定倾向特性，保证了算法满足伴随退火温度的不断下降，对恶化的新状态越来越难于接受这一模拟退火算法的最基本特征；另一方面其退火温度的连续性随机变化特性和隐含的“回火升温”过程，则有利于算法有效拒绝恶化解，加速算法收敛，能够更好地满足观察点设置问题对于算法收敛速度的要求。 再次，在分析地形数据的精度、误差等因素对观察点设置问题的解决准确性和解决效率影响程度的基础上，提出了一种基于离散余弦变换的地形数据内插方法(Discrete Cosine Transformation Interpolation method, DCTI)。新方法将传统空域上的地形内插转换到变换域上进行，同时充分利用了离散余弦变换的熵保持特性和能量压缩特性，简化了变换域上的内插过程，提高了地形数据内插的效率和精度。DCTI方法与其他现有典型地形数据内插方法相比，对地形可视性信息获取的准确性和效率影响最小，为平衡观察点设置问题解决过程中时间效率和准确度之间的关系，最终有效地解决观察点设置问题提供了数据基础。 最后，从智能算法和地形数据相结合的角度出发，提出了一种基于ISA和DCTI相结合的观察点设置问题多分辨率处理方法(Multi-Resolution Processing method, MRP)。新方法将模拟退火算法的逐次退火特点和地形数据的多分辨率表示充分结合，达到了发挥算法数据相结合的综合优势的目的。与现有单纯基于模拟退火算法的解决方法相比，在问题解决准确度保持不变的前提下，基于MRP方法的观察点设置问题解决的平均耗时减少85%~95%，为实际工程应用问题的解决提供了一条重要途径。

Three dimensional solution structure of neurotoxin CM-11 from the Ophiophagus hannah

The Ophiophagus hannah (King Cobra) neurotoxin CM-11 is a small protein with 72 amino acid residues, Based on complete assignments of H-1-NMR resonances and determination of secondary structures of CM-11, 349 distance and 27 dihedral angle constraints including 19 phi's and 8 chi's were collected from NOESY and DQF-COSY , and the chemical stereospecific assignment of beta(1)H was partially achieved, Twelve structures with lower energy was obtained via metric matrix distance geometry and refinement with simulated annealing, These structures have a low RMSD of 0.14 nm for backbone atoms and 0.20 nm for heavy atoms, with no distance constraint violation more than 0.05 nm, and no dihedral angle violation more than 3 degrees.

青藏高原及邻区面波频散反演--岩石圈-软流层体系结构研究

The Qinghai-Tibet Plateau lies in the place of the continent-continent collision between Indian and Eurasian plates. Because of their interaction the shallow and deep structures are very complicated. The force system forming the tectonic patterns and driving tectonic movements is effected together by the deep part of the lithosphere and the asthenosphere. It is important to study the 3-D velocity structures, the spheres and layers structures, material properties and states of the lithosphere and the asthenosphere for getting knowledge of their formation and evolution, dynamic process, layers coupling and exchange of material and energy. Based on the Rayleigh wave dispersion theory, we study the 3-D velocity structures, the depths of interfaces and thicknesses of different layers, including the crust, the lithosphere and the asthenosphere, the lithosphere-asthenosphere system in the Qinghai-Tibet Plateau and its adjacent areas. The following tasks include: (1)The digital seismic records of 221 seismic events have been collected, whose magnitudes are larger than 5.0 over the Qinghai-Tibet Plateau and its adjacent areas. These records come from 31 digital seismic stations of GSN , CDSN、NCDSN and part of Indian stations. After making instrument response calibration and filtering, group velocities of fundamental mode of Rayleigh waves are measured using the frequency-time analysis (FTAN) to get the observed dispersions. Furthermore, we strike cluster average for those similar ray paths. Finally, 819 dispersion curves (8-150s) are ready for dispersion inversion. (2)From these dispersion curves, pure dispersion data in 2°×2° cells of the areas (18°N-42°N, 70°E-106°E) are calculated by using function expansion method, proposed by Yanovskaya. The average initial model has been constructed by taking account of global AK135 model along with geodetic, geological, geophysical, receiving function and wide-angle reflection data. Then, initial S-wave velocity structures of the crust and upper mantle in the research areas have been obtained by using linear inversion (SVD) method. (3)Taking the results of the linear inversion as the initial model, we simultaneously invert the S wave velocities and thicknesses by using non-linear inversion (improved Simulated Annealing algorithm). Moreover, during the temperature dropping the variable-scale models are used. Comparing with the linear results, the spheres and layers by the non-linear inversion can be recognized better from the velocity value and offset. (4)The Moho discontinuity and top interface of the asthenosphere are recognized from the velocity value and offset of the layers. The thicknesses of the crust, lithosphere and asthenosphere are gained. These thicknesses are helpful to studying the structural differentia between the Qinghai-Tibet Plateau and its adjacent areas and among geologic units of the plateau. The results of the inversion will provide deep geophysical evidences for studying deep dynamical mechanism and exploring metal mineral resource and oil and gas resources. The following conclusions are reached by the distributions of the S wave velocities and thicknesses of the crust, lithosphere and asthenosphere, combining with previous researches. (1)The crust is very thick in the Qinghai-Tibet Plateau, varying from 60 km to 80 km. The lithospheric thickness in the Qinghai-Tibet Plateau is thinner (130-160 km) than its adjacent areas. Its asthenosphere is relatively thicker, varies from 150 km to 230 km, and the thickest area lies in the western Qiangtang. India located in south of Main Boundary thrust has a thinner crust (32-38 km), a thicker lithosphere of about 190 km and a rather thin asthenosphere of only 60 km. Sichuan and Tarim basins have the crust thickness less than 50km. Their lithospheres are thicker than the Qinghai-Tibet Plateau, and their asthenospheres are thinner. (2)The S-wave velocity variation pattern in the lithosphere-asthenosphere system has band-belted distribution along east-westward. These variations correlate with geology structures sketched by sutures and major faults. These sutures include Main Boundary thrust (MBT), Yarlung-Zangbo River suture (YZS), Bangong Lake-Nujiang suture (BNS), Jinshajiang suture (JSJS), Kunlun edge suture (KL). In the velocity maps of the upper and middle crust, these sutures can be sketched. In velocity maps of 250-300 km depth, MBT, BNS and JSJS can be sketched. In maps of the crustal thickness, the lithospheric thickness and the asthenospheric thickness, these sutures can be still sketched. In particular, MBT can be obviously resolved in these velocity maps and thickness maps. (3)Since the collision between India and Eurasian plate, the “loss” of surface material arising from crustal shortening is caused not only by crustal thickening but also by lateral extrusion material. The source of lateral extrusion lies in the Qiangtang block. These materials extrude along the JSJS and BNS with both rotation and dispersion in Daguaiwan. Finally, it extends toward southeast direction. (4)There is the crust-mantle transition zone of no distinct velocity jump in the lithosphere beneath the Qiangtang Terrane. It has thinner lithosphere and developed thicker asthenosphere. It implies that the crust-mantle transition zone of partial melting is connected with the developed asthenosphere. The underplating of asthenosphere may thin the lithosphere. This buoyancy might be the main mechanism and deep dynamics of the uplift of the Qinghai-Tibet hinterland. At the same time, the transport of hot material with low velocity intrudes into the upper mantle and the lower crust along cracks and faults forming the crust-mantle transition zone.

非均质储层叠前地震反演技术研究及应用

(1) I research on the relationship between elastic parameters, lithology and liquid. It is a physical base for pre-stack seismic inversion. I research all kinds of approximate expressions of Zoeppritz function. Then the relation of all kinds of approximate expressions can be confirmed. The geological model of water sand and gas sand in different depth was designed. Moreover I research on precision of all kinds of approximate expressions. (2) In process of seismic data which aim at amplitude recovery and apply in pre-stack seismic inversion, I advance to adopt double flow chart for different aim. Pre-stack noise elimination, real amplitude recovery and NMO correction of long offset are the key taches. (3) I made a systemic expatiate for the thinking and applicability about all kinds of expressions of elastic impedance. And mathematical model was applied to compare the precision with all kinds of expressions of elastic impedance. I propose a new pre-stack simultaneous inversion which is based on the Zoppritz function and simulated annealing algorithm. This method can ensure calculation precision of reflection coefficient from different incident angle and get a global optimum solution. Therefore this method improves the precision of pre-stack seismic inversion. (4) The object function of P-S wave pre-stack simultaneous inversion was established. I compared the precision and convergence between simultaneous inversion and P-wave inversion. And the results show that simultaneous inversion is superior to P-wave inversion. Through the study of AVO event of transformed wave, AVO characters of different kinds of gas sand were analyzed. (5) I carried out the study work of pre-stack seismic inversion for carbonate reservoir in middle of Tarim basin and sand shale reservoir in Sulige Area of Erdos Basin. The method and technology in this paper was applied to practical work. And I made a prediction for heterogeneous reservoir. Moreover it acquires a good application effect. Key Word: reflection coefficient， amplitude recovery， pre-stack seismic inversion， Heterogeneous reservoir，prediction.

Fourier有限差分法深度偏移与正演模拟

With the development of seismic exploration, the target becomes more and more complex, which leads to a higher demand for the accuracy and efficiency in 3D exploration. Fourier finite-difference (FFD) method is one of the most valuable methods in complex structure exploration, which keeps the ability of finite-differenc method in dealing with laterally varing media and inherits the predominance of the phase-screen method in stablility and efficiency. In this thesis, the accuracy of the FFD operator is highly improved by using simulated annealing algorithm. This method takes the extrapolation step and band width into account, which is more suitable to various band width and discrete scale than the commonely-used optimized method based on velocity contrast alone. In this thesis, the FFD method is extended to viscoacoustic modeling. Based on one-way wave equation, the presented method is implemented in frequency domain; thus, it is more efficient than two-way methods, and is more convenient than time domain methods in handling attenuation and dispersion effects. The proposed method can handle large velocity contrast and has a high efficiency, which is helpful to further research on earth absorption and seismic resolution. Starting from the frequency dispersion of the acoustic VTI wave equation, this thesis extends the FFD migration method to the acoustic VTI media. Compared with the convetional FFD method, the presented method has a similar computational efficiency, and keeps the abilities of dealing with large velocity contrasts and steep dips. The numerical experiments based on the SEG salt model show that the presented method is a practical migration method for complex acoustical VTI media, because it can handle both large velocity contrasts and large anisotropy variations, and its accuracy is relatively high even in strong anisotropic media. In 3D case, the two-way splitting technique of FFD operator causes artificial azimuthal anisotropy. These artifacts become apparent with increasing dip angles and velocity contrasts, which prevent the application of the FFD method in 3D complex media. The current methods proposed to reduce the azimuthal anisotropy significantly increase the computational cost. In this thesis, the alternating-direction-implicit plus interpolation scheme is incorporated into the 3D FFD method to reduce the azimuthal anisotropy. By subtly utilizing the Fourier based scheme of the FFD method, the improved fast algorithm takes approximately no extra computation time. The resulting operator keeps both the accuracy and the efficiency of the FFD method, which is helpful to the inhancements of both the accuracy and the efficiency for prestack depth migration. The general comparison is presented between the FFD operator and the generalized-screen operator, which is valuable to choose the suitable method in practice. The percentage relative error curves and migration impulse responses show that the generalized-screen operator is much sensiutive to the velocity contrasts than the FFD operator. The FFD operator can handle various velocity contrasts, while the generalized-screen operator can only handle some range of the velocity contrasts. Both in large and weak velocity contrasts, the higher order term of the generalized-screen operator has little effect on improving accuracy. The FFD operator is more suitable to large velocity contrasts, while the generalized-screen operator is more suitable to middle velocity contrasts. Both the one-way implicit finite-difference migration and the two-way explicit finite-differenc modeling have been implemented, and then they are compared with the corresponding FFD methods respectively. This work gives a reference to the choosen of proper method. The FFD migration is illustrated to be more attractive in accuracy, efficiency and frequency dispertion than the widely-used implicit finite-difference migration. The FFD modeling can handle relatively coarse grids than the commonly-used explicit finite-differenc modeling, thus it is much faster in 3D modeling, especially for large-scale complex media.

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

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.

砂泥岩薄互层储层预测方法与应用研究

With the deeply development of exploration and development in petroleum in China, new increasing reserves are found in old oil fields and the verge of the old ones through re-study of geological property. It is more and more important to discovery and develop thin layer or thin inter-bedded layers reservoirs. All of the targets are thin sand-shale inter-bedded reservoirs and the core technology is reservoir predictions between wells in thin sand-shale inter-bedded layers. The continuity of the thin sand-shale inter-bedded layers in space or separating and heterogeneity is the key of reservoir geology research. The seismic reflection, high resolution analysis method and inversion method to thin sand-shale inter-bedded layers are thorough discussed and deeply studied in this paper to try to find the methods and resolutions of reservoir geology research. The below is followed. 1. Based on the pre-research of other people, five models are created: the sand sphenoid body, interlay sandstone and interlay shale of the equal thickness, interlay sandstone of the equal thickness and interlay shale of the unequal thickness, interlay sandstone of the unequal thickness and interlay shale of the unequal thickness, interlay sandstone of the changing thickness in sequence and interlay shale of the changing thickness in sequence. Then the study of the forward modeling are conducted on the thin layer and thin inter-bedded layers geological characters and seismic reflections including amplitude, frequency, phase, wave shape and time-frequency responding in the domains of time and frequency. The affect of petro-physics difference of layers, single thin layer thickness, thickness of inter-bedded, layer number of inter-bedded, incident wavelet domain frequency and types, sample interval to seismic reflection characters, frequency spectrum and time-frequency respond of reflectivity is theoretically discussed. 2. Qualitatively analyzing the sedimentary rhythm of the thin inter-bedded layers in vertical orientation and computing the single layer thickness or the average thickness with the method of generalized S transform. Identifying the reflecting interface or lithology interface using the amplitude value of amplitude spectrum domain frequency. 3. Based on the seismic respond of thin sand-shale inter-bedded layers, bring out the high resolution analysis method of seismic data in thin sand-shale inter-bedded layers using wavelet analysis and the idea of affecting low and high frequency with middle frequency. Then analyzing the effect to the method and testing some wavelets in the method. This method is applied to the theoretical models and the field data. 4. Bring forward one improved very fast simulated annealing method (IVFSA) to resolve the problem nonlinearity and multi-parameters of the inversion in thin inter-bedded layers. And IVFSA is more productive and higher precision than general ways. 5. New target constrained function is used in the inversion based on the property of the inversion in thin inter-bedded layers. 6. Making the full use of geological and logging information, IVFSA and the new function are applied in the non-linear inversion to improve reservoir prediction and evaluation in thin inter-bedded formations combined with the idea of logging and seismic inversion. This method was applied to the field data and got good results.

全局优化的Fourier有限差分偏移方法研究

Seismic exploration is the main method of seeking oil and gas. With the development of seismic exploration, the target becomes more and more complex, which leads to a higher demand for the accuracy and efficiency in seismic exploration. Fourier finite-difference (FFD) method is one of the most valuable methods in complex structure exploration, which has obtained good effect. However, in complex media with wider angles, the effect of FFD method is not satisfactory. Based on the FFD operator, we extend the two coefficients to be optimized to four coefficients, then optimize them globally using simulated annealing algorithm. Our optimization method select the solution of one-way wave equation as the objective function. Except the velocity contrast, we consider the effects of both frequency and depth interval. The proposed method can improve the angle of FFD method without additional computation time, which can reach 75° in complex media with large lateral velocity contrasts and wider propagation angles. In this thesis, combinating the FFD method and alternative-direction-implicit plus interpolation(ADIPI） method, we obtain 3D FFD with higher accuracy. On the premise of keeping the efficiency of the FFD method, this method not only removes the azimuthal anisotropy but also optimizes the FFD mehod, which is helpful to 3D seismic exploration. We use the multi-parameter global optimization method to optimize the high order term of FFD method. Using lower-order equation to obtain the approximation effect of higher-order equation, not only decreases the computational cost result from higher-order term, but also obviously improves the accuracy of FFD method. We compare the FFD, SAFFD（multi-parameter simulated annealing globally optimized FFD）, PFFD, phase-shift method(PS）, globally optimized FFD （GOFFD）, and higher-order term optimized FFD method. The theoretical analyses and the impulse responses demonstrate that higher-order term optimized FFD method significantly extends the accurate propagation angle of the FFD method, which is useful to complex media with wider propagation angles.