Statistical distribution of depth-integrated local horizontal momentum for second-order random ocean waves in finite water depth

Based on the second-order random wave solutions of water wave equations in finite water depth, statistical distributions of the depth- integrated local horizontal momentum components are derived by use of the characteristic function expansion method. The parameters involved in the distributions can be all determined by the water depth and the wave-number spectrum of ocean waves. As an illustrative example, a fully developed wind-generated sea is considered and the parameters are calculated for typical wind speeds and water depths by means of the Donelan and Pierson spectrum. The effects of nonlinearity and water depth on the distributions are also investigated.

Statistical distribution of wave-surface elevation for second-order random directional ocean waves in finite water depth

Based on the second-order random wave solutions of water wave equations in finite water depth, a statistical distribution of the wave-surface elevation is derived by using the characteristic function expansion method. It is found that the distribution, after normalization of the wave-surface elevation, depends only on two parameters. One parameter describes the small mean bias of the surface produced by the second-order wave-wave interactions. Another one is approximately proportional to the skewness of the distribution. Both of these two parameters can be determined by the water depth and the wave-number spectrum of ocean waves. As an illustrative example, we consider a fully developed wind-generated sea and the parameters are calculated for various wind speeds and water depths by using Donelan and Pierson spectrum. It is also found that, for deep water, the dimensionless distribution reduces to the third-order Gram-Charlier series obtained by Longuet-Higgins [J. Fluid Mech. 17 (1963) 459]. The newly proposed distribution is compared with the data of Bitner [Appl. Ocean Res. 2 (1980) 63], Gaussian distribution and the fourth-order Gram-Charlier series, and found our distribution gives a more reasonable fit to the data. (C) 2002 Elsevier Science B.V. All rights reserved.

A new method for the estimation of oceanic mixed-layer depth using shipboard X-band radar images

Shipboard X-band radar images acquired on 24 June 2009 are used to study nonlinear internal wave characteristics in the northeastern South China Sea. The studied images show three nonlinear internal waves in a packet. A method based on the Radon Transform technique is introduced to calculate internal wave parameters such as the direction of propagation and internal wave velocity from backscatter images. Assuming that the ocean is a two-layer finite depth system, we can derive the mixed-layer depth by applying the internal wave velocity to the mixed-layer depth formula. Results show reasonably good agreement with in-situ thermistor chain and conductivity-temperature-depth data sets.

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.

Current status of bacteriological parameters and DOC/POC in Xiamen coastal waters

The surface and bottom waters samples were collected from six locations in Xiamen western sea. The quantified estimation of bacterial production (H-3-thymidine method) and observation of bacterial heterotrophic activity (C-14-glucose method) have been made in order to have a better understanding of the role of marine bacteria and their activities. The results showed that the mean value of bacterial heterotrophic activity was 9 X 10(8) cells/(L. h) in the surface waters and 2.6 X 10(8) cells/(L. h) in the bottom waters. The mean value of bacterial production was 38 X 108 cells/( L. h) in the surface waters and 7.1 X 10(8) cells/(L. h) in the bottom waters. The relationship between bacterial production, heterotrophic activity, PCC and DOC measured during this survey were discussed. The good understanding of the relationship between bacteria activity and total coliform was addressed.

Studies on two coordination polymers [M(mu(4)-pz25dc)](n) (M = Cd or Zn, pz25dc=pyrazine-2,5-dicarboxylato) with three-dimensional pillared-layer three-nodal framework: Synthesis, structural characterization, strong optical non-linearities and optical limiting properties

Two isomorphous new candidates [M(mu(4)-pz25dc)](n) (M = Cd, 1; Zn, 2; pz25dc = pyrazine-2,5-dicarboxylato)for nonlinear optical (NLO) materials have been synthesized hydrothermally and characterized crystallographically as pillared-layer three-nodal frameworks with one four-connected metal nodes and two crystallographically different four-connected ligand nodes. Their optical non-linearities are measured by the Z-Scan technique with an 8 ns pulsed laser at 532 nm. These two coordination polymers both exhibit strong NLO absorptive abilities [alpha(2) = (63 +/- 6) x 10 (12) mW (1) 1, ( 46 +/- 6) x 10 (11) mW (1) 2] and effective self-focusing performance [n(2) = (67 +/- 5) x 10 (18) 1, (13 +/- 3) x 10 (18) m(2) W (1) 2] in 1.02 x 10 (4) 1 and 1.05 x 10 (4) mol dm (3) 2 DMF solution separately. The values of the limiting threshold are also measured from the optical limiting experimental data. The heavy atom effect plays important role in the enhancement of optical non-linearities and optical limiting properties. (C) 2009 Elsevier B. V. All rights reserved.

Spatial variations of size-fractionated chlorophyll, cyanobacteria and heterotrophic bacteria in the Central and Western Pacific

Geographic and vertical variations of size-fractionated (0.2-1 mu m, 1-10 mu m, and >10 mu m) Chlorophyll a (Chl.a) concentration, cyanobacteria abundance and heterotrophic bacteria abundance were investigated at 13 stations from 4 degrees S, 160 degrees W to 30 degrees N, 140 degrees E in November 1993. The results indicated a geographic distribution pattern of these parameters with instances of high values occurring in the equatorial region and offshore areas, and with instance of low values occurring in the oligotrophic regions where nutrients were almost undetectable. Cyanobacteria showed the highest geographic variation (ranging from 27x10(3) to 16,582x10(3) cell l(-1)), followed by Chl.a (ranging from 0.048 to 0.178 mu g l(-1)), and heterotrophic bacteria (ranging from 2.84x10(3) to 6.50 x 10(5) cell l(-1)). Positive correlations were observed between nutrients and Chl.a abundance. Correspondences of cyanobacteria and heterotrophic bacteria abundances to nutrients were less significant than that of Chl.a. The total Chl.a was accounted for 1.0-30.9%, 35.9-53.7%, and 28.1-57.3% by the >10 mu m, 1-10 mu m and 0.2-1 mu m fractions respectively. Correlation between size-fractionated Chl.a and nutrients suggest that the larger the cell size, the more nutrient-dependent growth and production of the organism. The ratio of pheophytin to chlorophyll implys that more than half of the > 10 mu m and about one third of the 1-10 mu m pigment-containing particles in the oligotrophic region were non-living fragments, while most of the 1-10 mu m fraction was living cells. In the depth profiles, cyanobacteria were distributed mainly in the surface layer, whereas heterotrophic bacteria were abundant from surface to below the euphotic zone. Chl.a peaked at the surface layer (0-20 m) in the equatorial area and at the nitracline (75-100 m) in the oligotrophic regions. Cyanobacteria were not the principle component of the picoplankton. The carbon biomass ratio of heterotroph to phytoplankton was greater than 1 in the eutrophic area and lower than 1 in oligotrophic waters.

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

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

地质导向钻井关键技术及方法

On the basis of analyzing the principle and realization of geo-steering drilling system, the key technologies and methods in it are systematically studied in this paper. In order to recognize lithology, distinguish stratum and track reservoirs, the techniques of MWD and data process about natural gamma, resistivity, inductive density and porosity are researched. The methods for pre-processing and standardizing MWD data and for converting geological data in directional and horizontal drilling are discussed, consequently the methods of data conversion between MD and TVD and those of formation description and adjacent well contrast are proposed. Researching the method of identifying sub-layer yields the techniques of single well explanation, multi-well evaluation and oil reservoir description. Using the extremum and variance clustering analysis realizes logging phase analysis and stratum subdivision and explanation, which provides a theoretical method and lays a technical basis for tracing oil reservoirs and achieving geo-steering drilling. Researching the technique for exploring the reservoir top with a holdup section provides a planning method of wellpath control scheme to trace oil and gas reservoir dynamically, which solves the problem of how to control well trajectory on condition that the layer’s TVD is uncertain. The control scheme and planning method of well path for meeting the demands of target hitting, soft landing and continuous steering respectively provide the technological guarantee to land safely and drill successfully for horizontal, extended-reach and multi-target wells. The integrative design and control technologies are researched based on geology, reservoir and drilling considering reservoir disclosing ratio as a primary index, and the methods for planning and control optimum wellpath under multi-target restriction, thus which lets the target wellpath lie the favorite position in oil reservoir during the process of geo-steering drilling. The BHA (bottomhole assembly) mechanical model is discussed using the finite element method, and the BHA design methods are given on the basis of mechanical analyses according to the shape of well trajectory and the characteristics of BHA’s structure and deformation. The methods for predicting the deflection rate of bent housing motors and designing their assemblies are proposed based on the principle of minimum potential energy, which can clearly show the relation between the BHA’s structure parameters and deflection rate, especially the key factors’ effect to the deflection rate. Moreover, the interaction model between bit and formation is discussed through the process of equivalent formation and equivalent bit considering the formation anisotropy and bit anisotropy on the basis of analyzing the influence factors of well trajectory. Accordingly, the inherence relationship among well trajectory, formation, bit and drilling direction is revealed, which lays the theory basis and technique for predicting and controlling well trajectory.

悬索桥隧道式锚碇与围岩系统的破坏模式研究

The anchorages are unparalleled structures only in a suspension bridge, and as main bearing facilities, play an important role in connecting the superstructures and the ground. The tunnel anchorage, as one alternative type of the anchorages, has more advantages over its counterpart, the gravity anchorage. With the tunnel anchorages adopted, not only can surface excavation be reduced to protect the environment, and natural condition of the rock be utilized and potential bearing capacity of surrounding rock be mobilized to save engineering cost, but also the technological predominance of auxiliary engineering measures, such as prestressed concrete, anchoring piles, rock anchors and collar beam between the two separated anchorages, can be easily cooperated to work together harmoniously under the circumstances of poor rock quality. There are plentiful high mountains and deep canyons in west part of China, and long-span bridge construction is inevitably encountered in order to realize leapfrogging development of the transportation infrastructure. Western mountainous areas usually possess the conditions for constructing tunnel anchorages, and therefore, the tunnel anchorages, which are conformed to the conception of resource conservative and sustainable society, extremely have application and popularization value in western underdeveloped region. The scientific and technological problem about the design, construction and operation of tunnel anchorages should be further investigated. Combining the engineering of western tunnel anchorages for the Balinghe Suspension Bridge, this paper probed into the survey method and in-situ test method for tunnel anchorages, scientific rock quality evaluation of surrounding rock to provide reasonable physical and mechanical parameters for design, construction and operation of tunnel anchorages, bearing capacity estimation for tunnel anchorage, deformation prediction of the anchorage-rockmass system, tunnel-anchorage slope stability analysis and the evaluation of excavation stability and degree of safety of the anchorage tunnel. The following outcomes were obtained: 1. Materials of tunnel anchorages of suspension bridge built (and in progress) at home and abroad were systematically sorted out, with the engineering geological condition and geomechanical property of surrounding rock around the anchorage tunnel, the design size of anchorages and the construction method of anchorage tunnel paid more emphasis on, to unveil the internal relationship between the engineering geological conditions of surrounding rock and the design size and axis angle of anchorages and provide references for future design, construction and study of tunnel anchorages. 2. Physical and mechanical parameters were recommended based on three domestic and foreign methods of rock quality evaluation. 3. In-situ tests, adopting the back-thrust method, of two kinds of reduced scale model, 1/30 and 1/20, for the tunnel anchorages were conducted in the declining exploration drift with rock mass at the test depth being the same as surrounding rock around real anchorages, and reliable field rockmass displacement data were acquired. Attenuation relation between the increment of distance from the anchorage and the decrement of rockmass displacement under maximum test load, and influential scope suffered by anchorage load were obtained. 4. Using similarity theory, the magnitude of real anchorage and rockmass displacement under design load and degree of safety of the anchorage system were deduced. Furthermore, inversion analysis to deformation modulus of slightly weathered dolomite rock, the surrounding rock of anchorage tunnel, was performed by the means of numerical simulation. 5. The influential law of the geometrical size to the limit bearing capacity of tunnel anchorage was studied. 6. Based on engineering geological survey data, accounting for the combination of strata layer and adverse discontinuities, the failure patterns of tunnel anchorage slope were divided into three modes: sliding of splay saddle pier slope, superficial-layer slippage, and deep-layer slippage. Using virtual work principle and taking anchorage load in account, the stability of the three kinds of failure patterns were analyzed in detail. 7. The step-by-step excavation of anchorage tunnel, the numerical overload and the staged decrement of rock strength parameters were numerically simulated to evaluate the excavation stability of surrounding rock around anchorage tunnel, the overload performance of tunnel anchorage, and the safety margin of strength parameters of the surrounding rock.

基于随机介质理论的多尺度储层非均质性定量描述

Application of long-term exploration for oil and gas shows that the reservoir technology of prediction is one of the most valuable methods. Quantitative analysis of reservoir complexity is also a key technology of reservoir prediction. The current reservoir technologies of prediction are based on the linear assumption of various physical relationships. Therefore, these technologies cannot handle complex reservoirs with thin sands, high heterogeneities in lithological composition and strong varieties in petrophysical properties. Based on the above-mentioned complex reservoir, this paper conducts a series of researches. Both the comprehending and the quantitative analysis of reservoir heterogeneities have been implemented using statistical and non-linear theories of geophysics. At the beginning, the research of random media theories about reservoir heterogeneities was researched in this thesis. One-dimensional (1-D) and two-dimensional (2-D) random medium models were constructed. The autocorrelation lengths of random medium described the mean scale of heterogeneous anomaly in horizontal and deep directions, respectively. The characteristic of random medium models were analyzed. We also studied the corresponding relationship between the reservoir heterogeneities and autocorrelation lengths. Because heterogeneity of reservoir has fractal nature, we described heterogeneity of reservoir by fractal theory based on analyzing of the one-dimensional (1-D) and two-dimensional (2-D) random medium models. We simulated two-dimensional (2-D) random fluctuation medium in different parameters. From the simulated results, we can know that the main features of the two-dimensional (2-D) random medium mode. With autocorrelation lengths becoming larger, scales of heterogeneous geologic bodies in models became bigger. In addition, with the autocorrelation lengths becoming very larger, the layer characteristic of the models is very obvious. It would be difficult to identify sandstone such as gritstone, clay, dense sandstone and gas sandstone and so on in the reservoir with traditional impedance inversion. According to the obvious difference between different lithologic and petrophysical impedance, we studied multi-scale reservoir heterogeneities and developed new technologies. The distribution features of reservoir lithological and petrophysical heterogeneities along vertical and transverse directions were described quantitatively using multi-scale power spectrum and heterogeneity spectrum methods in this paper. Power spectrum (P spectrum) describes the manner of the vertical distribution of reservoir lithologic and petrophysical parameters and the large-scale and small-scale heterogeneities along vertical direction. Heterogeneity spectrum (H spectrum) describes the structure of the reservoir lithologic and petrophysical parameters mainly, that is to say, proportional composition of each lithological and petrophysical heterogeneities are calculated in this formation. The method is more reasonable to describe the degree of transverse multi-scale heterogeneities in reservoir lithological and petrophysical parameters. Using information of sonic logs in Sulige oil field, two spectral methods have been applied to the oil field, and good analytic results have been obtained. In order to contrast the former researches, the last part is the multi-scale character analysis of reservoir based on the transmission character of wave using the wavelet transform. We discussed the method applied to demarcate sequence stratigraphy and also analyzed the reservoir interlayer heterogeneity.

二维自组织介质中地震波的传播特性

The topic of this study is simulation in the two dimensional self-organized media. The study in complexity of the earth plays an important role in structures, sources and energy seismic detection. And it mainly focuses on vertical or horizontal heterogeneous, anisotropic and linear media. Based on 2D self-organized velocity model and four-order finite-difference method, we simulate different types self-organized media and the same type mode with various parameters such as horizontal relative length, vertical relative length, variations, and velocity background gradient. Also we analyze the seismograms with complexity methods with instant information including amplitude, energy and frequency. The results can be summarized as the fallows: (1) The waveforms fluctuate with the velocity variations; (2) Different type self-organized media bring different effects on the amplitudes, energy and waveforms; (3) Different parameters also produce various influences to seismograms. (4) The layer contains their self-organized features, from which we can investigate the quality of the earth.

地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.

波动方程模拟与偏移的实用化方法研究

Along with the widespread and in-depth applications in petroleum prospecting and development, the seismic modeling and migration technologies are proposed with a higher requirement by oil industrial, and the related practical demand is getting more and more urgent. Based on theories of modeling and migration methods for wave equation, both related with velocity model, I thoroughly research and develop some methods for the goal of highly effective and practical in this dissertation. In the first part, this dissertation probes into the layout designing by wave equations modeling, focusing on the target-oriented layout designing method guided by wave equation modeling in complicated structure areas. It is implemented by using the fourth order staggered grid finite difference (FD) method in velocity-stress 2D acoustic wave equations plus perfectly matched layer (PML) absorbing boundary condition. To design target-oriented layout: (a) match the synthetic record on the surface with events of subsurface structures by analyzing the snapshots of theoretical model; (b) determine the shot-gather distance by tracking the events of target areas and measuring the receiving range when it reaches the surface; (c) restrict the range of valid shot-gather distance by drawing seismic windows in single shot records; (d) choose the best trace distance by comparing the resolution of prospecting targets from the simulated records with different trace distance. Eventually, we obtained the observation system parameters, which achieve the design requirements. In the second part, this dissertation presents the practical method to improve the 3D Fourier Finite Difference (FFD) migration, and carefully analyzes all the factors which influence 3D FFD migration’s efficiency. In which, one of the most important parameters of migration is the extrapolating step. This dissertation presents an efficient 3D FFD migration algorithm, which use FFD propagator to extrapolate wavefields over big layers, and use Born-Kirchhoff interpolator to image wavefields over small layers between the big ones. Finally, I show the effectiveness of this hybrid migration method by comparing migration results from 3D SEG/EAGE model with different methods.

喜马拉雅西构造结区岩石圈演化三维有限元动力学数值模拟

The past three decades have seen numerous attempts to numerically model stress and strain patterns in the lithosphere of the Earth on both global and regional scales. This efforts have been indispensable in identifying the features we need to include in our endeavour to develop better models of our planet’s lithosphere and they have also raised our awareness for the many unresolved issue in the deep geodynamical issues that need to be addressed in the future. Nonetheless, in most models, the lithosphere is treated as a single layer with depth-averaged properties, and as the same distribution in the stress and strain fields, and as deforming under plane strain. All these above make a great hander for its reality and degree of recognition. As the beginning in this paper, some principal numerical models and results on the evolution of Tibetan plateau are reviewed and analyzed. Then, the geological and geophysical expedition on the Western Himalayan Syntaxis is briefly reviewed. Furthermore, we analysis the feature in deep geophysical field studies in this area and adjacent regions. Because, for most continents, stress models driven by plate boundary forces have successfully reproduced the main characteristics of the stress and strain field, we present a set of three-dimensional models of lithosphere system for a simplified geometry of the Western Himalayan Syntaxis area and its adjacent regions, where we try to match the first-order characteristics of the stress and strain fields of lithosphere since 10 Ma, and deformation and geodynamical evolution process in former 2Ma. Of course, the kinematic boundary conditions of the stress models driven by plate boundary forces were applied. The rheology plays a significant role in the lithospheric tectonics, which lead to different rheological parameters were used in different works although the have the same constitutive equations in models. So, in this paper we do not aim to produce all characteristics of the Western Himalayan Syntaxis areas’ stress and strain fields by the choices of various parameters, but rather the dynamic response between various rheological parameters and stress and strain fields. We have chosen to concentrate on the importance of rheology and lateral strength variations for lithospheric stress and strain patterns and use our findings to build a model of the Western Himalayan Syntaxis areas. In doing so, we want to go beyond purely elastic models or purely viscoelastic models. Compared the results of the crust viscosity in the Western Himalayan Syntaxis areas, we believed that, when various viscoelastic models are adopted, the selection of the coefficient of viscosity in the Western Syntaxis area has important influence on the its uplifts and evolutions. A wider uplift ranges and gently elevation was observed at the same time when a lower viscosity was used in our models, and vice versa. Data of stress magnitudes are not available, but it is clear that the stress levels must be at or below the failure threshold of rock under compression. Under these criteria, the calculation results show that the viscosity in the Western Syntaxis area should be smaller than 1023Pa.s When elastic model is adopted in relatively rigid Tarim basin, obvious changes are induced to the stress and strain fields of the whole Western Syntaxis area. We found that rigid block of lithosphere reduced stress levels within its interior and that, at the edges of such regions, stress orientation can change. Furthermore there is no evidence that such rigid regions act as stress barriers in that they shield areas in opposite sides of the structure from the influence of one another. In our models, the upper crustal material of the Western Syntaxis area does not turns to move westward. Whereas, because of the stress and strain fields have been decoupling at the interior of the lithosphere, we can get the results that the deep material must not move westward.