波流联合作用下掺氢输气海管泄漏扩散研究

由于氢气密度低,泄漏后的掺氢天然气在海洋波流作用下的浮升扩散行为有别于纯天然气,其泄漏扩散规律有待明晰。利用计算流体力学方法数值模拟了波流联合作用下掺氢输气海管泄漏扩散的过程,结果表明气体泄漏的过程可分为泄漏初期、向上浮升以及横向迁移3个阶段。当掺氢比φH2 <50%时,氢气泡的运动轨迹受天然气泡的影响显著。气体浮升高度和扩散直径的变化与时间成正相关,且随着掺氢比的增大,泄漏气体到达液面所需的时间延长。天然气泡和氢气泡直径在上升的过程中都逐渐增大,两者的浮升速度随浮升高度的增加先增大后减小,天然气泡的浮升速度衰减更快。氢气泡直径随掺氢量的增加而增大,天然气泡直径随掺氢量的增大而减小,两者的上升速度随着掺氢量的增大都表现出先增大后减小的趋势,且氢气泡的上升速度大于天然气泡的上升速度。波长和海流流速越大,泄漏气体的扩散直径越大。     

幔源CO_2充注驱油机理——以松辽盆地南部长岭凹陷为例

分析天然气(CH4)驱油、原油驱水原理,建立了幔源CO2流体的充注驱油模型,在模型中CO2能否形成足够驱动力是驱油的关键。松辽盆地南部长岭凹陷泉四段储层与该模型相符合,在储层中发现幔源CO2与油气混层现象,并且CO2充注时间晚于油气注入时间。根据研究区地质条件,对幔源CO2驱油动力和阻力以及影响其大小的参数(CO2与原油的密度、界面张力、孔喉半径和CO2柱高度)进行分析,得出幔源CO2与原油所产生的浮力足可以突破油气运移阻力(毛细管阻力)。从物理和数学的角度证明幔源成因CO2能够对油气运移起到推动作用。     

流体孔隙剪切变形与岩石圈层状运动机理

在岩石圈层状滑移系统中,如大陆地壳俯冲叠置构造、逆冲推覆断层、拆离断层、滑坡等系统,近于水平的断裂带上存在许多有一定压力的流体孔隙。如果这些流体孔隙是密闭的,便可看成液压千斤顶的液缸。那么,断层面以上滑移体的重量S_T不仅被断裂带上的固体支撑,也受孔隙液体举托力F_f支撑。设F_f与S_T之比为R,那么R就是一个描述滑移体滑移运动难易程度的重要参数。R值越大,滑移体越容易被移动（推动）。而岩石圈中层状滑移系统滑脱面的剪切变形可以显著地改变该参数的大小,从而对层状滑移运动产生深刻的影响。设剪切面上孔隙流体平均压强初始值为P₀,剪切角为α,剪切变形后的平均压强为P（α）,R（α）-α间存在下式的函数关系:R（α）=（R₀P（α））/（P₀cosα）由此可见,R（α）与P（α）正相关,与剪切角α也正相关。P（α）与剪切角α的函数关系目前还难以给出,但剪切过程中,流体孔隙压强保持不变的恒压过程是常见的,即P（α）=P₀。将该等式代入上式（1）,即可获得下式:R（α）=（R₀）/（cosα）,[P（α）=P₀,即恒压过程]由此作出恒压剪切过程的R（α）-α图解。从而发现,恒压剪切过程中,只要滑脱面上密闭流体孔隙存在初始压强,即R₀>0,那么,R（α）随α增长,当剪切角α达到一定程度,R（α）迅速超过1;R₀越大,实现R（α）>1所需的α值越小;反之,所需的α值越大。由此证明,密闭流体孔隙经过恒压剪切变形,导致R（α）>1的基本规律。非恒压过程总体规律不变。R（α）>1意味着层状滑移体完全被孔隙流体的浮托力支持。从而揭示岩石圈俯冲——推覆体及拆离机理和滑坡爆发的可能机制。     

Ascent of migmatites of a high-temperature metamorphic complex due to buoyancy beneath a volcanic arc: a mid-Cretaceous example from the eastern margin of Eurasia

ABSTRACT

Mixtures of melt and residue in a high-T metamorphic complex have a lower density and viscosity than the surrounding host crust, and the mixtures should ascend due to buoyancy. The mixtures are recognized as migmatites in the high-T metamorphic complex. To confirm ascent of migmatites, we conducted numerical simulations of ascent of a model migmatite (buoyant viscous fluid). The numerical simulations show that the model migmatite could rise to shallow levels of a model crust so long as it is continuously produced at the bottom of the model crust. Otherwise it ceases to rise at depth due to loss of buoyancy by cooling. The numerical simulations also show that the model migmatite experiences vertical thinning during the ascent. The ascent mechanism proposed in this paper requires the continuous production of partially melted rocks at the base of the crust, which is provided by a continuous input of energy into the crust from the mantle. Given that high-T metamorphic complexes are associated with igneous activity beneath a volcanic arc, the igneous activity reflects the energy input into the lower crust from the mantle. A high-grade part (migmatites) of a high-T metamorphic complex in the Omuta district of northern Kyushu, southwest Japan, experienced thinning during ascent. Large amount of igneous rocks, such as plutonic and volcanic rocks, are also distributed in northern Kyushu. Zircon U–Pb ages of igneous rocks from northern Kyushu reveal that igneous activity continued from 115 to 93 Ma, and that peak igneous activity at 110–100 Ma was synchronous with the ascent of migmatites of the high-T metamorphic complex in northern Kyushu. Therefore, the numerical simulations may provide an appropriate model of the ascent of migmatites of the high-T metamorphic complex beneath a volcanic arc, at the eastern margin of Eurasia in the mid-Cretaceous.     

Richardson number profiles in laboratory experiments applied to shallow seas

Abstract

A unified analysis is given of the critical conditions for the onset of stratification due to either a vertical or a horizontal buoyancy flux, with tidal or wind stirring.

The critical conditions for the onset of stratification with a horizontal buoyancy flux are found to be of the form of ratios of the tidal slope, or wind setup, to the equivalent surface slope due to the lateral density gradient. These ratios, which are easily determined from sea data, indicate that the profiles of critical flux Richardson Number, averaged over the stirring cycle, are similar to those inferred from the laboratory experiments of Hopfinger and Linden (1982) in which there is zero mean shear turbulence with a stabilising buoyancy flux, and also that the efficiency for the conversion of kinetic energy to potential energy for tidal stirring is similar to that for wind stirring.

The observed much greater efficiency for wind stirring, compared with tidal stirring with a vertical buoyancy flux, is also consistent with the existence of flux Richardson Number profiles in the sea similar to those occurring in the corresponding laboratory experiments. Using the solution of the turbulent kinetic energy equation for the water column, the relative importance of the production of turbulent kinetic energy, and its diffusion by turbulence are assessed, and the critical conditions for the onset of stratification with a vertical buoyancy flux are shown to reduce the classical Simpson—Hunter form.     

Patterns of deposition from experimental turbidity currents with reversing buoyancy

Turbidity currents are turbulent, sediment‐laden gravity currents which can be generated in relatively shallow shelf settings and travel downslope before spreading out across deep‐water abyssal plains. Because of the natural stratification of the oceans and/or fresh water river inputs to the source area, the interstitial fluid within which the particles are suspended will often be less dense than the deep‐water ambient fluid. Consequently, a turbidity current may initially be denser than the ambient sea water and propagate as a ground‐hugging flow, but later reverse in buoyancy as its bulk density decreases through sedimentation to become lower than that of the ambient sea water. When this occurs, all or part of the turbidity current lofts to form a buoyant sediment‐laden cloud from which further deposition occurs. Deposition from such lofting turbidity currents, containing a mixture of fine and coarse sediment suspended in light interstitial fluid, is explored through analogue laboratory experiments complemented by theoretical analysis using a ‘box and cloud’ model. Particular attention is paid to the overall deposit geometry and to the distributions of fine and coarse material within the deposit. A range of beds can be deposited by bimodal lofting turbidity currents. Lofting may encourage the formation of tabular beds with a rapid pinch‐out rather than the gradually tapering beds more typical of waning turbidity currents. Lofting may also decouple the fates of the finer and coarser sediment: depending on the initial flow composition, the coarse fraction can be deposited prior to or during buoyancy reversal, while the fine fraction can be swept upwards and away by the lofting cloud. An important feature of the results is the non‐uniqueness of the deposit architecture: different initial current compositions can generate deposits with very similar bed profiles and grading characteristics, highlighting the difficulty of reconstructing the nature of the parent flow from field data. It is proposed that deposit emplacement by lofting turbidity currents is common in the geological record and may explain a range of features observed in deep‐water massive sands, thinly bedded turbidite sequences and linked debrites, depending on the parent flow and its subsequent development. For example, a lofting flow may lead to a well sorted, largely ungraded or weakly graded bed if the fines are transported away by the cloud. However, a poorly sorted, largely ungraded region may form if, during buoyancy reversal, high local concentrations and associated hindered settling effects develop at the base of the cloud.     

A Case Study on the Effects of Heterogeneous Soil Moisture on Mesoscale Boundary-Layer Structure in the Southern Great Plains, U.S.A. Part I: Simple Prognostic Model

The atmospheric boundary-layer (ABL) depth was observed by airborne lidar and balloon soundings during the Southern Great Plains 1997 field study (SGP97). This paper is Part I of a two-part case study examining the relationship of surface heterogeneity to observed ABL structure. Part I focuses on observations. During two days (12–13 July 1997) following rain, midday convective ABL depth varied by as much as 1.5 km across 400 km, even with moderate winds. Variability in ABL depth was driven primarily by the spatial variation in surface buoyancy flux as measured from short towers and aircraft within the SGP97 domain. Strong correlation was found between time-integrated buoyancy flux and airborne remotely sensed surface soil moisture for the two case-study days, but only a weak correlation was found between surface energy fluxes and vegetation greenness as measured by satellite. A simple prognostic one-dimensional ABL model was applied to test to what extent the soil moisture spatial heterogeneity explained the variation in north–south ABL depth across the SGP97 domain. The model was able to better predict mean ABL depth and variations on horizontal scales of approximately 100 km using observed soil moisture instead of constant soil moisture. Subsidence, advection, convergence/divergence and spatial variability of temperature inversion strength also contributed to ABL depth variations. In Part II, assimilation of high-resolution soil moisture into a three-dimensional mesoscale model (MM5) is discussed and shown to improve predictions of ABL structure. These results have implications for ABL models and the influence of soil moisture on mesoscale meteorology     

Study on Physical Mechanism of Influence on Atmospheric Boundary Layer Depth in the Arid Regions of Northwest China

Using the sounding data of wind, temperature, and humidity in the boundary layer and micrometeorological data on the earth's surface observed in the same period in Dunhuang arid region of Northwest China,this paper researches characteristics of potential temperature, wind, and humidity profiles, confirms the structure and depth of thermodynamic boundary layer in Dunhuang region, and analyzses the relationship of depth of thermodynamic boundary layer with surface radiation, buoyancy flux as well as wind speed and wind direction shear in the boundary layer. The results show that the maximum depth of diurnal convective boundary layer is basically above 2000 m during the observational period, many times even in excess of 3000 m and sometimes up to 4000 m; the depth of nocturnal stable boundary layer basically maintains within a range of 1000-1500 m. As a whole, the depth of atmospheric boundary layer is obviously bigger than those results observed in other regions before. By analyzing, a preliminary judgement is that the depth of atmospheric thermodynamic boundary layer in Dunhuang region may relate to local especial radiation characteristics, surface properties (soil moisture content and heat capacity) as well as wind velocity shear of boundary layer, and these properties have formed strong buoyancy flux and dynamic forcing in a local region which are fundamental causes for producing a super deep atmospheric boundary layer.