Fluid flow equations governing primary oil migration as a separate phase

A mathematical model of primary oil migration as a separate phase out of compacting shales is presented. During burial and oil generation, source rock porosity decreases and oil saturation increases until residual oil saturation is reached. At this stage oil is expelled out by capillary and excess fluid pressure gradients. The model is a system of differential equations which relate changes in oil and water saturation in time to water and oil flow out of the source rock during burial. An additional set of equations for periods of erosion of overburden are also provided. The equations can be numerically solved by finite difference method. If oil and water flow is to be simulated during oil generation, then at each time step, changes by oil generation in oil and water saturations and porosity must be calculated. The solution procedure is briefly outlined.     

Simulation of green water flow on deck using non-linear dam breaking theory

Experiments carried out with models of floating production, storage and offloading platforms (FPSOs) showed that the flow of water over the deck edge, onto the deck resembled a suddenly released wall of water rather than a breaking wave. Therefore green water flow onto the deck was simulated using dam breaking theory, but the theory’s shallow-water assumptions may be limiting. In this paper a non-linear dam breaking problem is formulated. Equations of motion in the Lagrangian form are used and the solution is sought as an infinite series in time. Comparisons with the shallow water approximation are carried out.     

A Computational Model for Velocity Separation in Shallow Sea

Based on the hydrodynamical feature and the theoretical velocity profiles of tidal flow and vvind-induced flow in shal-low sea, a computational model is established for the first time, which can separate observed velocity into tidal velocity and wind-induced velocity by use of the least square method. With the model, not only the surface velocities of tidal flow and vvind-induced flow are obtained, but also the bed roughness height is determined and the wind velocity above the wa-ter surface is estimated. For verification of the model, the observed velocity in the Yellow River Estuary and the laborato-ry test is separated, then it is applied to the Yangtze River Estuary. All the results are satisfactory. The research results show that the model is simple in method, feasible in process and reasonable in result. The model is a valid approach to analysis and computation of field dala, and can be applied to separate the observed velocity in shallow sea; at the same time, reasonable boundary conditions of th     

作用于平台桩腿的海冰屈曲破坏冰压力随机过程模型及其参数确定

冰对结构的作用过程是典型的随机过程。本文在冰压力过程为平稳过程的假设下，从渤海海冰对平台桩腿作用的大量实测冰压力数据中，选取了２１条冰屈曲破坏时的冰压力时程曲线，对它们进行谱分析后，得到了单点冰屈曲破坏的压力随机过程的谱密度，并确定了谱参数及其跟环境要素的关系，依据文献［１］中冰压力沿圆柱面的空间分布，建立了绕桩腿的冰压力随机场模型，并得到了作用于桩腿的总冰力随机过程及其谱密度。本文的研究成果为平台结构冰激随机振动和疲劳累积损伤分析提供了荷载基础     

沃克环流上升区赤道低空强西风带的分析

采用ＴＯＧＡ／ＣＯＡＲＥ国际合作考察期间（１９９２年１０月─１９９３年２月）获得的２６７次定点（２°Ｓ，１５５°Ｅ）定时高空大气探测资料，进行计算分析，发现１９９２年１２月─１９９３年２月低空存在一支西风急流，有两次分别持续半月之久，而１９９２年１１月的两次赤道西风急流，因无赤道高空急流配合，仅持续２ｄ便消失。本文还指出：（１）赤道低空西风急流是各种海－气指数产生ＥＮＳＯ异常的重要信息；（２）赤道高空急流有滞后赤道低空急流２ｄ左右的响应关系；（３）赤道低空急流是纬向水汽输送的狭窄通道。     

Numerical Study on Instantaneous Discharge of Unsorted Particle Cloud in Cross Flow

The mixing characteristics of particles such as dredged sediment of variable size discharged into cross flow are studied by a 3D numerical model, which is developed to model the particle-fluid two-phase flow. The Eulerian method with the modified k-ε parameterization of turbulence for the fluid phase is used to solve fluid phase, while a Lagrangian method for the solid phase (particles), both the processes are coupled through the momentum sources. In the model the wake turbulence induced by particles has been included as additional source term in the k-ε model; and the variable drift velocities of the particles are treated efficiently by the Lagrangian method in which the particles are tracked explicitly and the diffusion process is approximated by a random walk model. The hydrodynamic behavior of dumping a cloud of particles is governed by the total buoyancy of the cloud, the drag force on each particle and the velocity of cross-flow. The computed results show a roughly linear relationship between the displacement of the frontal position and the longitudinal width of the particle cloud. The particle size in the cloud and the velocity of cross flow dominate the flow behavior. The computed results are compared with the results of laboratory experiments and satisfactory agreement is obtained.     

The Mathematical and Parameterized Expressions for Wave Induced Excess Flow of Momentum

Wave induced excess flow of momentum(WIEFM)is the averaged flow of momentum over a wave period due to wave presence,which may also be called 3-D radiation stress.In this paper,the 3-D current equations with WIEFM are derived from the averaged Navier-Stokes equations over a wave period,in which the velocity is separated into the large-scale background velocity,the wave particle velocity and the turbulent fluctuation velocity.A concept of wave fluctuating layer(WFL)is put forward,which is the vertical column from the wave trough to wave ridge.The mathematical expressions of WIEFM in WFL and below WFL are given separately.The parameterized expressions of WIEFM are set up according to the linear wave theory.The integration of WIEFM in the vertical direction equals the traditional radiation stress(namely 2-D radiation stress)given by Longuet-Higgins and Stewart.     

The Levant Slumps and the Phoenician Structures: collapse features along the continental margin of the southeastern Mediterranean Sea

Two distinct series of slumps deform the upper part of the sedimentary sequence along the continental margin of the Levant. One series is found along the base of the continental slope, where it overlies the disrupted eastern edge of the Messinian evaporites. The second series of slumps transects the continental margin from the shelf break to the Levant Basin. It seemed that the two series were triggered by two unrelated, though contemporaneous, processes. The shore-parallel slumps were initiated by basinwards flow of the Messinian salt, that carried along the overlying Plio-Quaternary sediments. Seawater that percolated along the detachment faults dissolved the underlying salt to form distinctly disrupted structures. The slope-normal slumps are located on top of large canyons that cut into the pre-Messinian sedimentary rocks. A layer of salt is found in the canyons, and the Plio-Quaternary sediments were deposited on that layer. The slumps are bounded by large, NW-trending faults where post-Messinian faulted offset was measured. We presume that the flow of the salt in the canyons also drives the slope-normal slumps. Thus thin-skinned halokynetic processes generated the composite post-Tortonian structural patterns of the Levant margin. The Phoenician Structures are a prime example of the collapse of a distal continental margin due to the dissolution of a massive salt layer.     

Research on Characteristics of Density Current Under the Action of Waves

In this paper,the characteristics of density current under the action of waves are describedwith the help of flume experiment and theoretical analysis.The study shows that turbid water under the ac-tion of the waves can present three types of motion,i.e.significant stratification,fragile stratification andstrong mixing.The motion of turbid water presents significant stratification when(H/D)/△ρ/ρ~(1/2)≤4.5,generally this state is known as density current.The formulas of motionvelocity,thickness,and discharge of density current moving on horizontal bottom are derived by use of ba-sic equations such as momemtum equation,equation of energy conservation and continuity equation offluid.The time-average velocity and the thickness of density current under the action of waves have a rela-tionship with such parameters as relative density(△ρ/ρ),wave height(H),and water depth(D).Whenthese parameters are determined,the time-average thickness and motion velocity of density current are al-so determined.The relat