THE ALGEBRAIC-METHOD ESTIMATING +/-PI/2 TYPE OF PHASES OF STRUCTURE FACTORS

The algebraic formulas of 1.5 and 2.5 rank which can be applied to estimating +/- pi/2 type of phases for P2(1)2(1)2(1) space group were derived using the method of structure factor algebra. Both types of the formulas are satisfactory for two known crystal structures in estimating their +/- pi/2 type of phases.

Determination of fractional energy loss of waves in nearshore waters using an improved high-order Boussinesq-type model

Fractional energy losses of waves due to wave breaking when passing over a submerged bar are studied systematically using a modified numerical code that is based on the high-order Boussinesq-type equations. The model is first tested by the additional experimental data, and the model's capability of simulating the wave transformation over both gentle slope and steep slope is demonstrated. Then, the model's breaking index is replaced and tested. The new breaking index, which is optimized from the several breaking indices, is not sensitive to the spatial grid length and includes the bottom slopes. Numerical tests show that the modified model with the new breaking index is more stable and efficient for the shallow-water wave breaking. Finally, the modified model is used to study the fractional energy losses for the regular waves propagating and breaking over a submerged bar. Our results have revealed that how the nonlinearity and the dispersion of the incident waves as well as the dimensionless bar height (normalized by water depth) dominate the fractional energy losses. It is also found that the bar slope (limited to gentle slopes that less than 1:10) and the dimensionless bar length (normalized by incident wave length) have negligible effects on the fractional energy losses.

Improvement of different source function expressions in SWAN model for the Bohai Sea

The simulating wave nearshore (SWAN) wave model has been widely used in coastal areas, lakes and estuaries. However, we found a poor agreement between modeling results and measurements in analyzing the chosen four typical cases when we used the default parameters of the source function formulas of the SWAN to make wave simulation for the Bohai Sea. Also, it was found that at the same wind process the simulated results of two wind generation expressions (Komen, Janssen) demonstrated a large difference. Further study showed that the proportionality coefficient alpha in linear growth term of wave growth source function plays an unperceived role in the process of wave development. Based on experiments and analysis, we thought that the coefficient alpha should change rather than be a constant. Therefore, the coefficient alpha changing with the variation of friction velocity U (*) was introduced into the linear growth term of wave growth source function. Four weather processes were adopted to validate the improvement in the linear growth term. The results from the improved coefficient alpha agree much better with the measurements than those from the default constant coefficient alpha. Furthermore, the large differences of results between Komen wind generation expression and Janssen wind generation expression were eliminated. We also experimented with the four weather processes to test the new white-capping mechanisms based on the cumulative steepness method. It was found that the parameters of the new white-capping mechanisms are not suitable for the Bohai Sea, but Alkyon's white-capping mechanisms can be applicable to the Bohai Sea after amendments, demonstrating that this improvement of parameter alpha can improve the simulated results of the Bohai Sea.

Effective dielectric response of nonlinear composites with external ac and dc electric field

The perturbation method is developed to investigate the effective nonlinear dielectric response of Kerr composites when the external ac and dc electric field is applied. Under the external ac and dc electric field E-app=E-a(1+sin omegat), the effective coupling nonlinear response can be induced by the cubic nonlinearity of Kerr nonlinear materials at the zero frequency, the finite basic frequency omega, the second and the third harmonics, 2omega and 3omega, and so on. As an example, we have investigated the cylindrical inclusions randomly embedded in a host and derived the formulas of the effective nonlinear dielectric response at harmonics in dilute limit. For a higher concentration of inclusions, we have proposed a nonlinear effective-medium approximation by introducing the general effective nonlinear response. With the relationships between the effective nonlinear response at harmonics and the general effective nonlinear response, we have derived a set of formulas of the effective nonlinear dielectric responses at harmonics for a larger volume fraction. (C) 2004 American Institute of Physics.

Effective dielectric response of composites witli graded material

The effective dielectric response of linear composites containing graded material is investigated under an applied electric field Eo. For the cylindrical inclusion with gradient dielectric function, epsilon(i)(r) = b + cr, randomly embedded in a host with dielectric constant epsilon(m), we have obtained the exact solution of local electric potential of the composite media regions, which obeys a linear constitutive relation D = epsilonE, using hypergeometric function. In dilute limit, we have derived the effective dielectric response of the linear composite media. Furthermore, for larger volume fraction, the formulas of effective dielectric response of the graded composite media, are given.

Effective properties of graded elliptical cylindrical composites

The effective property has been investigated theoretically in graded elliptical cylindrical composite's consisting of inhomogeneous graded elliptical cylinders and an isotropic matrix under external uniform electric field. As a theoretical model, the dielectric gradient profile in the elliptical cylinder is modeled by a power-law function of short semi-axis variable parameter (xi(2) - 1) in the elliptical cylindrical coordinates, namely epsilon(i)(xi) = c(k) (xi(2) - 1)(k), where c(k) and k are the parameters, and xi is the long semi-axis space variable in an elliptical cylindrical inclusion region. In the dilute limit, the local analytical potentials in inclusion and matrix regions are derived exactly by means of the hyper-geometric function, and the formulas are given for estimating the effective dielectric responses under the external lfield along (x) over cap- and (y) over cap -directions, respectively. Furthermore, we have demonstrated that our effective response formulas can be reduced to the well-known results of homogeneous isotropic elliptical cylindrical composites if we take the limit k -> 0 in graded elliptical cylindrical composites. (c) 2006 Elsevier B.V. All rights reserved.

Effective properties of spherically anisotropic piezoelectric composites

Piezoelectric composites consisting of spherically anisotropic piezoelectric inclusions (i.e., piezoceramic material) in an infinite nonpiezoelectric matrix under a uniform electric field are theoretically investigated. Analytical solutions for the elastic displacements and the electric potentials are derived exactly. Taking account of the coupling effects of elasticity, permittivity, and piezoelectricity, formulas are derived for the effective dielectric and piezoelectric responses in the dilute limit. A piezoelectric response mechanism is revealed, in which the effective piezoelectric response vanishes irrespective of how much spherically anisotropic piezoelectric inclusions are inside. Moreover, the effective coupled responses of the piezoelectric composites show that the effective dielectric responses decrease (increase) as the inclusion elastic (piezoelectric) constants increase.

Dielectric response of spherically anisotropic graded piezoelectric composites

A graded piezoelectric composite consisting of a spherically anisotropic graded piezoelectric inclusion imbedded in an infinite nonpiezoelectric matrix, with the physical properties of the graded spherical inclusion having a power-law profile with respect to the radial variable r, is studied theoretically. Under an external uniform electric field, the electric displacement field and the elastic stress tensor field of this spherically anisotropic graded piezoelectric composite are derived exactly by means of displacement separation technique, based on the governing equations in the dilute limit. A piezoelectric response mechanism, in which the effective piezoelectric response vanishes along the z direction (or x,y directions), is revealed in this kind of graded piezoelectric composites. Furthermore, it is found that the effective dielectric constant decreases (or increases) with the volume fraction p of the inclusions if the exponent parameter k of the grading profile is larger (or smaller) than a critical value. (C) 2007 American Institute of Physics.

基于AFM的机器人化纳米操作中纳观力的初步研究

针对基于原子力显微镜(AFM)的机器人化纳米操作,对探针作用下探针—基片—微粒之间纳观力的作用规律进行了初步分析.指出起主要作用的纳观力为范德华力、接触斥力、纳米摩擦力、毛细作用力以及纳米静电力等五种,并初步推导了各种纳观力的表达形式.通过力—距离曲线仿真与实验验证了所进行分析的合理性;该分析有助于进行纳米操作的精确控制.

原子力显微镜管式扫描器运动学建模与误差分析

针对样本扫描模式原子力显微镜,对其管式扫描器-样本-探针系统进行了运动学分析,建立了该系统的运动学模型,该模型表明:对于给定原子力显微镜扫描器,样本上与探针接触点的横向和纵向位移取决于探针尖端相对于扫描管轴心的偏置量、所加电压(或名义扫描范围)及样本厚度。据此模型,对由于弯曲运动模式所产生的两种重要误差—交叉耦合误差及扫描范围误差进行了定量分析,分析表明:扫描范围误差主要受样本厚度及名义扫描范围影响,而Z向交叉耦合误差主要受探针偏置量及名义扫描范围影响,实验验证了所建立的运动学模型和误差计算公式的正确性;另外,还提出了相应的减小误差的方法。

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

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.

旋转轴对称介质中的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

扩散波勘探方法基础及其在石油油气藏开发中的应用研究

As we know, the essence of exploration is objective body determined by getting the information. Such as seismic、electrical and electromagnetic prospecting, they are the common methods of the exploration. Therefore, They have a complete set of theory now. In fact, the effective information can also be got by the diffusion way, it is called diffusion prospecting. The diffusion way prospecting is necessary and important. The way of diffusion prospecting is studied in the paper and main works include below: (1) On the basis of studying basic law of the diffusion, the paper gives the idea of diffusion wave and the formulas of computing diffusion wave function. (2) The paper studies the way of the diffusion prospecting and the methods of data processing. At the same time, it also expounds the characteristics and the applied foreground of the diffusion prospecting. (3) The paper gives the tomography idea and the basic method of diffusion CT. Meanwhile, it also expounds the foreground that the diffusion CT is applied in oil development prospecting. (4) As the inversion of the diffusion equation is a part of the diffusion prospecting way, the methods of diffusion equation inversion are studied and the two formulas are deduced --Laplace transform and polynomial fitting inversion formulas. As the other important result of diffusion equation inversion, the inversion can offer a new analysis method for well Testing in oil development. In order to show a set of methods in the paper feasible, forward、inversion and CT numerical simulation are done in the paper.

双聚焦偏移速度估计与可分表示叠前弹性反演

In exploration seismology, the geologic target of oil and gas reservoir in complex medium request the high accuracy image of the structure and lithology of the medium. So the study of the prestack image and the elastic inversion of seismic wave in the complex medium come to the leading edge. The seismic response measured at the surface carries two fundamental pieces of information: the propagation effects of the medium and the reflections from the different layer boundaries in the medium. The propagation represent the low-wavenumber component of the medium, it is so-called the trend or macro layering, whereas the reflections represent the high-wavenumber component of the medium, it is called the detailed or fine layering. The result of migration velocity analysis is the resolution of the low-wavenumber component of the medium, but the prestack elastic inversion provided the resolution of the high-wavvenumber component the medium. In the dissertation, the two aspects about the migration velocity estimation and the elastic inversion have been studied.Firstly, any migration velocity analysis methods must include two basic elements: the criterion that tell us how to know whether the model parameters are correct and the updating that tell us how to update the model parameters when they are incorrect, which are effected on the properties and efficiency of the velocity estimation method. In the dissertation, a migration velocity analysis method based on the CFP technology has been presented in which the strategy of the top-down layer stripping approach are adapted to avoid the difficult of the selecting reduce .The proposed method has a advantage that the travel time errors obtained from the DTS panel are defined directly in time which is the difference with the method based on common image gather in which the residual curvature measured in depth should be converted to travel time errors.In the proposed migration velocity analysis method, the four aspects have been improved as follow:? The new parameterization of velocity model is provided in which the boundaries of layers are interpolated with the cubic spline of the control location and the velocity with a layer may change along with lateral position but the value is calculated as a segmented linear function of the velocity of the lateral control points. The proposed parameterization is suitable to updating procedure.? The analytical formulas to represent the travel time errors and the model parameters updates in the t-p domain are derived under local lateral homogeneous. The velocity estimations are iteratively computed as parametric inversion. The zero differential time shift in the DTS panel for each layer show the convergence of the velocity estimation.? The method of building initial model using the priori information is provided to improve the efficiency of velocity analysis. In the proposed method, Picking interesting events in the stacked section to define the boundaries of the layers and the results of conventional velocity analysis are used to define the velocity value of the layers? An interactive integrate software environment with the migration velocity analysis and prestack migration is built.The proposed method is firstly used to the synthetic data. The results of velocity estimation show both properties and efficiency of the velocity estimation are very good.The proposed method is also used to the field data which is the marine data set. In this example, the prestack and poststack depth migration of the data are completed using the different velocity models built with different method. The comparison between them shows that the model from the proposed method is better and improves obviously the quality of migration.In terms of the theoretical method of expressing a multi-variable function by products of single-variable functions which is suggested by Song Jian (2001), the separable expression of one-way wave operator has been studied. A optimization approximation with separable expression of the one-way wave operator is presented which easily deal with the lateral change of velocity in space and wave number domain respectively and has good approach accuracy. A new prestack depth migration algorithm based on the optimization approximation separable expression is developed and used to testing the results of velocity estimation.Secondly, according to the theory of the seismic wave reflection and transmission, the change of the amplitude via the incident angle is related to the elasticity of medium in the subsurface two-side. In the conventional inversion with poststack datum, only the information of the reflection operator at the zero incident angles can be used. If the more robust resolutions are requested, the amplitudes of all incident angles should be used.A natural separable expression of the reflection/transmission operator is represented, which is the sum of the products of two group functions. One group function vary with phase space whereas other group function is related to elastic parameters of the medium and geological structure.By employing the natural separable expression of the reflection/transmission operator, the method of seismic wave modeling with the one-way wave equation is developed to model the primary reflected waves, it is adapt to a certain extent heterogeneous media and confirms the accuracy of AVA of the reflections when the incident angle is less than 45'. The computational efficiency of the scheme is greatly high.The natural separable expression of the reflection/transmission operator is also used to construct prestack elastic inversion algorithm. Being different from the AVO analysis and inversion in which the angle gathers formed during the prstack migration are used, the proposed algorithm construct a linear equations during the prestack migration by the separable expression of the reflection/transmission operator. The unknowns of the linear equations are related to the elasticity of the medium, so the resolutions of them provided the elastic information of the medium.The proposed method of inversion is the same as AVO inversion in , the difference between them is only the method processing the amplitude via the incident angle and computational domain.

句容海安区块热史、成烃史与二次生烃

Jurong -Hai'an block of lower Yangtz area is one of the important petroleum exploration area among the residual marine basins in the south China. In the history of the basin's evolution, the strongly compressing, napping, folding and deforming during Indosinian to Yanshan epoch resulted in destruction of the early formed petroleum pool. Therefore, the strategy exploration of the secondary hydrocarbon generation and later formation of petroleum pool was brought forward for petroleum exploration and planning in the study area..On the basis of tectonic and sedimentary as well as present-day thermal regime evolution, using vitrinite reflectance and apatite fission track data and the paleo-heat flow based method, the basin's thermal history is reconstructed and hence the strata's temperature history are obtained. In addition, the maturation histories of the main four sets of marine hydrocarbon source rocks in the block are calculated. Furthermore, taking the maturity in the end of early Cretaceous as the original maturity and according to the formulas fitted by the secondary hydrocarbon generation model, the secondary hydrocarbon generation potential of the four sets of source rock is evaluated.The results of thermal history reconstruction show that Jurong-Hai'an block was under an uniform thermal setting during the Caledonian to Hercynian period and characterized by middle heat flow (52~57rnW/ m2). The uniform thermal setting was divided during and after Indosinian to Yanshan epoch. Wuwei area of southern Anhui province was under the high heat flow setting(~90 mW/m2)between 236Ma and 137Ma; Jurong area of southern Jiangsu ,Huangqiao area and Subei basin reached its maximum heat flow of 90,84 and 78-82 mW/m2 at 101Ma,157Ma and 56Ma respectivelyThe study of secondary hydrocarbon generation shows that the upper Paleozoic and Triassic source rocks have excellent secondary hydrocarbon generation potential. Silurian source rock posses some secondary hydrocarbon generation. Cambrian source rock, however, nearly has no secondary hydrocarbon generation. Overall there are no advantages of secondary hydrocarbon generation in the southern area of Jiansu. The intensity of secondary hydrocarbon generation in North Jiangsu basin is definitely better than that in the southern area of Jiangsu.