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%，为实际工程应用问题的解决提供了一条重要途径。

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

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.

全局优化的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.

超参数控制下的波阻抗约束反演

The seismic survey is the most effective geophysical method during exploration and development of oil/gas. As a main means in processing and interpreting seismic data, impedance inversion takes up a special position in seismic survey. This is because the impedance parameter is a ligament which connects seismic data with well-logging and geological information, while it is also essential in predicting reservoir properties and sand-body. In fact, the result of traditional impedance inversion is not ideal. This is because the mathematical inverse problem of impedance is poor-pose so that the inverse result has instability and multi-result, so it is necessary to introduce regularization. Most simple regularizations are presented in existent literature, there is a premise that the image(or model) is globally smooth. In fact, as an actual geological model, it not only has made of smooth region but also be separated by the obvious edge, the edge is very important attribute of geological model. It's difficult to preserve these characteristics of the model and to avoid an edge too smooth to clear. Thereby, in this paper, we propose a impedance inverse method controlled by hyperparameters with edge-preserving regularization, the inverse convergence speed and result would be improved. In order to preserve the edge, the potential function of regularization should satisfy nine conditions such as basic assumptions edge preservation and convergence assumptions etc. Eventually, a model with clear background and edge-abnormity can be acquired. The several potential functions and the corresponding weight functions are presented in this paper. The potential functionφLφHL andφGM can meet the need of inverse precision by calculating the models. For the local constant planar and quadric models, we respectively present the neighborhood system of Markov random field corresponding to the regularization term. We linearity nonlinear regularization by using half-quadratic regularization, it not only preserve the edge, and but also simplify the inversion, and can use some linear methods. We introduced two regularization parameters (or hyperparameters) λ2 and δ in the regularization term. λ2 is used to balance the influence between the data term and the transcendental term; δ is a calibrating parameter used to adjust the gradient value at the discontinuous position(or formation interface). Meanwhile, in the inverse procedure, it is important to select the initial value of hyperparameters and to change hyperparameters, these will then have influence on convergence speed and inverse effect. In this paper, we roughly give the initial value of hyperparameters by using a trend- curve of φ-(λ2, δ) and by a method of calculating the upper limit value of hyperparameters. At one time, we change hyperparameters by using a certain coefficient or Maximum Likelihood method, this can be simultaneously fulfilled with the inverse procedure. Actually, we used the Fast Simulated Annealing algorithm in the inverse procedure. This method overcame restrictions from the local extremum without depending on the initial value, and got a global optimal result. Meanwhile, we expound in detail the convergence condition of FSA, the metropolis receiving probability form Metropolis-Hasting, the thermal procession based on the Gibbs sample and other methods integrated with FSA. These content can help us to understand and improve FSA. Through calculating in the theoretic model and applying it to the field data, it is proved that the impedance inverse method in this paper has the advantage of high precision practicability and obvious effect.

Application of the Conjugate Gradient Optimization Algorithm of Neural Network on Laser Machining

在应用激光技术加工复杂曲面时,通常以采样点集为插值点来建立曲面函数,然后实现曲面上任意坐标点的精确定位。人工神经网络的BP算法能实现函数插值,但计算精度偏低,往往达不到插值精确要求,造成较大的加工误差。提出人工神经网络的共轭梯度最优化插值新算法,并通过实例仿真,证明了这种曲面精确定位方法的可行性,从而为激光加工的三维精确定位提供了一种良好解决方案。这种方法已经应用在实际中。

Simulating shock to detonation transition: Algorithm and results

An algorithm based on flux-corrected transport and the Lagrangian finite element method is presented for solving the problem of shock dynamics. It is verified through the model problem of one-dimensional strain elastoplastic shock wave propagation that the algorithm leads to stable, non-oscillatory results. Shock initiation and detonation wave propagation is simulated using the algorithm, and some interesting results are obtained. (C) 1999 Academic Press.

Effect of Cold Rolling Ratio and Heating Rate of Batch Annealing on Textures of Micro-Carbon Steel Sheet

在四辊冷轧试验机和Gleeble-1500试验机上进行了热轧微碳钢板的冷轧和退火试验。用D/max-RC衍射仪测量了试样的,/”层织构,并用Roe软件进行了ODF分析。研究表明,所研究的热轧微碳深冲板压下率约为75%,退火升温速度为20-40℃/h时,试样为{111}织构特征；压下率较大(80%)时,退火织构为较弱的{111}组分。无论{111}织构还是非{111}织构都是在形核阶段开始形成,在晶粒长大优先长大,受到定向形核和选择生长双重机制的作用。