断阶型潜山带地质模型的建立及储层预测--以CG20潜山为例

以 CG2 0潜山为例 ,从建立地质模型入手 ,包括地层模型、构造模型、储集模型、储盖组合模型、速度模型等 ,认识到各套地层分布和储层物性的差异均与地震响应密切相关 ,因此可以利用地震波的信息 ,预测潜山储层的发育及分布情况。在对 CG2 0潜山进行精细全三维构造解释的基础上 ,探讨性地应用了测井约束反演、吸收系数、相干分析及三维模式识别等技术 ,对潜山储层进行了预测 ,从而提高了潜山勘探的效益 ,并为类似断阶型潜山带的勘探提供了成功的经验 ,具有一定的指导意义。     

A Reconstruction of Observed Profiles in the Sea East of Japan Using Vertical Coupled Temperature-Salinity EOF Modes

It is important to estimate hard-to-observe parameters in the ocean interior from easy-to-observe parameters. This study therefore demostrates a reconstruction of observed temperature and salinity profiles of the sea east of Japan (30°≈40°N, 140°≈150°E). The reconstruction was done by estimating suboptimal state from several values of the observed profiles and/or sea surface dynamic height (SDH) calculated from the profiles. The estimation used a variational method with vertical coupled temperature-salinity empirical orthogonal function (EOF) modes. Profiles of temperature and salinity in the subtropical region are effectively reconstructed from in situ temperature profile data, or sea surface temperature (SST) and SDH. For example, the analyzed temperature field from SST and SDH has an accuracy to within 1°C in the subtropical region. Salinity in the sea north of Kuroshio, however, is difficult to estimate because of its complex variability which is less correlated with temperature than in the subtropical region. Sea surface salinity is useful to estimate the subsurface structure. We also show the possibility that the estimation is improved by considering nonlinearity in the equation calculating SDH from temperature and salinity analysis values in order to examine the misfit between analysis and observation. Analysis using TOPEX/POSEIDON altimetry data instead of SDH was also performed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

P矢量方法在南海夏季环流诊断计算中的应用

基于1998年6～7月南海调查航次的CTD资料,对南海环流采用最近发展的P矢量方法进行诊断计算.计算结果:黑潮向西入侵南海,然后做反气旋弯曲向东北方向流动,最终有通过巴士海峡流出南海的趋势.在南海北部存在一个气旋性环流,这个环流的强度和范围随深度增加而减小.该环流的冷中心位置随深度增加稍向南移.南海中部、越南以东海域存在一个明显的气旋涡和反气旋涡,尤其在200m及其以上水层均相当稳定,反气旋涡位于越南以东,其中心位置在11°53'N,111°50'E,气旋涡的中心位置在13°17'N,112°55'E,两者的尺度皆约为250km.吕宋岛西侧存在一个反气旋涡.在计算海区南部、巴拉望岛西南海域,100m以上层存在一个反气旋式涡.从各层流场分布均可以显示海流在西部强化的现象.     

Modified mild slope equation and open boundary conditions

A numerical model to compute wave field is developed. It is based on the Berkhoff diffraction-refraction equation, in which an energy dissipation term is added, to take into account the breaking and the bottom friction phenomena. The energy dissipation function, by breaking and by bottom friction, is introduced in the Berkhoff equation to obtain a new equation of propagation.The resolution is done with the hybrid finite element method, where lagrangians elements are used.     

Analysis of seabed instability using element free Galerkin method

Wave-induced seabed instability, either momentary liquefaction or shear failure, is an important topic in ocean and coastal engineering. Many factors, such as seabed properties and wave parameters, affect the seabed instability. A non-dimensional parameter is proposed in this paper to evaluate the occurrence of momentary liquefaction. This parameter includes the properties of the soil and the wave. The determination of the wave-induced liquefaction depth is also suggested based on this non-dimensional parameter. As an example, a two-dimensional seabed with finite thickness is numerically treated with the EFGM meshless method developed early for wave-induced seabed responses. Parametric study is carried out to investigate the effect of wavelength, compressibility of pore fluid, permeability and stiffness of porous media, and variable stiffness with depth on the seabed response with three criteria for liquefaction. It is found that this non-dimensional parameter is a good index for identifying the momentary liquefaction qualitatively, and the criterion of liquefaction with seepage force can be used to predict the deepest liquefaction depth.     

Facies recognition of some Tertiary coals applied to prediction of oil source rock occurrence

Coals are oil source rocks in many of the Tertiary basins of Southeast Asia. The precursors of these hydrogen rich and oxygen poor coals are coastal plain peats which have mainly developed in an everwet and tropical climate. In these environments water flow and reworking can concentrate liptinitic kerogen in preference to vitrinitic kerogen. The distribution, petrography and chemistry of the coaly Miocene source rocks present in the Kutai Basin are described. The recognition of environmental controls on the accumulation of potentially oil-prone coals and coaly shales in deltaic environments is an aid to predictive source bed recognition in petroleum exploration. Comments on the environment of deposition of coaly sediments in the basins of the Norwegian Sea are discussed with reference to their possible oil and/or gas sourcing potential. The Triassic - Jurassic coals of the Haltenbanken area may become more oil-prone towards the delta margins, and facies mapping could aid oil exploration in this area.     

Baroclinic Buoyancy-Inertia Joint Stability Parameter

At present, the barotropic buoyant stability parameter has been derived from a vertical virtual displacement of a water parcel. The barotropic inertial stability parameter in the eccentrically cyclogeostrophic, basic current field was derived in 2003 from a horizontal cross-stream virtual displacement of a parcel. By expressing acceleration of a parcel due to a virtual displacement, which is arbitrarily sloping within a vertical section across the basic current, in terms of natural coordinates, we derived the vertical component of baroclinic buoyant stability parameter B ₂ ², the horizontal component of baroclinic inertial stability parameter I ₂ ², the baroclinic joint stability parameter J ², its buoyant component B ² and its inertial component I ². B ² is far greater than I ₂ ², and when neglecting relative vorticity except for vertical shear, a downward convex curve of J ² plotted against the slope of a virtual displacement follows a trend of B ² curve. If a parcel displaces along a horizontal surface or an isopycnal surface, however, B ² vanishes, and J ² becomes equal to I ². Actual parcel is apt to displace not only along the bottom slope, but also along the sea surface and an isopycnal interfacial surface, which is approximately equivalent to an isentropic surface, preferred by lateral mixing and exchange of momentum. Such actual displacement makes B ² vanishing, and grants I ² an important role. The present analysis of I ² examining effects due to curvature and horizontal and vertical shear vorticities are useful in deepening our understanding of baroclinic instability in actual oceanic streams.     

Preconditioned Optimizing Utility for Large-dimensional analyses (POpULar)

I present the derivation of the Preconditioned Optimizing Utility for Large-dimensional analyses (POpULar), which is developed for adopting a non-diagonal background error covariance matrix in nonlinear variational analyses (i.e., analyses employing a non-quadratic cost function). POpULar is based on the idea of a linear preconditioned conjugate gradient method widely adopted in ocean data assimilation systems. POpULar uses the background error covariance matrix as a preconditioner without any decomposition of the matrix. This preconditioning accelerates the convergence. Moreover, the inverse of the matrix is not required. POpULar therefore allows us easily to handle the correlations among deviations of control variables (i.e., the variables which will be analyzed) from their background in nonlinear problems. In order to demonstrate the usefulness of POpULar, we illustrate two effects which are often neglected in studies of ocean data assimilation before. One is the effect of correlations among the deviations of control variables in an adjoint analysis. The other is the nonlinear effect of sea surface dynamic height calculation required when sea surface height observation is employed in a three-dimensional ocean analysis. As the results, these effects are not so small to neglect.