Second-order Stokes solutions for internal waves in three-layer density-stratified fluid

In this paper, internal waves in three-layer stratified fluid are investigated by using a perturbation method, and the second-order asymptotic solutions of the velocity potentials and the second-order Stokes solutions of the associated elevations of the interfacial waves are presented based on the small amplitude wave theory. As expected, the first-order solutions are consistent with ordinary linear theoretical results, and the second-order solutions describe the second-order modification on the linear theory and the interactions between the two interfacial waves. Both the first-order and second-order solutions derived depend on the depths and densities of the three-layer fluid. It is also noted that the solutions obtained from the present work include the theoretical results derived by Umeyama as special cases.

Second-order Stokes wave solutions for intefacial waves in three-layer stratified fluid with background current

In this paper, interfacial waves in three-layer stratified fluid with background current are investigated using a perturbation method, and the second-order asymptotic solutions of the velocity potentials and the second-order Stokes wave solutions of the associated elevations of the interfacial waves are presented based on the small amplitude wave theory, and the Kelvin-Helmholtz instability of interfacial waves is studied. As expected, for three-layer stratified fluid with background current, the first-order asymptotic solutions (linear wave solutions), dispersion relation and the second-order asymptotic solutions derived depend on not only the depths and densities of the three-layer fluid but also the background current of the fluids, and the second-order Stokes wave solutions of the associated elevations of the interfacial waves describe not only the second-order nonlinear wave-wave interactions between the interfacial waves but also the second-order nonlinear interactions between the interfacial waves and currents. It is also noted that the solutions obtained from the present work include the theoretical results derived by Chen et al (2005) as a special case. It also shows that with the given wave number k (real number) the interfacial waves may show Kelvin-Helmholtz instability.

A kind of extended Korteweg-de Vries equation and solitary wave solutions for interfacial waves in a two-fluid system

This paper considers interfacial waves propagating along the interface between a two-dimensional two-fluid with a flat bottom and a rigid upper boundary. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. It just focuses on the weakly non-linear small amplitude waves by introducing two small independent parameters: the nonlinearity ratio epsilon, represented by the ratio of amplitude to depth, and the dispersion ratio mu, represented by the square of the ratio of depth to wave length, which quantify the relative importance of nonlinearity and dispersion. It derives an extended KdV equation of the interfacial waves using the method adopted by Dullin et al in the study of the surface waves when considering the order up to O(mu(2)). As expected, the equation derived from the present work includes, as special cases, those obtained by Dullin et al for surface waves when the surface tension is neglected. The equation derived using an alternative method here is the same as the equation presented by Choi and Camassa. Also it solves the equation by borrowing the method presented by Marchant used for surface waves, and obtains its asymptotic solitary wave solutions when the weakly nonlinear and weakly dispersive terms are balanced in the extended KdV equation.

基于解耦控制的SCARA机器人轨迹跟踪控制器设计

SCARA型机器人的控制问题由于其动力学模型中没有重力矩项的作用而得以简化,由于在实际应用中经常要求其高速运动,则对具有强耦合的哥氏力与向心力的控制就成为制约其系统性能的重要问题。提出通过线性变换对机器人系统解耦,将高阶系统转化为解耦的低阶系统进行控制的方法,并且应用极点配置对解耦的系统求解机器人控制器。该方法无需测量关节速度和加速度,只需要测量关节位置信号。所提出的控制器既能保证闭环系统全局渐进稳定,又能通过对线性化系统闭环极点的配置来获得期望的闭环系统响应性能。仿真实验证明了所提出的控制器设计方法的可行性。

网络控制系统稳定性的图理论

针对具有有界时延和数据包丢失的网络控制系统,提出了一种新的稳定性判据.基于Lyapunov方法和图论理论,给出非线性离散和连续网络控制系统渐近稳定的充分条件,获得保持这两类系统稳定的最大允许时延界,得到控制器设计方法.并且,利用区间矩阵的谱特征,给出网络控制系统区间稳定的充分条件.设计算法,获得比例积分反馈控制器增益.算例表明所提方法的有效性。

高阶广义屏波动方程叠前深度偏移

Attaining sufficient accuracy and efficiency of generalized screen propagator and improving the quality of input gathers are often problems　of wave equation presack depth migration,　in this paper,a high order formula of generalized screen propagator　for one-way wave equation is proposed by using the asymptotic expansion of single-square-root operator. Based on the formula，a new generalized screen propagator is developed ,which is composed of split-step Fourier propagator and high order correction terms，the new generalized screen propagator not only improving calculation precision without sharply increasing the quantity of computation，facilitates the suitability of generalized screen propagator to the media with strong lateral velocity variation.　As wave-equation prestack depth migration is sensitive to the quality of input gathers, which greatly affect the output，and the available seismic data processing system has inability to obtain traveltimes corresponding to the multiple arrivals, to estimate of great residual statics， to merge seismic datum from different projects and to design inverse Q filter, we establish　difference equations with an embodiment of Huygens’s principle for obtaining traveltimes corresponding to the multiple arrivals,bring forward a time variable matching filter for seismic datum merging by using the fast algorithm called Mallat tree for wavelet transformations， put forward a method for estimation of residual statics by applying the optimum model parameters estimated by iterative inversion with three organized algorithm,i.e,the CMP intertrace cross-correlation algorithm,the Laplacian image edge extraction algorithm,and the DFP algorithm, and present phase-shift inverse Q filter based on Futterman’s amplitude and phase-velocity dispersion formula and wave field extrapolation theory. All of their numerical and real data calculating results shows that our theory and method are practical and efficient. Key words: prestack depth migration, generalized screen propagator, residual statics，inverse Q filter ,traveltime，3D seismic datum mergence

黄土-红粘土化学风化与晚新生代气候变化——甘肃平凉白水剖面为例

Late Cenozoic has witnessed a series of climate-environmental change which ends with a transitional shift from greenhouse to icehouse conditions. In last two decades, scientists began to employ the tectonic uplift and its weathering effect to interpret the climatic changes during the late Cenozoic. However, this endeavor has partly been restricted by the lacking of regional and global chemical weathering data. The loess-red clay deposit in the Loess Plateau may record the weathering features of the detritus material from the wide range upwind of the Loess Plateau. Therefore geochemistry of the loess-plaeosol and red clay sequences may provide insights into the regional chemical weathering regime and the connection between the chemical weathering and the late Cenozoic climate-environmental change Here we selected 319 samples from the Baishui section near the Pingliang City, Gansu Province, and analyzed them with X-ray fluorescence. Based on the result, we reconstruct the chemical weathering history of the Baishui section since 6Ma. We chose CIA as the proxy for chemical weathering intensity. The CIA ratio in soil units is higher than in adjacent loess horizons, but lower than in the red clay, in good agreement with the field observation. The CIA ratios of the Baishui section correlates well with the global ice volume fluctuations, indicating that the global cooling may contribute a lot to the chemical weathering variations in Chinese Loess Plateau. There are at least 3 million-year time scale variations that can be identified in the chemical weathering intensity curve, i.e., between 3.3 to 2.1 Ma, 1.7 to 0.9 Ma and from 0.9Ma. We think these may reflect the combined effect of the tectonic uplift and ice sheets on monsoon intensity. Other time scale variations can be also observed. In the period between 2.4 and 0.8 Ma, the CIA record display the 400,000 years cycle, which may be resulted from the Tibetan uplift during the Pliocene-early Pleistocene which have significantly amplified the monsoon response sensitivity to the orbital-scale variations in insolation. From 1.2 Ma, the 100,000 years period became intensifying, and particularly after 0.8 Ma, the earlier monsoon response at 400,000 year periodicity was overwhelmed by the ice sheet forcing at 100,000 year periodicity. These may indicate that the expansions of the Northern Hemisphere ice sheets may have crossed a threshold, which enforce the monsoon responding at the 100,000 year periodicity.

波前重建法射线追踪及走时表压缩

As a fast and effective method for approximate calculation of seismic numerical simulation, ray tracing method, which has important theory and practical application value, in terms of seismic theory and seismic simulation, inversion, migration, imaging, simplified from seismic theory according to geometric seismic, means that the main energy of seismic wave field propagates along ray paths in condition of high-frequency asymptotic approximation. Calculation of ray paths and traveltimes is one of key steps in seismic simulation, inversion, migration, and imaging. Integrated triangular grids layout on wavefront with wavefront reconstruction ray tracing method, the thesis puts forward wavefront reconstruction ray tracing method based on triangular grids layout on wavefront, achieves accurate and fast calculation of ray paths and traveltimes. This method has stable and reasonable ray distribution, and overcomes problems caused by shadows in conventional ray tracing methods. The application of triangular grids layout on wavefront, keeps all the triangular grids stable, and makes the division of grids and interpolation of a new ray convenient. This technology reduces grids and memory, and then improves calculation efficiency. It enhances calculation accuracy by accurate and effective description and division on wavefront. Ray tracing traveltime table, which shares the character of 2-D or 3-D scatter data, has great amount of data points in process of seismic simulation, inversion, migration, and imaging. Therefore the traveltime table file will be frequently read, and the calculation efficiency is very low. Due to these reasons, reasonable traveltime table compression will be very necessary. This thesis proposes surface fitting and scattered data compression with B-spline function method, applies to 2-D and 3-D traveltime table compression. In order to compress 2-D (3-D) traveltime table, first we need construct a smallest rectangular (cuboidal) region with regular grids to cover all the traveltime data points, through the coordinate range of them in 2-D surface (3-D space). Then the value of finite regular grids, which are stored in memory, can be calculated using least square method. The traveltime table can be decompressed when necessary, according to liner interpolation method of 2-D (3-D) B-spline function. In the above calculation, the coefficient matrix is stored using sparse method and the liner system equations are solved using LU decomposition based on the multi-frontal method according to the sparse character of the least square method matrix. This method is practiced successfully in several models, and the cubic B-spline function can be the best basal function for surface fitting. It make the construction surface smooth, has stable and effective compression with high approximate accuracy using regular grids. In this way, through constructing reasonable regular grids to insure the calculation efficiency and accuracy of compression and surface fitting, we achieved the aim of traveltime table compression. This greatly improves calculation efficiency in process of seismic simulation, inversion, migration, and imaging.

高斯束射线追踪

Using the approximate high-frequency asymptotic methods to solve the scalar wave equation, we can get the eikonal equation and transport equation. Solving the eikonal equation by the method of characteristics provides a mathematical derivation of ray tracing equations. So, the ray tracing system is folly based on the approximate high-frequency asymptotic methods. If the eikonal is complex, more strictly, the eikonal is real value at the rays and complex outside rays, we can derive the Gaussian beam. This article mainly concentrates on the theory of Gaussian beam. To classical ray tracing theory, the Gaussina beam method (GBM) has many advantages. First, rays are no longer required to stop at the exact position of the receivers; thus time-consuming two-point ray tracing can be avoided. Second, the GBM yields stable results in regions of the wavefield where the standard ray theory fails (e.g., caustics, shadows zones and critical distance). Third, unlike seismograms computed by conventional ray tracing techniques, the GBM synthetic data are less influenced by minor details in the model representation. Here, I realize kinematical and dynamical system, and based on this, realize the GBM. Also, I give some mathematical examples. From these examples, we can find the importance and feasibility of the ray tracing system. Besides, I've studied about the reflection coefficient of inhomogeneous S-electromagnetic wave at the interface of conductive media. Basing on the difference of directions of phase shift constant and attenuation constant when the electromagnetic wave propagates in conductive medium, and using the boundary conditions of electromagnetic wave at the interface of conductive media, we derive the reflection coefficient of inhomogeneous S-electromagnetic wave, and draw the curves of it. The curves show that the quasi total reflection will occur when the electromagnetic wave incident from the medium with greater conductivity to the medium with smaller conductivity. There are two peak, values at the points of the critical angles of phase shift constant and attenuation constant, and the reflection coefficient is smaller than 1. This conclusion is different from that of total reflection light obviously.

高斯束偏移的方法及应用

Gaussian beam is the asymptotic solution of wave equation concentred at the central ray. The Gaussian beam ray tracing method has many advantages over ray tracing method. Because of the prevalence of multipath and caustics in complex media, Kirchhoff migration usually can not get satisfactory images, but Gaussian beam migration can get better results.The Runge-Kutta method is used to carry out the raytracing, and the wavefront construction method is used to calculate the multipath wavefield. In this thesis, a new method to determine the starting point and initial direction of a new ray is proposed take advantage of the radius of curvature calculated by dynamic ray tracing method.The propagation characters of Gaussian beam in complex media are investigated. When Gaussian beam is used to calculate the Green function, the wave field near the source was decomposed in Gaussian beam in different direction, then the wave field at a point is the superposition of individual Gaussian beams.Migration aperture is the key factor for Kirchhoff migration. In this thesis, the criterion for the choice of optimum aperture is discussed taking advantage of stationary phase analysis. Two equivalent methods are proposed, but the second is more preferable.Gaussian beam migration based on dip scanning and its procedure are developed. Take advantage of the travel time, amplitude, and takeoff angle calculated by Gaussian beam method, the migration is accomplished.Using the proposed migration method, I carry out the numerical calculation of simple theoretical model, Marmousi model and field data, and compare the results with that of Kirchhoff migration. The comparison shows that the new Gaussian beam migration method can get a better result over Kirchhoff migration, with fewer migration noise and clearer image at complex structures.

序列视觉信息整合的特性与机制

One of the great puzzles in the psychology of visual perception is that the visual world appears to be a coherent whole despite our viewing it through temporally discontinuous series of eye fixations. The investigators attempted to explain this puzzle from the perspective of sequential visual information integration. In recent years, investigators hypothesized that information maintained in the visual short-term memory (VSTM) could become visual mental images gradually during time delay in visual buffer and integrated with information perceived currently. Some elementary studies had been carried out to investigate the integration between VSTM and visual percepts, but further research is required to account for several questions on the spatial-temporal characteristics, information representation and mechanism of integrating sequential visual information. Based on the theory of similarity between visual mental image and visual perception, this research (including three studies) employed the temporal integration paradigm and empty cell localization task to further explore the spatial-temporal characteristics, information representation and mechanism of integrating sequential visual information (sequential arrays). The purpose of study 1 was to further explore the temporal characteristics of sequential visual information integration by examining the effects of encoding time of sequential stimuli on the integration of sequential visual information. The purpose of study 2 was to further explore the spatial characteristics of sequential visual information integration by investigating the effects of spatial characteristics change on the integration of sequential visual information. The purpose of study 3 was to explore the information representation of information maintained in the VSTM and integration mechanism in the process of integrating sequential visual information by employing the behavioral experiments and eye tracking technology. The results indicated that: (1) Sequential arrays could be integrated without strategic instruction. Increasing the duration of the first array could cause improvement in performance and increasing the duration of the second array could not improve the performance. Temporal correlation model was not fit to explain the sequential array integration under long-ISI conditions. (2) Stimuli complexity influenced not only the overall performance of sequential arrays but also the values of ISI at asymptotic level of performance. Sequential arrays still could be integrated when the spatial characteristics of sequential arrays changed. During ISI, constructing and manipulating of visual mental image of array 1 were two separate processing phases. (3) During integrating sequential arrays, people represented the pattern constituted by the objects' image maintained in the VSTM and the topological characteristics of the objects' image had some impact on fixation location. The image-perception integration hypothesis was supported when the number of dots in array 1 was less than empty cells, and the convert-and-compare hypothesis was supported when the number of the dot in array 1 was equal to or more than empty cells. These findings not only contribute to make people understand the process of sequential visual information integration better, but also have significant practical application in the design of visual interface.