EQUATION OF STATE FOR POLYMERS

In the theoretical study on equation of state for polymers, much attention has been paid to the polymer in liquid state, but less to that in solid state. Therefore, some empirical and semi-empirical equations of state have been used to describe its pressure-volume-temperature (P-V-T) relations.

Modification to the Hardisty equation, regarding the relationship between sediment transport rate and particle size

Bagnold-type bed-load equations are widely used for the determination of sediment transport rate in marine environments. The accuracy of these equations depends upon the definition of the coefficient k(1) in the equations, which is a function of particle size. Hardisty (1983) has attempted to establish the relationship between k(1) and particle size, but there is an error in his analytical result. Our reanalysis of the original flume data results in new formulae for the coefficient. Furthermore, we found that the k(1) values should be derived using u(1) and u(1cr) data; the use of the vertical mean velocity in flumes to replace u(1) will lead to considerably higher k(1) values and overestimation of sediment transport rates.

Low-amplitude BSRs and gas hydrate concentration on the northern margin of the South China Sea

The passive northern continental margin of the South China Sea is rich in gas hydrates, as inferred from the occurrence of bottom-simulating reflectors (BSR) and from well logging data at Ocean Drilling Program (ODP) drill sites. Nonetheless, BSRs on new 2D multichannel seismic reflection data from the area around the Dongsha Islands (the Dongsha Rise) are not ubiquitous. They are confined to complex diapiric structures and active fault zones located between the Dongsha Rise and the surrounding depressions, implying that gas hydrate occurrence is likewise limited to these areas. Most of the BSRs have low amplitude and are therefore not clearly recognizable. Acoustic impedance provides information on rock properties and has been used to estimate gas hydrate concentration. Gas hydrate-bearing sediments have acoustic impedance that is higher than that of the surrounding sediments devoid of hydrates. Based on well logging data, the relationship between acoustic impedance and porosity can be obtained by a linear regression, and the degree of gas hydrate saturation can be determined using Archie's equation. By applying these methods to multichannel seismic data and well logging data from the northern South China Sea, the gas hydrate concentration is found to be 3-25% of the pore space at ODP Site 1148 depending on sub-surface depth, and is estimated to be less than values of 5% estimated along seismic profile 0101. Our results suggest that saturation of gas hydrate in the northern South China Sea is higher than that estimated from well resistivity log data in the gas hydrate stability zone, but that free gas is scarce beneath this zone. It is probably the scarcity of free gas that is responsible for the low amplitudes of the BSRs.

Numerical wave flume study on wave motion around submerged plates

Nonlinear interaction between surface waves and a submerged horizontal plate is investigated in the absorbed numerical wave flume developed based on the volume of fluid (VOF) method. The governing equations of the numerical model are the continuity equation and the Reynolds-Averaged Navier-Stokes (RANS) equations with the k-epsilon turbulence equations. Incident waves are generated by an absorbing wave-maker that eliminates the waves reflected from structures. Results are obtained for a range of parameters, with consideration of the condition under which the reflection coefficient becomes maximal and the transmission coefficient minimal. Wave breaking over the plate, vortex shedding downwave, and pulsating flow below the plate are observed. Time-averaged hydrodynamic force reveals a negative drift force. All these characteristics provide a reference for construction of submerged plate breakwaters.

The nonlinear evolution equation on the shallow water wave

Based on the variation principle, the nonlinear evolution model for the shallow water waves is established. The research shows the Duffing equation can be introduced to the evolution model of water wave with time.

Asymptotic periodicity of the Volterra equation with infinite delay

For some species, hereditary factors have great effects on their population evolution, which can be described by the well-known Volterra model. A model developed is investigated in this article, considering the seasonal variation of the environment, where the diffusive effect of the population is also considered. The main approaches employed here are the upper-lower solution method and the monotone iteration technique. The results show that whether the species dies out or not depends on the relations among the birth rate, the death rate, the competition rate, the diffusivity and the hereditary effects. The evolution of the population may show asymptotic periodicity, provided a certain condition is satisfied for the above factors. (c) 2006 Elsevier Ltd. All rights reserved.

Asymptotic weighted-periodicity of the impulsive parabolic equation with time delay

The main aim of this paper is to investigate the effects of the impulse and time delay on a type of parabolic equations. In view of the characteristics of the equation, a particular iteration scheme is adopted. The results show that Under certain conditions on the coefficients of the equation and the impulse, the solution oscillates in a particular manner-called "asymptotic weighted-periodicity".

Traveling wave front for the Fisher equation on an infinite band region

Instead of discussing the existence of a one-dimensional traveling wave front solution which connects two constant steady states, the present work deals with the case connecting a constant and a nonhomogeneous steady state on an infinite band region. The corresponding model is the well-known Fisher equation with variational coefficient and Dirichlet boundary condition. (c) 2006 Elsevier Ltd. All rights reserved.

Roles of equatorial waves and western boundary reflection in the seasonal circulation of the equatorial Indian Ocean

An ocean general circulation model (OGCM) is used to study the roles of equatorial waves and western boundary reflection in the seasonal circulation of the equatorial Indian Ocean. The western boundary reflection is defined as the total Kelvin waves leaving the western boundary, which include the reflection of the equatorial Rossby waves as well as the effects of alongshore winds, off-equatorial Rossby waves, and nonlinear processes near the western boundary. The evaluation of the reflection is based on a wave decomposition of the OGCM results and experiments with linear models. It is found that the alongshore winds along the east coast of Africa and the Rossby waves in the off-equatorial areas contribute significantly to the annual harmonics of the equatorial Kelvin waves at the western boundary. The semiannual harmonics of the Kelvin waves, on the other hand, originate primarily from a linear reflection of the equatorial Rossby waves. The dynamics of a dominant annual oscillation of sea level coexisting with the dominant semiannual oscillations of surface zonal currents in the central equatorial Indian Ocean are investigated. These sea level and zonal current patterns are found to be closely related to the linear reflections of the semiannual harmonics at the meridional boundaries. Because of the reflections, the second baroclinic mode resonates with the semiannual wind forcing; that is, the semiannual zonal currents carried by the reflected waves enhance the wind-forced currents at the central basin. Because of the different behavior of the zonal current and sea level during the reflections, the semiannual sea levels of the directly forced and reflected waves cancel each other significantly at the central basin. In the meantime, the annual harmonic of the sea level remains large, producing a dominant annual oscillation of sea level in the central equatorial Indian Ocean. The linear reflection causes the semiannual harmonics of the incoming and reflected sea levels to enhance each other at the meridional boundaries. In addition, the weak annual harmonics of sea level in the western basin, resulting from a combined effect of the western boundary reflection and the equatorial zonal wind forcing, facilitate the dominance by the semiannual harmonics near the western boundary despite the strong local wind forcing at the annual period. The Rossby waves are found to have a much larger contribution to the observed equatorial semiannual oscillations of surface zonal currents than the Kelvin waves. The westward progressive reversal of seasonal surface zonal currents along the equator in the observations is primarily due to the Rossby wave propagation.

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.

随钻地震技术的研究及应用

Seismic While Drilling (SWD) is a new wellbore seismic technique. It uses the vibrations produced by a drill-bit while drilling as a downhole seismic energy source. The continuous signals generated by the drill bit are recorded by a pilot sensor attached to the top of the drill-string. Seismic wave receivers positioned in the earth near its surface receive the seismic waves both directly and reflection from the geologic formations. The pilot signal is cross-correlated with the receiver signals to compute travel-times of the arrivals (direct arrival and reflected arrival) and attenuate incoherent noise. No downhole intrusmentation is required to obtain the data and the data recording does not interfere with the drilling process. These characteristics offer a method by which borehole seismic data can be acquired, processed, and interpreted while drilling. As a Measure-While-Drill technique. SWD provides real-time seismic data for use at the well site . This can aid the engineer or driller by indicating the position of the drill-bit and providing a look at reflecting horizons yet to be encountered by the drill-bit. Furthermore, the ease with which surface receivers can be deployed makes multi-offset VSP economically feasible. First, this paper is theoretically studying drill-bit wavefield， interaction mode between drill-bit and formation below drill-bit , the new technique of modern signal process was applied to seismic data, the seismic body wave radiation pattern of a working roller-cone drill-bit can be characterized by theoretical modeling. Then , a systematical analysis about the drill-bit wave was done, time-distance equation of seismic wave traveling was established, the process of seismic while drilling was simulated using the computer software adaptive modeling of SWD was done . In order to spread this technique, I have made trial SWD modeling during drilling. the paper sketches out the procedure for trial SWD modeling during drilling , the involved instruments and their functions, and the trial effect. Subsurface condition ahead of the drill-bit can be predicted drillstring velocity was obtained by polit sensor autocorrelation. Reference decovolution, the drillstring multiples in the polit signal are removed by reference deconvolution, the crosscorrelation process enhance the signal-to-noise power ratio, lithologies. Final, SWD provides real-time seismic data for use at the well site well trajectory control exploratory well find out and preserve reservoirs. intervel velocity was computed by the traveltime The results of the interval velocity determination reflects the pore-pressure present in the subsurface units ahead of the drill-bit. the presences of fractures in subsurface formation was detected by shear wave. et al.

辽东湾隐蔽油气藏地球物理识别技术研究

In China and world, more than half the recent basin discovered reserves involve lithologic hydrocarbon reservoir reserves. The major target for further hydrocarbon basin exploration is the subtle reservoir. The Liaodong Bay prospect is much important in Bohai Sea, which includes Liaoxi low uplift, Liaodong uplift, Liaoxi sag and Liaozhong sag. After dozens years’ exploration in Liaodong Bay, few unexplored big-and-middle-sized favorable structural traps are remained and most of the stock structure targets are bad for fragmentary. Thus seeking for new prospect area and making a breakthrough, have become the unique way to relieve the severe exploration condition in Liaodong Bay. Technique Route Based on the petrophysical property of target area, the seismic forward inference of typical subtle trap model is expanded with analysis of logging, seismic and geologic data. According to petrophysical characteristics and forward inference and research on seismic response of actual seismic data in target area, the optimization of geophysical technique is used in subtle trap identification and the geophysical identification technique system of subtle reservoir is formed. The Key Research ① Petrophysical Model The petrophysical parameter is the basic parameter for seismic wave simulation. The seismic response difference of rocks bearing different fluids is required. With the crossplot of log data, the influence of petrophysical parameters on rock elastic properties of target area is analyzed, such as porosity, shale index, fluid property and saturation. Based on the current research on Biot-Gassmann and Kuster-Toksoz model, the petrophysical parameter calculator program which can be used for fluid substitution is established. ② S-wave evaluation based on conventional log data The shear velocity is needed during forward inference of AVO or other elastic wave field. But most of the recent conventional log data is lack of shear wave. Thus according to the research on petrophysical model, the rock S-wave parameter can be evaluated from conventional log data with probability inverse method. ③ AVO forward modeling based on well data For 6 wells in JZ31-6 block and 9 wells in LD22-1 block, the AVO forward modeling recording is made by log curve. The classification of AVO characteristics in objective interval is made by the lithologic information. ④ The 2D parameter model building and forward modeling of subtle hydrocarbon trap in target area. According to the formation interpretation of ESS03D seismic area, the 2D parameter model building and seismic wave field forward modeling are carried on the given and predicted subtle hydrocarbon trap with log curve. ⑤ The lithology and fluid identification of subtle trap in target area After study the seismic response characteristics of lithology and fluid in given target area, the optimization of geophysical technique is used for lithology identification and fluid forecast. ⑥The geophysical identification technique system of subtle reservoir The Innovative Points of this Paper ① Based on laboratory measurement and petrophysical model theory, the rock S-wave parameter can be evaluated from conventional log data with probability inverse method. Then the fluid substitution method based on B-G and K-T theory is provided. ② The method and workflow for simulating seismic wave field property of subtle hydrocarbon trap are established based on the petrophysical model building and forward modeling of wave equation. ③ The description of subtle trap structural feature is launched. According to the different reflection of frequency wave field structural attribute, the fluid property of subtle trap can be identified by wave field attenuation attribute and absorption analysis. ④ It’s the first time to identify subtle trap by geophysical technique and provide exploration drilling well location. ⑤ The technique system of subtle reservoir geophysical identification is formed to provide available workflow and research ideas for other region of interest.

基于弹性参数分析的碳酸盐岩储层预测方法研究

Prediction of Carbonate Reservoir Based on the Elastic Parameter Analysis Zhang Guangzhi (Solid Geophysics) Directed by Professor Liu Hong Abstract With the exploration and development of Puguang Oilfield, oil-gas exploration of carbonate rock in China has shown good prospects. Research on earthquake prediction methods for carbonate reservoir becomes the key of oil and gas exploration. Starting with analysis of geological characteristics of carbonate rock, prestack AVO inversion method, prestack elastic impedance inversion and parameter calculation method and seismic attribute extraction and optimization method were studied based on the analysis of rock physics in this work. First, variation characteristic and law of carbonate rock reservoir parameters were studied based on experimental data of rock physics, log data, analysis assay data, mud logging data and seismic data, so as to lay a foundation for the further reservoir identification and description. Then, the structure, type and propagation law of seismic wave field were analyzed through seismic forward modeling of the reservoir, and contact between information from log and geology data with elastic parameters, such as compressional wave and shear wave velocity and density were established, so as to provide a standard for reservoir identification and hydrocarbon detection using seismic reflection characteristics of the research area. Starting with the general concept of inverse problem, through analysis of Zoeppritz equation, three kinds of pre-stack inversion methods were derived and analyzed in detail, the AVO 3-parameter inversion based on Bayesian theory, the prestack AVO waveform inversion method and the simultaneous inversion method, based on the statistical hypothesis of inversion parameters and observation data and the Gauss distribution assumption of noise. The three methods were validated by model data and real data. Then, the elastic wave impedance inversion method of carbonate reservoir was investigated and the method of elastic parameter extraction from elastic impedance data was put forward. Based on the analysis of conventional methods of seismic attribute extraction and optimization, the time-frequency attributes and the wavelet attributes with time and amplitude feature were presented, and the prestack seismic attribute calculation method which can characterize the reservoir rock and fluid characteristic was presented. And the optimization of seismic attribute using the nonlinear KPCA method was also put forward. A series of seismic prediction technologies for carbonate reservoir were presented based on analysis of rock physics and seismic forward simulation technology. Practical application of these technologies was implemented in A oil field of Southern China and good effect has been achieved. Key words: carbonate rock; reservoir prediction; rock physics, prestack seismic inversion; seismic attribute

地震数据波动方程反演方法研究

Geophysical inversion is a theory that transforms the observation data into corresponding geophysical models. The goal of seismic inversion is not only wave velocity models, but also the fine structures and dynamic process of interior of the earth, expanding to more parameters such as density, aeolotropism, viscosity and so on. As is known to all, Inversion theory is divided to linear and non-linear inversion theories. In rencent 40 years linear inversion theory has formed into a complete and systematic theory and found extensive applications in practice. While there are still many urgent problems to be solved in non-linear inversion theory and practice. Based on wave equation, this dissertation has been mainly involved in the theoretical research of several non-linear inversion methods: waveform inversion, traveltime inversion and the joint inversion about two methods. The objective of gradient waveform inversion is to find a geologic model, thus synthetic seismograms generated by this geologic model are best fitted to observed seismograms. Contrasting with other inverse methods, waveform inversion uses all characteristics of waveform and has high resolution capacity. But waveform inversion is an interface by interface method. An artificial parameter limit should be provided in each inversion iteration. In addition, waveform information will tend to get stuck in local minima if the starting model is too far from the actual model. Based on velocity scanning in traditional seismic data processing, a layer-by-layer waveform inversion method is developed in this dissertation to deal with weaknesses of waveform inversion. Wave equation is used to calculate the traveltime and derivative (perturbation of traveltime with respect to velocity) in wave-equation traveltime inversion (WT). Unlike traditional ray-based travetime inversion, WT has many advantages. No ray tracing or traveltime picking and no high frequency assumption is necessary and good result can be got while starting model is far from real model. But, comparing with waveform inversion, WT has low resolution. Waveform inversion and WT have complementary advantages and similar algorithm, which proves that the joint inversion is a better inversion method. And another key point which this dissertation emphasizes is how to give fullest play to their complementary advantages on the premise of no increase of storage spaces and amount of calculation. Numerical tests are implemented to prove the feasibility of inversion methods mentioned above in this dissertation. Especially for gradient waveform inversion, field data are inversed. This field data are acquired by our group in Wali park and Shunyi district. Real data processing shows there are many problems for waveform inversion to deal with real data. The matching of synthetic seismograms with observed seismograms and noise cancellation are two primary problems. In conclusion, on the foundation of the former experiences, this dissertation has implemented waveform inversions on the basis of acoustic wave equation and elastic wave equation, traveltime inversion on the basis of acoustic wave equation and traditional combined waveform traveltime inversion. Besides the traditional analysis of inversion theory, there are two innovations: layer by layer inversion of seimic reflection data inversion and rapid method for acoustic wave-equation joint inversion.