The effects of random surface waves on the steady Ekman current solutions

The response of near-surface current profiles to wind and random surface waves are studied based on the approach of Jenkins [1989. The use of a wave prediction model for driving a near surface current model. Dtsch. Hydrogr. Z. 42,134-149] and Tang et al. [2007. Observation and modeling of surface currents on the Grand Banks: a study of the wave effects on surface currents. J. Geophys. Res. 112, C10025, doi:10.1029/2006JC004028]. Analytic steady solutions are presented for wave-modified Ekman equations resulting from Stokes drift, wind input and wave dissipation for a depth-independent constant eddy viscosity coefficient and one that varies linearly with depth. The parameters involved in the solutions can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water, and the solutions reduce to those of Lewis and Belcher [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans. 37, 313-351] when only the effects of Stokes drift are included. As illustrative examples, for a fully developed wind-generated sea with different wind speeds, wave-modified current profiles are calculated and compared with the classical Ekman theory and Lewis and Belcher's [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans 37, 313-351] modification by using the Donelan and Pierson [1987. Radar scattering and equilibrium ranges in wind-generated waves with application to scatterometry. J. Geophys. Res. 92, 4971-5029] wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. Illustrative examples for a fully developed sea and the comparisons between observations and the theoretical predictions demonstrate that the effects of the random surface waves on the classical Ekman current are important, as they change qualitatively the nature of the Ekman layer. But the effects of the wind input and wave dissipation on surface current are small, relative to the impact of the Stokes drift. (C) 2008 Elsevier Ltd. All rights reserved.

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.

Ship-board measurements and estimations of air-sea fluxes in the western tropical Pacific during TOGA COARE

Direct air-sea flux measurements were made on RN Kexue #1 at 40 degrees S, 156 degrees E during the Tropical Ocean Global Atmosphere (TOGA) Coupled Ocean-Atmospheric Response Experiment (COARE) Intensive Observation Period (IOP). An array of six accelerometers was used to measure the motion of the anchored ship, and a sonic anemometer and Lyman-alpha hygrometer were used to measure the turbulent wind vector and specific humidity. The contamination of the turbulent wind components by ship motion was largely removed by an improvement of a procedure due to Shao based on the acceleration signals. The scheme of the wind correction for ship motion is briefly outlined. Results are presented from data for the best wind direction relative to the ship to minimize flow distortion effects. Both the time series and the power spectra of the sonic-measured wind components show swell-induced ship motion contamination, which is largely removed by the accelerometer correction scheme, There was less contamination in the longitudinal wind component than in the vertical and transverse components. The spectral characteristics of the surface-layer turbulence properties are compared with those from previous land and ocean results, Momentum and latent heat fluxes were calculated by eddy correlation and compared to those estimated by the inertial dissipation method and the TOGA COARE bulk formula. The estimations of wind stress determined by eddy correlation are smaller than those from the TOGA COARE bulk formula, especially for higher wind speeds, while those from the bulk formula and inertial dissipation technique are generally in agreement. The estimations of latent heal flux from the three different methods are in reasonable agreement. The effect of the correction for ship motion on latent heat fluxes is not as large as on momentum fluxes.

Effect of Stokes drift on upper ocean mixing

Stokes drift is the main source of vertical vorticity in the ocean mixed layer. In the ways of Coriolis - Stokes forcing and Langmuir circulations, Stokes drift can substantially affect the whole mixed layer. A modified Mellor-Yamada 2.5 level turbulence closure model is used to parameterize its effect on upper ocean mixing conventionally. Results show that comparing surface heating with wave breaking, Stokes drift plays the most important role in the entire ocean mixed layer, especially in the subsurface layer. As expected, Stokes drift elevates both the dissipation rate and the turbulence energy in the upper ocean mixing. Also, influence of the surface heating, wave breaking and wind speed on Stokes drift is investigated respectively. Research shows that it is significant and important to assessing the Stokes drift into ocean mixed layer studying. The laboratory observations are supporting numerical experiments quantitatively.

Wave breaking on turbulent energy budget in the ocean surface mixed layer

As an important physical process at the air-sea interface, wave movement and breaking have a significant effect on the ocean surface mixed layer (OSML). When breaking waves occur at the ocean surface, turbulent kinetic energy (TKE) is input downwards, and a sublayer is formed near the surface and turbulence vertical mixing is intensively enhanced. A one-dimensional ocean model including the Mellor-Yamada level 2.5 turbulence closure equations was employed in our research on variations in turbulent energy budget within OSML. The influence of wave breaking could be introduced into the model by modifying an existing surface boundary condition of the TKE equation and specifying its input. The vertical diffusion and dissipation of TKE were effectively enhanced in the sublayer when wave breaking was considered. Turbulent energy dissipated in the sublayer was about 92.0% of the total depth-integrated dissipated TKE, which is twice higher than that of non-wave breaking. The shear production of TKE decreased by 3.5% because the mean flow fields tended to be uniform due to wave-enhanced turbulent mixing. As a result, a new local equilibrium between diffusion and dissipation of TKE was reached in the wave-enhanced layer. Below the sublayer, the local equilibrium between shear production and dissipation of TKE agreed with the conclusion drawn from the classical law-of-the-wall (Craig and Banner, 1994).

Submesoscale activity over the shelf of the northern South China Sea in summer: simulation with an embedded model

We applied a primitive equation ocean model to simulate submesoscale activities and processes over the shelf of the northern South China Sea (NSCS) with a one-way nesting technology for downscaling. The temperature and density fields showed that submesoscale activities were ubiquitous in the NSCS shelf. The vertical velocity was considerably enhanced in submesoscale processes and could reach an average of 58 m per day in the subsurface. At this point, the mixed layer depth also was deepened along the front, and the surface kinetic energy also increased with the intense vertical movement induced by submesoscale activity. Thus, submesoscale stirring/mixing is important for tracers, such as temperature, salinity, nutrients, dissolved organic, and inorganic carbon. This result may have implication for climate and biogeochemical investigations.

潮波模型的优化开边界方法研究

数值模式是潮波研究的一种有利手段，但在研究中会面临各种具体问题，包括开边界条件的确定、底摩擦系数和耗散系数的选取等。数据同化是解决这些问题的一种途径，即利用有限数量的潮汐观测资料对潮波进行最优估计，其根本目的是迫使模型预报值逼近观测值，使模式不要偏离实际情况太远。本文采用了一种优化开边界方法，沿着数值模型的开边界优化潮汐水位信息，目的是设法使数值解在动力约束的意义下接近观测值，获得研究区域的潮汐结果。边界值由指定优化问题的解来定，以提高模拟区域的潮汐精度，最优问题的解是基于通过开边界的能量通量的变化，处理开边界处的观测值与计算值之差的最小化。这里提供了辐射型边界条件，由Reid 和Bodine（本文简称为RB）推导，我们将采用的优化后的RB方法（称为ORB）是优化开边界的特殊情况。 本文对理想矩形海域（ E- E， N- N， 分辨率 ）进行了潮波模拟，有东部开边界，模式采用ECOM3D模式。对数据结果的误差分析采用，振幅平均偏差，平均绝对偏差，平均相对误差和均方根偏差四个值来衡量模拟结果的好坏程度。 需要优化入开边界的解析潮汐值本文采用的解析解由方国洪《海湾的潮汐与潮流》（1966年）方法提供，为验证本文所做的解析解和方文的一致，本文做了其第一个例子的关键值a,b,z,结果与其结果吻合的相当好。但略有差别，分析的可能原因是两法在具体迭代方案和计算机保留小数上有区别造成微小误差。另外，我们取m=20，得到更精确的数值，我们发现对前十项的各项参数值，取m=10,m=20各项参数略有改进。当然我们可以获得m更大的各项参数值。 同时为了检验解析解的正确性讨论m和l变化对边界值的影响，结果指出，增大m，m=20时，u的模最大在本身u1或u2的模的6%；m=100时，u的模最大在本身u1或u2的模的4%；m再增大，m=1000时，u的模最大在本身u1或u2的模的4%，改变不大。当l<1时， =0处u的模最大为2。当l=1时， =0处u的模最大为0.1，当l>1时，l越大，u的模越小，当l=10时，u的模最大为0.001，可以认为为0。 为检验该优化方法的应用情况，我们对理想矩形区域进行模拟，首先将本文所采用的优化开边界方法应用于30m的情况，在开边界优化入开边界得出模式解，所得模拟结果与解析解吻合得相当好，该模式解和解析解在整个区域上，振幅平均绝对偏差为9.9cm,相位平均绝对偏差只有4.0 ，均方根偏差只有13.3cm，说明该优化方法在潮波模型中有效。 为验证该优化方法在各种条件下的模拟结果情况，在下面我们做了三类敏感性试验： 第一类试验：为证明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，我们来比较ORB条件与RB条件的优劣，我们模拟用了两个不同的摩擦系数，k分别为：0,0.00006。 结果显示，针对不同摩擦系数，显示在开边界上使用ORB条件的解比使用RB条件的解无论是振幅还是相位都有显著改善，两个试验均方根偏差优化程度分别为84.3%，83.7%。说明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，大大提高了模拟水平。上述的两个试验得出， k=0.00006优化结果比k=0的好。 第二类试验，使用ORB条件确定优化开边界情况下，在东西边界加入出入流的情况，流考虑线性和非线性情况，结果显示，加入流的情况，潮汐模拟的效果降低不少，流为1Sv的情况要比5Sv的情况均方根偏差相差20cm,而不加流的情况只有0.2cm。线性流和非线性流情况两者模式解相差不大，振幅，相位各项指数都相近, 说明流的线性与否对结果影响不大。 第三类试验，不仅在开边界使用ORB条件，在模式内部也使用ORB条件，比较了内部优化和不优化情况与解析解的偏差。结果显示，选用不同的k,振幅都能得到很好的模拟，而相位相对较差。另外，在内部优化的情况下，考虑不同的k的模式解， 我们选用了与解析解相近的6个模式解的k,结果显示，不同的k,振幅都能得到很好的模拟，而相位较差。 总之，在开边界使用ORB条件比使用RB条件好，振幅相位都有大幅度改进，在加入出入流情况下，流的大小对模拟结果有影响，但线形流和非线性流差别不大。内部优化的结果显示，模式采用不同的k都能很好模拟解析解的振幅。

Response of chlorophyll fluorescence to dynamic light in three alpine species differing in plant architecture

A study was carried out to examine the effect of dynamic photosynthetically active photon flux density (PPFD) on photoinhibition and energy use in three herbaceous species, prostrate Saussurea superba, erect-leaved S. katochaete, and half-erect-leaved Gentiana straminea, from the Qinghai-Tibet Plateau. Chlorophyll fluorescence response was measured under each of three sets of high-low PPFD combinations: 1700-0, 1400-300, and 1200-500 mu mol m(-2) s(-1), illuminating in four dynamic frequencies: 1, 5, 15, and 60 cycles per 2 h. The total light exposure time was 2h and the integrated PPFD was the same in all treatments. The highest frequency of PPFD fluctuation resulted in the lowest photochemical activity, the highest level of non-photochemical quenching, and the greatest decrease of F-v/F-m (maximal photochemical efficiency of PSII). The 5 and 15 cycles per 2h treatments resulted in higher photochemical activity than the 1 cycle per 2h treatment. The 1700-0 PPFD combination led to the lowest photochemical activity and more serious photoinhibition in all species. S. superba usually exhibited the highest photochemical activity and CO2 uptake rate, the lowest reduction of F-v/F-m,F- and the smallest fraction of energy in thermal dissipation. With similar fractions of thermal dissipation, S. katochaete had relatively less photoinhibition than G. straminea owing to effective F-o quenching. The results suggest that high frequency of fluctuating PPFD generally results in photoinhibition, which is more serious under periods of irradiation with high light intensity. (c) 2005 Elsevier B.V. All rights reserved.

基于SIM1-2的现场总线仪表电路设计

媒体结合单元(MAU)是现场总线仪表中的重要部件。介绍一种新型的媒体结合单元电路器件SIM1-2,该器件符合IEC61158-2数据链路层协议规范。在详细分析该器件的系统结构和主要功能的基础上,介绍其在总线仪表系统中的应用。并从资源需要,处理速度和功耗等方面分析单片机、通讯控制器、MAU电路的选择等关键问题。

边坡工程的极限分析上限法

Theory of limit analysis include upper bound theorem and lower bound theorem. To deal with slope stability analysis by limit analysis is to approximate the real solution from upper limit and lower limit. The most used method of limit analysis is upper bound theorem, therefore it is often applied to slope engineering in many cases. Although upper bound approach of limit analysis can keep away from vague constitutive relation and complex stress analyses, it also can obtain rigorous result. Assuming the critical surface is circular slip surface, two kinematically admissible velocity fields for perpendicular slice method and radial slice method can be established according to the limit analysis of upper bound theorem. By means of virtual work rate equation and strength reduction method, the upper-bound solution of limit analysis for homogeneous soil slope can be obtained. A log-spiral rotational failure mechanism for homogeneous slope is discussed from two different conditions which represent the position of shear crack passing the toe and below the toe. In the dissertition, the author also establishes a rotational failure mechanics with combination of different logarithmic spiral arcs. Furthermore, the calculation formula of upper bound solution for inhomogeneous soil slope stability problem can be deduced based on the upper bound approach of rigid elements. Through calculating the external work rate caused by soil nail, anti-slide pile, geotechnological grid and retaining wall, the upper bound solution of safety factor of soil nail structure slope, slip resistance of anti-slide pile, critical height of reinforced soil slope and active earth pressure of retaining wall can be obtained by upper bound limit analysis method. Taking accumulated body slope as subject investigated, with study on the limit analysis method to calculate slope safety factor, the kinematically admissible velocity fields of perpendicular slice method for slope with broken slip surface is proposed. Through calculating not only the energy dissipation rate produced in the broken slip surfaces and the vertical velocity discontinuity, but also the work rate produced by self-weight and external load, the upper bound solution of slope with broken slip surface is deduced. As a case study, the slope stability of the Sanmashan landslide in the area of the Three Gorges reservoir is analyzed. Based on the theory of limit analysis, the upper bound solution for rock slope with planar failure surface is obtained. By means of virtual work-rate equation, energy dissipation caused by dislocation of thin-layer and terrane can be calculated; furthermore, the formulas of safety factor for upper bound approach of limit analysis can be deduced. In the end, a new computational model of stability analysis for anchored rock slope is presented after taking into consideration the supporting effect of rock-bolts, the action of seismic force and fissure water pressure. By using the model, not only the external woke-rate done by self-weight, seismic force, fissure water pressure and anchorage force but also the internal energy dissipation produced in the slip surface and structural planes can be totally calculated. According to the condition of virtual work rate equation in limit state, the formula of safety factor for upper bound limit analysis can be deduced.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

基于工程地质力学的岩土工程加固理论、技术与实践

In the engineering reinforcement of-rock and soil mass, engineers must consider how to obtain better reinforcing effect at the cost of less reinforcing expense, which, in fact, is the aim of reinforcement design. In order to accomplish the purpose, they require not only researching the material used to reinforce and its structure, but also taking into account of several important geological factors, such as the structure and property of rock and soil mass. How to improve the reinforcing effect according to engineering geomechanical principle at the respect of the reinforcement of engineering soil and rock mass is studied and discussed in this paper. The author studies the theory, technology and practice of geotechnical reinforcement based on engineering geomechanics, taking example for the soil treatment of Zhengzhou Airport, the effect analysis of reinforcement to the slope on the left bank of Wuqiangxi Hydropower Station and the reinforcing design of the No. 102 Landslide and unique sand-slide slope on the Sichuan-Tibet Highway. The paper is comprised of two parts for the convenience of discussion. In the first part, from the first chapter to the fifth chapter, trying to perform the relevant research and application at the viewpoint of soil mass engineering geomechanics, the author mainly discusses the study of reinforcing soft ground soil through dynamical consolidation and its application. Then, in the second part, from the sixth chapter to the eleventh chapter, the study of new technologies in the rock slope reinforcement and their application are discussed. The author finds that not only better reinforcing effect can be gained in the research where the principle and method of rock mass engineering geomechanics is adopted, but also new reinforcing technologies can be put forward. Zhengzhou Airport is an important one in central plains. It lies on Yellow River alluvial deposit and the structure of stratum is complex and heterogeneous. The area of airport is very large, which can result in differential settlement easily, damage of airport and aircraft accident, whereas, there are no similar experiences to dispose the foundation, so the foundation treatment become a principal problem. During the process of treatment, the method of dynamic compaction was adopted after compared with other methods using the theory of synthetic integration. Dynamic compaction is an important method to consolidate foundation, which was successfully used in the foundation of Zhengzhou Airport. For fill foundation, controlling the thickness of fill so as to make the foundation treatment can reach the design demand and optimum thickness of the fill is a difficult problem. Considering this problem, the author proposed a calculation method to evaluate the thickness of fill. The method can consider not only the self-settlement of fill but also the settlement of the ground surface under applied load so as to ensure the settlement occurred during the using period can satisfy the design demand. It is proved that the method is correct after using it to choose reasonable energy of dynamic compaction to treat foundation. At the same time, in order to examine the effect of dynamic compaction, many monitor methods were adopted in the test such as static loading test, modulus of resilience test, deep pore pressure -test, static cone penetration test and the variation of the pore volume measurement. Through the tests, the author summarized the discipline of the accumulation and dissipation of pore pressure in Yellow River alluvial deposit under the action of dynamic compaction, gave a correct division of the property change of silt and clay under dynamic compaction, determined the bearing capacity of foundation after treatment and weighted the reinforcing effect of dynamic consolidation from the variation of the soil particle in microcosmic and the parameter of soil mass' density. It can be considered that the compactness of soil is in proportion to the energy of dynamic compaction. This conclusion provided a reference to the research of the "Problem of Soil Structure-the Central Problem of Soil Mechanics in 21 Century ". It is also important to strengthen rock mass for water conservancy and electric power engineering. Slip-resistance pile and anchoring adit full of reinforced concrete are usually adopted in engineering experience to strengthen rock mass and very important for engineering. But there also some deficiency such as the weakest section can't be highlighted, the monitor is inconvenient and the diameter of pile and adit is very large etc. The author and his supervisor professor Yangzhifa invented prestressed slip-resistance pile and prestressed anchoring adit full of reinforced concrete, utilizing the advantage that the prestressed structure has better anti-tensile characteristic (this invention is to be published). These inventions overcome the disadvantages of general slip-resistance pile and anchoring adit full of reinforced concrete and have the functions of engineering prospecting, strengthening, drainage and monitor simultaneous, so they have better strengthened effect and be more convenient for monitor and more economical than traditional methods. Drainage is an important factor in treatments of rock mass and slop. In view of the traditional drainage method that drainage pore often be clogged so as to resulted in incident, professor Yangzhifa invented the method and setting of guide penetration by fiber bundle. It would take good effect to use it in prestressed slip-resistance pile and anchoring adit full of reinforced concrete. In this paper, the author took example for anchoring adit full of reinforced concrete used to strengthen Wuqiangxi left bank to simulate the strengthened effect after consolidated by prestressed slip-resistance pile, took example for 102 landslide occurred along Sichuan-Tibet highway to simulate the application of slip-resistance pile and the new technology of drainage. At the same time the author proposed the treatment method of flowing sand in Sichuan-Tibet highway, which will benefit the study on strengthening similar engineering. There are five novelties in the paper with the author's theoretical study and engineering practice: 1. Summarizing the role of pore water pressure accumulation and dissipation of the Yellow River alluvial and diluvial soil under the action of dynamical consolidation, which has instructive significance in the engineering construction under the analogical engineering geological conditions in the future. It has not been researched by the predecessors. 2. Putting forward the concept of density D in microcosmic based on the microcosmical structure study of the soil sample. Adopting D to weight the reinforcing effect of dynamic consolidation is considered to be appropriate by the means of comparing the D values of Zhengzhou Airport's ground soil before with after dynamically consolidating reinforcement, so a more convenient balancing method can be provided for engineering practice. 3. According to the deep research into the soil mass engineering geology, engineering rock and soil science, soil mechanics, as well as considerable field experiments, improving the consolidating method in airport construction, from the conventional method, which is dynamically compactmg original ground surface firstly, then filling soil and dynamically layer-consolidating or layer-compacting at last to the upgraded method, which is performing dynamical consolidation after filling soil to place totally at the extent of the certain earth-filling depth. The result of the dynamical consolidation not only complies with the specifications, but also reduces the soil treatment investment by 10 million RMB. 4. Proposing the method for calculating the height of the filled soil by the means of estimating the potential displacement produced in the original ground surface and the filled earth soil under the possible load, selecting the appropriate dynamically-compacting power and determining the virtual height of the filled earth soil. The method is proved to be effective and scientific. 5. According to the thought of Engineering Geomechanics Metal-Synthetic Methodology (EGMS), patenting two inventions (to the stage of roclamation, with Professor Yang Zhi-fa, the cooperative tutor, and etc.) in which multi-functions, engineering geological investigation, reinforcement, drainage and strength remedy, are integrated all over in one body at the viewpoint of the breakage mechanism of the rock slope.

陡坡带砂砾岩体油气成藏模式研究——以东营凹陷西北部陡坡带为例

Based on the study of sequence stratigraphy, modern sedimentary, basin analysis, and petroleum system in abrupt slop of depression, this paper builds sedimentary system and model, sandy bodies distribution, and pool-forming mechanism of subtle trap. There are some conclusions and views as follows. By a lot of well logging and seismic analysis, the author founded up the sequence stratigraphic of the abrupt slope, systematically illustrated the abrupt slope constructive framework, and pointed out that there was a special characteristics which was that south-north could be divided to several fault block and east-west could be carved up groove and the bridge in studying area. Based all these, the author divided the studying area to 3 fault block zone in which because of the groove became the basement rock channel down which ancient rivers breathed into the lake, the alluvial fan or fan delta were formed. In the paper, the author illustrated the depositional system and depositional model of abrupt slope zone, and distinguished 16 kinds of lithofacies and 3 kinds of depositional systems which were the alluvial fan and fan-delta system, lake system and the turbidite fan or turbidity current deposition. It is first time to expound completely the genetic pattern and distributing rule of the abrupt slope sandy-conglomeratic fan bodies. The abrupt slope sandy-conglomeratic fan bodies distribute around the heaves showing itself circularity shape. In studying area, the sandy-conglomeratic fan bodies mainly distribute up the southern slope of Binxian heave and Chenjiazhuang heave. There mainly are these sandy-conglomeratic fan body colony which distributes at a wide rage including the alluvial fan, sub-water fluvial and the turbidite fan or the other turbidity current deposition in the I fault block of the Wangzhuang area. In the II fault block there are fan-delta front and sub-water fluvial. And in the Binnan area, there mainly are those the alluvial fan (down the basement rock channel) and the sandy-conglomeratic fan body which formed as narrowband sub-water fluvial (the position of bridge of a nose) in the I fault block, the fan-delta front sandy-conglomeratic fan body in the H fault block and the fan-delta front and the turbidity current deposition sandy-conglomeratic fan body in the m fault block. Based on the reservoir outstanding characteristics of complex classic composition and the low texture maturity, the author comparted the reservoir micro-structure of the Sha-III and Sha-IV member to 4 types including the viscous crude cementation type, the pad cementation type, the calcite pore-funds type and the complex filling type, and hereby synthetically evaluated 4 types sandy- conglomeratic fan body reservoir. In the west-north abrupt slope zone of Dongying Depression, the crude oil source is belonging to the Sha-III and Sha-IV member, the deep oil of Lijin oilfield respectively come from the Sha-III and Sha-IV member, which belongs to the autogeny and original deposition type; and the more crude oil producing by Sha-IV member was migrated to the Wangzhuan area and Zhengjia area. The crude oil of Binnan oil-field and Shanjiasi oil-field belongs to mixed genetic. It is the first time to illustrate systematically the genetic of the viscous crude that largely being in the studying area, which are that the dissipation of the light component after pool-forming, the biological gradation action and the bath-oxidation action, these oil accumulation belonging to the secondary viscous crude accumulation. It is also the first time to compart the studying area to 5 pool-forming dynamical system that have the characteristic including the common pressure and abnormal pressure system, the self-fountain and other-fountain system and the closing and half-closing system etc. The 5 dynamical systems reciprocally interconnected via the disappearance or merger of the Ethology and the fluid pressure compartment zone, the fault and the unconformity surface, hereby formed duplicated pattern oil-gas collecting zone. Three oil-gas pool-forming pattern were founded, which included the self-fountain side-direction migrated collecting pattern, the self-fountain side-direction ladder-shape pool-forming pattern and the other-fountain pressure releasing zone migrated collecting pattern. A series of systemic sandy-conglomeratic fan bodies oil-gas predicting theory and method was founded, based on the groove-fan corresponding relation to confirm the favorable aim area, according as the characteristic of seismic-facies to identify qualitatively the sandy-conglomeratic fan bodies or its scale, used the temporal and frequency analysis technique to score the interior structure of the sandy- conglomeratic fan bodies, applied for coherent-data system analysis technology to describe the boundary of the sandy-conglomeratic fan bodies, and utilized the well logging restriction inversion technique to trace quantificational and forecast the sandy-conglomeratic fan bodies. Applied this technique, totally 15 beneficial sandy-conglomeratic fan bodies were predicted, in studying area the exploration was preferably guided, and the larger economic benefit and social benefit was acquired.

微功耗海底地震仪的研制

Our motherland has large area of maritime space. Searching and developing ocean becomes more and more important. So Ocean Bottom Seismometer (OBS) as an absolutely necessary equipment can be used in many oceanic fields. OBS not only is an important instrument for discovering structure of lithosphere of ocean bottom, but also plays a main role of oceanic geophysical exploration. The paper introduces my relational work. The MCI micro-power broad frequency seismometer was developed independently. Its power dissipation is less than 300mW. It has some merits including miniature volumeN light mass and cheap price. It is an ideal device not only for the collection high-resolution natural seismic data, but also for the fields of seismic sounding and engineering seismology. Many new high technique were applied to develop this instrument including over-sampling A/D converter, high performance 32bit Micro Process Unit and Flash memory with smart-media interface. Base on the achievement, I have accomplished the showpiece of OBS, which is applied to the deepwater oil and gas geophysical exploration. Because of micro-power dissipation, the seismograph and the sonar releaser can be integrated into a sphere cabin. By this means, the instrument's frequency of resonance and frequency of couple are improved obviously. The data acquisition system of OBS is improved from MCI seismometer. The capacity of flash memory is enlarge from 1G bytes to 8G bytes. The advance MPU in data acquisition system is used to integrate other function modules such as sonar, GPS, compass and digital transmitter.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.