STM image of tungstosilicic acid molecular layer on highly oriented pyrolytic graphite

Molecular layer of tungstosilicic acid (H4SiW12O40) deposited on freshly-cleaved highly oriented pyrolytic graphite (HOPG) was observed by scanning tunneling microscopy (STM) in air at room temperature. The molecular dimension (11.5 Angstrom) of H4SiW12O40 measured by STM is consistent with known crystallographic parameter. We also imaged the boundary of H4SiW12O40 molecular layer on HOPG showing that molecular layer of H4SiW12O40 was formed. It has been proved that individual tungstosilicic acid species is imaged. The probable reason for the formation of the molecular layer is also discussed.

IN-SITU SCANNING-TUNNELING-MICROSCOPY IMAGE OF A HIGHLY ORIENTED PYROLYTIC-GRAPHITE ELECTRODE MODIFIED WITH POTASSIUM BIS(11-MOLYBDOSILICATO) DYSPROSATE(III) IN SODIUM-SULFATE SOLUTION

The variation in molecule adsorption mode on pretreated highly oriented pyrolytic graphite electrodes, modified with the title complex K10H3[Dy(SiMo11O39)(2)] by cyclic voltammetry in the title complex solution, was observed in situ by electrochemical scanning tunnelling microscopy (ECSTM) with molecular resolution in sodium sulphate solution. According to the ECSTM images and the known molecular structure we conclude that the adsorption mode of the title complex modified electrode changed during potential cycling from ''vertical'' to ''inclined'' and then ''horizontal'' or ''flat'' mode, i.e. the title complex adsorbed on the surface of electrode by one ligand of the complex at first, then began to incline and was finally adsorbed by two ligands of the complex. This result indicates that the adsorption mode on the modified electrode surface changed during potential cycling in the sulphate solution and a much more stable molecular layer was formed. The change in adlattice of adsorbates on the modified electrode surface from hexagonal to rectangular was also observed by ECSTM. A plausible model was given to explain this process.

Influence of mean water depth and a subsurface sandbar on the onset and strength of wave breaking

Wave breaking in the open ocean and coastal zones remains an intriguing yet incompletely understood process, with a strong observed association with wave groups. Recent numerical study of the evolution of fully nonlinear, two-dimensional deep water wave groups identified a robust threshold of a diagnostic growth-rate parameter that separated nonlinear wave groups that evolved to breaking from those that evolved with recurrence. This paper investigates whether these deep water wave-breaking results apply more generally, particularly in finite-water-depth conditions. For unforced nonlinear wave groups in intermediate water depths over a flat bottom, it was found that the upper bound of the diagnostic growth-rate threshold parameter established for deep water wave groups is also applicable in intermediate water depths, given by k(0) h greater than or equal to 2, where k(0) is the mean carrier wavenumber and h is the mean depth. For breaking onset over an idealized circular arc sandbar located on an otherwise flat, intermediate-depth (k(0) h greater than or equal to 2) environment, the deep water breaking diagnostic growth rate was found to be applicable provided that the height of the sandbar is less than one-quarter of the ambient mean water depth. Thus, for this range of intermediate-depth conditions, these two classes of bottom topography modify only marginally the diagnostic growth rate found for deep water waves. However, when intermediate-depth wave groups ( k(0) h greater than or equal to 2) shoal over a sandbar whose height exceeds one-half of the ambient water depth, the waves can steepen significantly without breaking. In such cases, the breaking threshold level and the maximum of the diagnostic growth rate increase systematically with the height of the sandbar. Also, the dimensions and position of the sandbar influenced the evolution and breaking threshold of wave groups. For sufficiently high sandbars, the effects of bottom topography can induce additional nonlinearity into the wave field geometry and associated dynamics that modifies the otherwise robust deep water breaking-threshold results.

Joint statistical distribution of two-point sea surface elevations in finite water depth

Based on the second-order random wave solutions of water wave equations in finite water depth, a joint statistical distribution of two-point sea surface elevations is derived by using the characteristic function expansion method. It is found that the joint distribution depends on five parameters. These five parameters can all be determined by the water depth, the relative position of two points and the wave-number spectrum of ocean waves. As an illustrative example, for fully developed wind-generated sea, the parameters that appeared in the joint distribution are calculated for various wind speeds, water depths and relative positions of two points by using the Donelan and Pierson spectrum and the nonlinear effects of sea waves on the joint distribution are studied. (C) 2003 Elsevier B.V. All rights reserved.

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.

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

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

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.

Role of snow depth in spring of Tibetan Plateau in onset of South China Sea summer monsoon

The role of snow depth of Tibetan Plateau in the onset of South China Sea summer monsoon and the influence of ENSO on snow depth of Tibetan Plateau are investigated with use of data from ECMWF reanalysis and NCEP/NCAR reanalysis. The results are as follows: (1) The snow depth data from ECMWF reanalysis are tested and reliable, and can be used to study the influence of snow depth of Tibetan Plateau on the onset of South China Sea summer monsoon; (2) Anomaly of snow depth of Tibetan Plateau causes anomaly in air temperature and its contrast between the Indian Ocean and the continent resulting in easterly wind anomaly over 500 hPa and hence as well as in the atmospheric circulation in the lower layer. For the year of negative anomaly of snow depth a westerly wind anomaly with a cyclone pair takes place, while for positive anomaly of snow depth an easterly anomaly occurs with an anticyclone pair; (3) While positive anomaly of SST occurs in the eastern Pacific Ocean, positive anomaly of air pressure also takes place over the eastern Indian Ocean and the South China Sea, causing stronger meridional pressure gradient between the ocean and continent and then westerly wind anomaly. At the same time, the atmospheric pressure increases in the northern Tibetan Plateau, northerly wind gets stronger, and subtropical front strengthens. All of these are favorable for snowfall over Tibetan Plateau.

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.

Factors influencing the climatological mixed layer depth in the South China Sea: numerical simulations

The mixed layer depth (MLD) in the upper ocean is an important physical parameter for describing the upper ocean mixed layer. We analyzed several major factors influencing the climatological mixed layer depth (CMLD), and established a numerical simulation in the South China Sea (SCS) using the Regional Ocean Model System (ROMS) with a high-resolution (1/12A degrees x1/12A degrees) grid nesting method and 50 vertical layers. Several ideal numerical experiments were tested by modifying the existing sea surface boundary conditions. Especially, we analyzed the sensitivity of the results simulated for the CMLD with factors of sea surface wind stress (SSWS), sea surface net heat flux (SSNHF), and the difference between evaporation and precipitation (DEP). The result shows that of the three factors that change the depth of the CMLD, SSWS is in the first place, when ignoring the impact of SSWS, CMLD will change by 26% on average, and its effect is always to deepen the CMLD; the next comes SSNHF (13%) for deepening the CMLD in October to January and shallowing the CMLD in February to September; and the DEP comes in the third (only 2%). Moreover, we analyzed the temporal and spatial characteristics of CMLD and compared the simulation result with the ARGO observational data. The results indicate that ROMS is applicable for studying CMLD in the SCS area.

高阶广义屏波动方程叠前深度偏移

Attaining sufficient accuracy and efficiency of generalized screen propagator and improving the quality of input gathers are often problems　of wave equation presack depth migration,　in this paper,a high order formula of generalized screen propagator　for one-way wave equation is proposed by using the asymptotic expansion of single-square-root operator. Based on the formula，a new generalized screen propagator is developed ,which is composed of split-step Fourier propagator and high order correction terms，the new generalized screen propagator not only improving calculation precision without sharply increasing the quantity of computation，facilitates the suitability of generalized screen propagator to the media with strong lateral velocity variation.　As wave-equation prestack depth migration is sensitive to the quality of input gathers, which greatly affect the output，and the available seismic data processing system has inability to obtain traveltimes corresponding to the multiple arrivals, to estimate of great residual statics， to merge seismic datum from different projects and to design inverse Q filter, we establish　difference equations with an embodiment of Huygens’s principle for obtaining traveltimes corresponding to the multiple arrivals,bring forward a time variable matching filter for seismic datum merging by using the fast algorithm called Mallat tree for wavelet transformations， put forward a method for estimation of residual statics by applying the optimum model parameters estimated by iterative inversion with three organized algorithm,i.e,the CMP intertrace cross-correlation algorithm,the Laplacian image edge extraction algorithm,and the DFP algorithm, and present phase-shift inverse Q filter based on Futterman’s amplitude and phase-velocity dispersion formula and wave field extrapolation theory. All of their numerical and real data calculating results shows that our theory and method are practical and efficient. Key words: prestack depth migration, generalized screen propagator, residual statics，inverse Q filter ,traveltime，3D seismic datum mergence

VTI各向异性叠前深度偏移及速度分析

Theoretical research, laboratory test and field observation show that most of sediment rock has anisotropic features. It will produce some notable errors when applying isotropic methods such as prestack depth migration and velocity analysis to dada acquired under anisotropic condition; it also has a bad effect on geologic interpretation. Generally speaking, the vertical transverse isotropic media is a good approximation to geologic structure, thus it has an important realistic meaning for anisotropic prestack depth migration theory researching and precise complex geologic imaging if considering anisotropic effect of seismic wave propagation. There are two indispensable parts in prestack depth migration of realistic records, one is proper prestack depth migration algorithm, and the other is velocity analysis using prestack seismic data. The paper consists of the two aspects. Based on implicit finite difference research proposed by Dietrich Ristow et al （1997） about VTI media prestack depth migration, the paper proposed split-step Fourier prestack depth migration algorithm （VTISSF） and Fourier finite difference algorithm （VTIFFD） based on wave equation for VTI media, program are designed and the depth migration method are tested using synthetic model. The result shows that VTISSF is a stable algorithm, it generally gets a good result if the reflector dip is not very steep, while undermigration phenomena appeared in steep dips case; the VTIFFD algorithm bring us better result in steep dips with lower efficiency and frequency dispersion. For anisotropic prestack depth migration velocity analysis of VTI media, The paper discussed the basic hypothesis of VTI model in velocity analysis algorithm, basis of anisotropic prestack depth migration velocity analysis and travel time table calculation of VTI media in integral prestack depth migration. Then , analyzed the P-wave common imaging gather in the case of homogeneous velocity and vertically variable velocity . studied the residual correction in common imaging gather produced by media parameter error, analyzed the condition of flat event and correct depth in common imaging gather . In this case, the anisotropic model parameter vector is , is vertical velocity of a point at top surface, is vertical velocity gradient, and are anisotropic parameter. We can get vertical velocity gradient from seismic data; then the P-wave common imaging gather of VTI media whose velocity varies in vertical and horizontal direction, the relationship between media parameter and event residual time shift of common image gather are studied. We got the condition of flattening common imaging gather with correct depth. In this case the anisotropic model parameter vector is , is velocity gradient in horizontal direction. As a result, the vertical velocity grads can be decided uniquely, but horizontal velocity grads and anisotropic parameter can’t be distinguished if no priori information available, our method is to supply parameter by velocity scanning; then, as soon as is supplied we can get another four parameters of VTI media from seismic data. Based on above analysis, the paper discussed the feasibility of migration velocity analysis in vertically and horizontally varied VTI media, synthetic record of three models are used to test the velocity analysis method . Firstly, anisotropic velocity analysis test is done using a simple model with one block, then we used a model with multiple blocks, thirdly, we analyzed the anisotropic velocity using a part of Marmousi model. The model results show that this velocity analysis method is feasible and correct.

构造复杂区变速成图方法研究

In the complex structure areas, velocity field building and structure mapping are important for seismic exploration. With the development of seismic exploration, the methods of structure mapping, reservoir prediction and reservoir description all require high precious velocity field. And more accurate depth-structure maps are required for well site design. Aiming at the problems and defects in velocity analysis and structure mapping in oil seismic exploration, the paper which is based on the studies of real data in several areas combines the theories with practical application, and analyzes the precision and applicability of several methods of velocity model building. After that, the following methods are mainly studied: the coherence inversion methods based on the pre-stack CMP gathers or stacking velocity; the interval velocity inversion methods constrained by multi-well; the Random Simulation method; 3D Image Ray Map Migration method and the structure mapping in floating datum and in fixed datum, and then we conclude the method of building high precious seismic velocity field and structure mapping with variable velocity. Firstly, the paper analyses the distributing rule of the velocity variation in the areas with complex structures in the northwest of China, then points out that velocity is a crucial factor which influences the precision of structure mapping, and the velocity variations have something to do with the shapes of the structures, the variety of lithology and so on. The key point of improving the precision of seismic velocity field is to obtain a structure mapping with high precision. We also describe the range and conditions of these methods. Secondly, by comparing many popular methods of velocity model building, we propose a new method in the use of velocity model building. The new method is more effective in velocity model building under every kind of complex condition and is worthy of spreading. At last, the paper fingers out that it is a system engineering to study variable velocity mapping in every kind of complex structure areas. Every step of the work can affect the final results. So it is important to build high efficient and practical velocity model and the flows of mapping processing. The paper builds the flows and gives some examples. The method has been applied in more than ten exploring surveys. The application proves that this method could bring good effect on researching on low-amplitude trap, reservoir prediction, reservoir description and the integrated research of oil&gas geology. Keywords: structure mapping velocity model building complex structure variable velocity media

裂隙/孔隙型弹性介质中地震波场模拟与叠前逆时偏移

With the development of oil/gas seismic exploration, seismic survey for fracture/porosity type reservoir is becoming more and more important. As for China, since it has over 60% store of low porosity and low permeability oil/gas reservoir, it’s more urgent to validly describe fracture/porosity type oil/gas trap and proposing the related, developed seismic technique. To achieve mapping fracture/porosity region and its development status, it demands profound understanding of seismic wave propagation discipline in complex fractured/pored media. Meanwhile, it has profound scientific significance and applied worth to study forward modeling of fracture/porosity type media and pre-stacked reverse time migration. Especially, pre-stacked reverse-time migration is the lead edge technique in the field of seismology and seismic exploration. In this paper, the author has summarized the meaning, history and the present state of numerical simulation of seismic propagation in fractured/pored media and seismic exploration of fractured/pored reservoirs. Extensive Dilatancy Anisotropy (EDA) model is selected as media object in this work. As to forward modeling, due to local limitation of solving spatial partial derivative when using finite-difference and finite-element method, the author turns to pseudo-spectral method (PSM), which is based on the global characteristic of Fourier transform to simulate three-component elastic wave-field. Artifact boundary effect reduction and simulation algorithm stability are also discussed in the work. The author has completed successfully forward modeling coding of elastic wave-field and numerical simulation of two-dimensional and three-dimensional EDA models with different symmetric axis. Seismic dynamic and kinematical properties of EDA media are analyzed from time slices and seismic records of wave propagation. As to pre-stacked reverse-time migration for elastic wave-field in fractured/pored media, based on the successful experience in forward modeling results with PSM, the author has studied pre-stacked reverse-time depth-domain migration technique using PSM of elastic wave-field in two dimensional EDA media induced by preferred fracture/pore distribution. At the same time, different image conditions will bring up what kind of migration result is detailed in this paper. The author has worded out software for pre-stacked reverse-time depth-domain migration of elastic wave-field in EDA media. After migration processing of a series of seismic shot gathers, influences to migration from different isotropic and anisotropy models are described in the paper. In summary, following creative research achievements are obtained:  Realizing two-dimensional and three-dimensional elastic wave-field modeling for fractured/pored media and related software has been completed.  Proposed pre-stacked reverse-time depth-domain migration technique using PSM of elastic wave-field.  Through analysis of the seismic dynamic and kinematical properties of EDA media, the author made a conclusion that collection of multi-component seismic data can provide important data basis for locating and describing the fracture/pore regions and their magnitudes and the preferred directions.  Pre-stacked reverse-time depth-domain migration technique has the ability to reconstruct complex geological object with steep formations and tilt fracture distribution. Neglecting seismic anisotropy induced by the preferred fracture/pore distribution, will lead to the disastrous imaging results.