Modeling vertical heat transport in the wave affected surface layer of the ocean

In considering the vertical heat transport problems in the upper ocean, the flat upper boundary approximation for the free surface and the horizontal homogenous hypothesis are usually applied. However, due to the existence of the wave motion, the application of this approximation may result in some errors to the solar irradiation since it decays quickly in respect to the actual thickness of the water layer below the surface; on the other hand, due to the fluctuation of the water layer depth, it is improper to neglect the effects of the horizontal advection and turbulent diffusion since they also contribute to the vertical heat transport. A new model is constructed in this study to reflect these effects. The corresponding numerical simulations show that the wave motion may remarkably accelerate the vertical heat transferring process and the variation of the temperature in the wave affected layer appears in an oscillating manner.

THE EFFECTS OF MONSOONS AND CONNECTIVITY OF SOUTH CHINA SEA ON THE SEASONAL VARIATIONS OF WATER EXCHANGE IN THE LUZON STRAIT

Seasonal variations of water exchange in the Luzon Strait are studied numerically using the improved Princeton Ocean Model (POM) with a consideration of the effects of connectivity of South China Sea (SCS) and monsoons. The numerical simulations are carried out with the strategy of variable grids, coarse grids for the Pacific basin and fine grids for the SCS. It. is shown that the Mindoro Strait plays an important role in adjusting the water balance between the Pacific and the SCS. The SCS monsoon in summer seasons hinders the entrance of the Pacific water into the SCS through the Luzon Strait while the SCS monsoon in winter seasons promotes the entrance of Pacific water into the SCS through the Luzon Strait. However, the SCS monsoon does not affect the annual mean Luzon Strait transport, as is mainly determined by the Pacific basin wind.

Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor-Yamada level 2.5 turbulence closure scheme, was implemented to investigate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic energy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corresponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the temperature gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.

On the nonlinear integral transform of an ocean wave spectrum into an along-track interferometric synthetic aperture radar image spectrum

We present a new nonlinear integral transform relating the ocean wave spectrum to the along-track interferometric synthetic aperture radar (AT-INSAR) image spectrum. The AT-INSAR, which is a synthetic aperture radar (SAR) employing two antennas displaced along the platform's flight direction, is considered to be a better instrument for imaging ocean waves than the SAR. This is because the AT-INSAR yields the phase spectrum and not only the amplitude spectrum as with the conventional SAR. While the SAR and AT-INSAR amplitude spectra depend strongly on the modulation of the normalized radar cross section (NRCS) by the long ocean waves, which is poorly known, the phase spectrum depends only weakly on this modulation. By measuring the phase difference between the signals received by both antennas, AT-INSAR measures the radial component of the orbital velocity associated with the ocean waves, which is related to the ocean wave height field by a well-known transfer function. The nonlinear integral transform derived in this paper differs from the one previously derived by Bao et al. [1999] by an additional term containing the derivative of the radial component of the orbital velocity associated with the long ocean waves. By carrying out numerical simulations, we show that, in general, this additional term cannot be neglected. Furthermore, we present two new quasi-linear approximations to the nonlinear integral transform relating the ocean wave spectrum to the AT-INSAR phase spectrum.

卫星高度计资料在海洋潮汐研究中的应用

本文依据收集到的392个地面验潮站8个主要分潮（M2、S2、K1、O1、N2、K2、P1及Q1）的调和常数，对现有7个全球大洋潮汐模式的准确度进行了检验，结果显示各模式在深海区域均达到了比较高的准确度，相互之间差别也不大。经验模式GOT00和CSR4.0、同化模式NAO99、反演同化模式TPXO7.0、数值同化模式FES2002和FES2004的M2分潮均方根偏差在3 cm左右，其它分潮（S2、K1、O1、N2、K2、P1及Q1）大约在1～2 cm。本文还依据中国近海18个岛屿的调和常数对其中的5个大洋潮汐模式的准确度进行了检验，结果表明，M2分潮均方根偏差在6～14 cm，明显高于大洋部分的偏差，其中日本国家天文台的潮汐模式NAO99在中国近海的结果相对较准确。 我们利用1992年8月至2008年8月的TOPEX/POSEIDON和JASON-1(T/P-J)卫星高度计资料，对沿卫星轨道的302816个站点进行了14个分潮的潮汐调和分析，得到了全球大洋潮汐的8个主要分潮以及2个气象分潮Sa、Ssa的经验同潮图。主要结果有：（1）各分潮在卫星上升轨道与下降轨道的交叉点(约7000个)相关性分析表明：M2分潮的振幅和迟角的相关系数很高（分别为0.9965和0.9961）；S2，K1，O1和Sa分潮也有较好的相关性（相关系数为0.94～0.99）；（2）该结果与392地面个验潮站吻合较好，其中M2分潮的振幅、迟角和向量的均方根偏差分别为：1.73 cm，2.340和2.93 cm；S2，K1和O1分潮的振幅、迟角和向量的均方根偏差为1 cm左右，5.250～7.270和1.5～2.1 cm，该精度与最近几年国际上的主要大洋潮汐模式的准确度相近；（3）首次通过卫星资料获得了Sa、Ssa分潮的同潮图。周期为1年的Sa分潮与大洋105个地面站相比，振幅、迟角和向量的均方根偏差分别为1.50 cm、18.360和2.16 cm。在此基础上，进一步分析了构成Sa、Ssa气象分潮的两个主要因素（海水密度以及海面气压）在全球的分布。 在T/P-J等卫星资料无法覆盖到南大洋和北冰洋，本文利用Princeton Ocean Model（POM）进行了数值模拟，模拟结果与162个地面实测站（其中南大洋30个，北冰洋132个）的观测比较一致。基于卫星资料分析的结果和数值模拟结果合并得到了全球大洋的8个主要分潮同潮图。在此基础上通过全球潮汐能量耗散的计算得到潮能通量的分布，并得到全球M2、S2、K1和O1分潮的潮汐能量耗散率为2.431TW、0.401TW、0.336TW和0.176TW。 本文还利用卫星资料对南海潮汐进行了研究，在中国南海，获得了主要的半日潮、全日潮、四分日分潮和长周期分潮（M2，S2，N2，K2，K1，O1，P1，Q1，M4, MS4，Sa, Ssa）的经验同潮图。与南海沿岸94个地面验潮站的数据符合得比较好，M2，S2，K1及O1等4个主要分潮的平均振幅差为2～4 cm，均方根偏差分别是9～11 cm.其它4个主要分潮N2，K2，P1，Q1的平均振幅差为1～2 cm，均方根偏差为2～4 cm。此外，本文还利用卫星高度计资料潮汐分析结果沿卫星轨道进行高通滤波，分离得出中国近海的M2，S2，K1及O1分潮的内潮信息。

模型直升机航向自适应保性能控制

提出了一种新的模型直升机航向控制算法。针对具有模型不确定性的直升机航向线性模型,提出了一种具有自适应机制的最优保性能控制器。该控制策略通过引入自适应机制降低固定增益控制器所固有的保守性,并且控制器的反馈增益应用线性矩阵不等式(LMIs)方法解得。理论分析和数字仿真表明所设计的控制器具有良好的鲁棒稳定性能。

一种改进的遗传算法

针对传统遗传算法在编码方案及交叉操作中存在的局限性，提出了一种新的遗传算法的改进 方法．该方法（１）以实数编码代替二进制编码，有效地解决了传统遗传算法中二进制编码串的长度与 计算精度、运算量之间的矛盾，（２）根据适应度值对父染色体进行重组操作，克服了传统遗传算法中 交叉操作所存在的盲目性．最后，以求解自然对数和神经网络的训练为例验证了所提出方法的有效 性．

不同基板预热温度对激光金属沉积成形过程热应力影响的研究

研究不同基板预热温度对激光金属沉积成形过程热应力的影响,对于降低成形过程的热应力,抑制成形过程裂缝的产生,减小成形过程试样和基板的翘曲变形具有非常重要的意义。根据有限元分析中的"单元生死"技术,编程建立了基板预热对激光金属沉积成形过程热应力影响的三维多道多层数值模拟模型,详细分析了基板未预热和分别预热到200℃、300℃、400℃、500℃、600℃时对沉积成形过程VonMise’s热应力、X方向、Y方向以及Z方向热应力的影响。在与模拟过程相同的参数下,利用镍基合金粉末分别在基板未预热和分别预热到300℃、400℃、500℃、600℃时进行了成形试验,试验的结果跟数值模拟结果吻合较好。

三峡库区巴东组泥灰岩溶蚀作用与岩体力学性质弱化研究

Abstract The karsrt erosion engineering geology became a highlight problem in recent years, in particularly, the karst erosion of marlite of Badong formation made the rock mechanics weaken in Three Gorges Reservoir area, which reduces the safety of slope. During the immigrant construction, many high slopes have been formed, whose instabilities problems pose serious threats to the safety of the people and properties. The accidents of the slope failure take place now and then. By testing, it has been found that the karst erosion pattern and dissolution rate of marlite are not weaker than that of the pure limestone. Furthermore, owning to the weathering and unloading, the karst erosion of the marlite will reach certain depth of the slope, which is named infiltrated karst erosion. The karst erosion made the rock mass quality of slope or foundation worse in a large scale. The karst erosion geological disasters, taken place or not, has become the main restrictive factors to the social stability and economic development. Thus the karst erosion process and mechanism of marlite of Badong formation are studied as the main content of this dissertation. The weakening characteristic of rock mass mechanics parameters are studied along with the rock mass structure deformation and failure processes in the course of the karst erosion. At first, the conditions and influencing factors of the karst erosion are analyzed in the investigative region, on the basis of different karst erosion phenomenon of the marlite and different failure modes of slope. Then via indoor the karst erosion tests, it is analyzed that the karst erosion will change the rock mass composition and its structure. Through test, the different karst erosion phenomena between micro and macro have been observed, and the karst erosion mechanism of the marlite has been summarized. Damage theory is introduced to explain the feature of dissolution pore and the law of crack propagation in the marlite. By microscope and the references data, it can be concluded that the karst erosion process can be divided into rock minerals damage and rock structural damage. And the percent of karst erosion volume is named damage factor, which can be used to describe the quantify karst erosion degree of marlite. Through test, the rock mechanical properties in the different period of karst erosion are studied. Based on the damage mechanics theory and the test result, the relation between the karst erosion degree of marlite and weakening degree of mechanical properties is summarized. By numerical simulations, the karst erosive rock mass mechanics is verified. The conclusion is drawn as below: to the rock mass of marlite, the karst erosion damage made mechanics parameters variation, the deformation modulus, cohesion, and inter friction angle reduce as the negative exponent with the increasing of the karst erosion volume, however, the Poisson ratio increases as the positive exponent with the karst erosion volume increasing. It should be noticed that the deduced formulations are limited to the test data and certain conditions. It is suitable to the rock mass parametric weakening process after the karst erosion of marlite in Three Gorges Reservoir area. Based on the failure types of marlite slope in the field, the karst erosion and weathering process of rock mass are analyzed. And the evolution law of deformation and failure of the marlite mass is studied. The main failure feature of the marlite slope is the karst erosive structure subsidence mode in Three Gorges Reservoir area. The karst erosive structure subsidence mode is explained as follows: the rock mass undergoes the synthetic influence, such as weathering, unloading, corrosion, and so on, many pores and cavities have been formed in the rock mass interior, the rock mass quality is worsen and the rock mass structure is changed, and then the inherent structure of rock mass is collapsed under its gravity, therefore, the failure mode of compaction and subsidence take place. Finally, two examples are used to verify the rock mass parameters in Three Gorges Reservoir area, and the relationship between the marlite slope stability and the time of karst erosion is proposed.

旋转轴对称介质中的qP波一阶射线追踪

Elastic anisotropy is a very common phenomenon in the Earth’s interior, especial for sedimentary rock as important gas and oil reservoirs. But in the processing and interpretation of seismic data, it is assumption that the media in the Earth’s interior is completely elastic and isotropic, and then the methods based on isotropy are used to deal with anisotropic seismic data, so it makes the seismic resolution lower and the error on images is caused. The research on seismic wave simulation technology can improve our understanding on the rules of seismic wave propagation in anisotropic media, and it can help us to resolve problems caused by anisotropy of media in the processing and interpretation of seismic data. So researching on weakly anisotropic media with rotated axis of symmetry, we study systematically the rules of seismic wave propagation in this kind of media, simulate the process with numerical calculation, and get the better research results. The first-order ray tracing (FORT) formulas of qP wave derived can adapt to every anisotropic media with arbitrary symmetry. The equations are considerably simpler than the exact ray tracing equations. The equations allow qP waves to be treated independently from qS waves, just as in isotropic media. They simplify considerably in media with higher symmetry anisotropy. In isotropic media, they reduce to the exact ray tracing equations. In contrast to other perturbation techniques used to trace rays in weakly anisotropic media, our approach does not require calculation of reference rays in a reference isotropic medium. The FORT-method rays are obtained directly. They are computationally more effective than standard ray tracing equations. Moreover the second-order travel time corrections formula derived can be used to reduce effectively the travel time error, and improve the accuracy of travel time calculation. The tensor transformation equations of weak-anisotropy parameters in media with rotated axis of symmetry derived from the Bond transformation equations resolve effectively the problems of coordinate transformation caused by the difference between global system of coordinate and local system of coordinate. The calculated weak-anisotropy parameters are completely suitable to the first-order ray tracing used in this paper, and their forms are simpler than those from the Bond transformation. In the numerical simulation on ray tracing, we use the travel time table calculation method that the locations of the grids in the ray beam are determined, then the travel times of the grids are obtained by the reversed distance interpolation. We get better calculation efficiency and accuracy by this method. Finally we verify the validity and adaptability of this method used in this paper with numerical simulations for the rotated TI model with anisotropy of about 8% and the rotated ORTHO model with anisotropy of about 20%. The results indicate that this method has better accuracy for both media with different types and different anisotropic strength. Keywords: weak-anisotropy, numerical simulation, ray tracing equation, travel time, inhomogeneity

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

边坡动力响应分析及应用研究

As a marginal subject, dynamic responses of slopes is not only an important problem of engineering geology (Geotechnical problem), but also of other subjects such as seismology, geophysics, seismic engineering and engineering seismic and so on. Owning to the gulf between different subjects, it is arduous to study dynamic responses of slopes and the study is far from ripeness. Studying on the dynamic responses of slopes is very important in theories as well as practices. Supported by hundreds of bibliographies, this paper systemically details the development process of this subject, introduces main means to analyze this subject, and then gives brief remarks to each means respectively. Engineering geology qualitative analysis is the base of slopes dynamic responses study. Because of complexity of geological conditions, engineering geology qualitative analysis is very important in slopes stability study, especially to rock slopes with complex engineering geology conditions. Based on research fruits of forerunners, this paper summarizes factors influencing slopes dynamic stability into five aspects as geology background, stratums, rock mass structure, and topography as well as hydrogeology condition. Based on rock mass structure controlling theory, engineering geology model of the slope is grouped into two typical classes, one is model with obvious controlling discontinuities, which includes horizontal bedded slope, bedding slope, anti-dip slope, slide as well as slope with base rock and weathered crust; the other is model without obvious controlling discontinuities, which includes homogeneous soil slope, joint rock mass slope. Study on slope failure mechanism under dynamic force, the paper concludes that there are two effects will appear in slope during strong earthquake, one is earthquake inertia force, the other is ultra pore pressure buildup. The two effects lead to failure of the slope. To different types of slope failure, the intensity of two effects acting on the slope is different too. To plastic flow failure, pore pressure buildup is dominant; to falling rock failure and toppling failure, earthquake inertia force is dominant in general. This paper briefly introduces the principle of Lagrangian element method. Through a lot of numerical simulations with FLAC3D, the paper comprehensively studies dynamic responses of slopes, and finds that: if the slope is low, displacement, velocity and acceleration are linear enlarging with elevation increasing in vertical direction; if the slope is high enough, displacement, velocity and acceleration are not linear with elevation any more, on the other hand, they fluctuate with certain rhythm. At the same time, the rhythm appears in the horizontal direction in the certain area near surface of the slope. The distribution form of isoline of displacement, velocity and acceleration in the section of the slope is remarkably affected by the slope angle. In the certain area near the slope surface, isoline of displacement, velocity and acceleration is parallel to the surface of the slope, in the mean time, the strike direction of the extreraum area is parallel to the surface of the slope too. Beyond this area, the isoline direction and the strike direction of the extremum area turn to horizontal with invariable distance. But the rhythm appearing or not has nothing to with the slope angle. The paper defines the high slope effect and the low slope effect of slopes dynamic responses, discusses the threshold height H^t of the dynamic high slope effect, and finds that AW is proportional to square root of the dynamic elastic moduli El P , at the same time, it is proportional to period Tof the dynamic input. Thus, the discriminant of H^t is achieved. The discriminant can tell us that to a slope, if its height is larger than one fifth of the wavelength, its response regular will be the dynamic high slope effect; on the other hand, its response regular will be the dynamic low slope effect. Based on these, the discriminant of different slopes taking on same response under the same dynamic input is put forward in this paper. At the same time, the paper studies distribution law of the rhythm extremum point of displacement, velocity and acceleration, and finds that there exists relationship of N = int among the slope height H, the number of the rhythm extremum VHlhro) point N and ffthre- Furthermore, the paper points out that if N^l, the response of the slope will be dynamic high slope effect; \fNnumerical simulation result. The Jinping hydropower station, located at the middle stream of Yalongjiang River and the west side of Jinpingdahewan, is the highest double-curvature arch dam planned in building in the world. However, deep fractures are well developed at the left abutment slope, so stability of the slope is a key engineering geological problem. The paper studies the dynamic stability of the slope. In the analysis, considering regional geology and regional seismology of the slope, combining with characteristics of the slope as well as the scale of the project, PGA(peak ground acceleration) of the site is decided (about 197.1cm/s2) by engineering seismology, and then the seismic input used in the slope dynamic analysis is determined. Comprehensive studies are carried out on the slope, especially to its deep fractures, and then the paper concludes that the deep fractures of slope are the result of the combination unloading effect of gravity and tectonic stress. At the same time, the failure model of the slope under dynamic input is attained. Based on these, the stability is comprehensively studied for section IV-IV with numerical simulation method as well as method putting forward in chapter 5. At last, the paper concludes that under dynamic input, the section will slide along fault f9, some deep fractures and fault/5 with certain permanent displacement, and this must be taken into consideration in the engineering.

高山峡谷地区大跨度桥梁岸坡稳定性分析、趋势预测及应用研究

In order to realize fast development of the national economy in a healthy way and coordinate progress with whole society, the country has implemented the strategy of development of the western region. An important action of finishing this strategic task is to accelerate the highway construction in the western region, join the western region and places along the coast, the river, the border with goods and materials, technology, and personnel interchanges, and then drive development of the local economy.The western region was influenced by the Himalaya Tectonization in Cenozoic, and the crust rose and became the plateau. In the course of rising, rivers cut down sharply to form a lot of high mountains and gorges.Because of topography and geomorphology, bridges in the traffic construction in the alpine gorge area are needed. Rivers have characteristics of large flow, fast velocity and high and steep river valley, so building a pier in the river is not only very difficult, but also making the cost increase. At the same time, the impact that the pier is corroded and the bridge base that is drawn to be empty by flow are apt to cause destruction of the pier. For those reasons, suspending bridge and cable-stay bridge are usually adopted with the single and large span. For the large span bridge, the pier foundation could receive ten thousand and more vertical strength, bending moment and near kiloton horizontal thrust.Because bank slope in the alpine gorge district is cut deeply and unsettled big, natural stability is worse under endogenic and exogenic force. When bank slope bears heavy vertical strength, bending moment and horizontal thrust facing the river, it will inevitably make the balance state of rock and soil mass change, bridge bank slope deform, and even destroyed. So the key problem at the time of the large span's bridge construction in the alpine gorge area is how to make it stable.So based on the spot investigation, the Engineering Geology Analysis Method is very important to grasp the bank slope stability. It can provide the bank slope stability macroscopic ally and qualitatively, and reference to the indoor calculation. The Engineering Geology Analysis Method is that by way of analyzing and investigating terms of bank slope instability, stability development trend, the ancient rock slide and devolution in the site, stability comprehensive evaluation primarily, current and future stability of bank slope is gotten, realizing the intention to serving the concrete engineering.After the Engineering Geology Analysis Method is applied to project instances of BeiPan River Bridge and BaLin River Bridge, results are accord with bank slope actual conditions, which proves sites are suited to building bridges from site stability.we often meet bank slope stability issues in the traffic construction in the alpine gorge areao Before the evaluation of the bank slope stability, the engineering geological condition is investigated first. After that, the next exploration target and geology measures are decided. So, the Engineering Geology Analysis Method that the investigation of the engineering geological condition is the main content is quite important in practice. The other evaluations of the bank slope stability are based on it. Because foundation receives very heavy load, for the big span's bridge in the alpine gorge area, a long pile of the large diameter (D^0.8m) is usually selected. In order to reflect rock mass's deformation properties under rock-socketed pile function, the author has used the FLAG30 software for rock and soil mass and done many numerical simulations. By them, the author launches the further investigation on deformation properties of bank slope under different slope angle, pile length, diameter, elastic modulus, load, bank slope's structure, etc. Some conclusion meaningful to the design and produce are obtained.

亚洲地形和海陆分布变化对东亚环境格局形成演化的作用：数值模拟与地质记录对比

Earlier studies on the distribution of geological environmental indicators in China revealed drastic changes from a zonal climate pattern (planetary-wave-dominant pattern) in the Paleogene to a monsoon-dominant one in the Neogene, which suggested an inception of the initial East-Asian summer monsoon. However, there are different views about the time and causes of the changes.Here, we attempt to compile a series of paleoenvironmental maps based on newly collected climate indicators from the literatures and chronologically constrained evidence of geological maps in order to re-examine the temporal and spatial evolution of climate belts in China during the Cenozoic with special emphasis on the changes of the arid belt. These indicators include mammalian fauna, coal, carbonate concretions, jarosite, salt, gypsum deposits and pollen assemblages etc, with chronological controls that we believe reliable. Pollen assemblages and mammalian fauna have been classified into three categories (arid, semi-arid/sub-humid, humid) to reflect the intensity of aridity/humidity. Salt, jarosite and gypsum deposits are classified as the arid indicators. Carbonate concretions and coal are classified into the semi-arid/sub-humid and humid one respectively. Paleoenvironmental maps at 8 time slices have been reconstructed. They are the Paleocene, Eocene, Oligocene, Miocene, Early Miocene, Middle Miocene, Late Miocene and Pliocene.And furthermore, we attempt to use IAP^AGCM to simulate the evolution of climate belts in emphasizing on the changes of the rain band, and compare the results with the paleoenvironmental maps in order to examine the causes of the drastic paleoenvironmental changes near the Oligocene/Miocene boundary. 36 sensitive numerical experiments are carried out using the IAP__AGCM to analyze the impacts of the uplift of the Himalayan-Tibetan complex, shrinkage of the Paratethys Sea, expansion of the South China Sea and the development of the polar ice sheets on rain band in China.The main conclusions are as follows:The obtained results essentially confirm the earlier conclusions about a zonal climate pattern in the Paleogene and a different pattern in the Neogene, and illustrate that a monsoon-dominant environmental pattern with inland aridity formed by the Early Miocene, which is temporally consist with the onset of eolian deposits in China.Cenozoic cooling and the formation of polar ice sheets are unlikely the main causes to the changes of environmental patterns mentioned above in China. But northern hemispheric cooling and the ice-sheets can intensify the Siberian High Pressure, and strengthen the winter monsoon circulations and enhance the aridity in the west part of China. These results support the earlier studies.Shrinkage of the Paratethys Sea and uplift of the Himalayan-Tibetan complex played important roles in strengthening the East Asian monsoon and induceing the above changes of environmental pattern, which is consistent with the earlier studies. Furthermore, "the monsoon-dominant pattern" appears when the Himalayan-Tibetan complex reaches to about 1000-2000 meters high and the Paratethys Sea retreats to the Turan Plate.4) Expansion of the South China Sea is another significant factor that drives the evolution of environmental patterns. We believe that the above three factors co-act and drive the change of the environmental patterns from a planetary-wave-dominant one to a monsoon-dominant one. However, the impacts of each factor vary by regions. The uplift mainly increases the humidity in Southwestern China and the aridity in northwestern country. The shrinkage mainly increases the humidity in Northern China and also enhances the aridity in the northwestern country. The expansion greatly increases the humidity in the south part of China.