油气微渗漏的垂向结构

<正> 油气微渗漏在垂向上具有一定的分区结构,这种结构是由油气藏和围岩的孔隙结掏、空间结构、地质结构和物理化学环境所决定的。油气向上渗漏的微泡迁移机制被人们普遍地肯定和接受。这种理论认为,烃类气体呈胶粒气泡型在浮力作用下以布朗运动方式沿着油气藏上覆岩层中充满水的、相互连通的“微裂隙”网络向地表迁移;而油气藏中的微量元素(油气灰分)也可以在气水界面张力的作用下被吸附在微泡表面,与微泡一同渗漏迁移。在这种机制中,微渗漏过程主     

波生沿岸流数值模拟

为了更好地研究近岸海域波生沿岸流,建立了基于高阶Boussinesq水波方程的波生沿岸流时域数值模型。控制方程在中等水深范围内具有较好的色散性和变浅作用性能,同时具有二阶完全非线性特征,适合描述近岸区域波浪强非线性运动。通过采用松弛造波方法实现了非线性波浪的无反射入射,采用周期性侧边界条件模拟开敞边界。通过数值试验,讨论了模型中主要参数对数值结果的影响。利用率定后的参数模拟了均匀坡度海岸上产生的沿岸流,通过和实验数据的对比验证了模型的准确性和适用性。利用模型数值模拟了不同波浪入射条件(包括周期、波高和波浪入射角度)对波生沿岸流的影响。     

近岸波生流运动三维数值模拟及验证

开发建立了近岸波生流运动三维数值计算模式。模式中,引入了三维时均剩余动量、破波表面水滚、波浪水平与垂向紊动作为主要驱动力,同时考虑了波流共同作用的底部剪切力。推导了可综合反映底坡、能量传递率和密度影响的水滚能量传输方程;将Larson-Kraus的二维波浪水平紊动系数表达式拓展至三维。采用大量实测数据和文献资料测试验证了所建模式,表明所建模式可有效模拟波浪增减水、底部离岸流、沿岸流、裂流、堤后环流等不同维度的波生流现象。此外,研究也表明破波水滚效应可解释波生流峰值向岸推移的物理现象,从而在模拟中不能忽略;破波带内沿岸流速垂向较为均匀的现象与波浪附加垂向紊动有关。     

垂向二维非恒定流及悬浮物分布模型研究

采用σ坐标变换拟合复杂边界、VOF方法追踪自由水面变化和SIMPLE算法求解控制方程等建立垂向二维非恒定水流及悬浮物分布数学模型。该模型适用于变量在横向上分布基本均匀的水流流动特性及悬浮物扩散迁移过程的模拟,与刚盖假定或静水压强假定等条件下建立的垂向二维数学模型相比更符合实际情况。分别对概化水库水流结构的沿程变化、明渠均匀流及带有槽沟的明渠非均匀流及悬浮物浓度分布进行了模拟,计算结果与实测值符合良好。模型的应用可为河道或水库的洪水调度、水污染控制、取水工程优化等提供参考。     

磨刀门河口洪季波生流及其泄洪影响

磨刀门已由"径流型"向"径流-波浪型"河口转变,波浪已是该河口主要动力之一,但波浪对河口洪季水流及泄洪的影响缺少研究。在2-D潮流数学模型中添加随潮位实时变化的波浪辐射应力,建立波浪潮流耦合数学模型;波浪求解采用缓坡方程,背景水深由潮流模型实时提供,可通过比较考虑和未考虑波浪影响的河口流场来分析波浪对泄洪的影响。在年均常浪作用下,磨刀门河口洪季涨落潮阶段均有明显的波生环流结构。由于波浪作用方向向陆,波生流减弱了浅滩区的向海余流,增大了浅滩向陆余流;受浅滩向海余流减弱影响,河口动力自调整后形成归槽水流,促使深槽内向海余流增大。波浪有顶托河口泄洪之势,可改变滩槽泄洪分配比例;年均常浪的波高较小,其对潮流及泄洪的影响区域限制在浅水区,故对泄洪的负面影响有限。     

波浪载荷导致黄河口潮坪沉积物垂向运移现场观测研究

2005年8月7~8日,现代黄河三角洲刁口地区经历了一次由台风麦莎引起的风暴潮过程。通过对比分析风暴潮前后受保护潮坪滩面沉积物的粒度特征变化情况,发现了风暴潮期间在强烈的波浪载荷导致的渗流作用下,粉质土海床出现了极细粉粒由下向上运移并输出滩面的现象。结合现场试验期间采集到的孔隙水压力监测数据,本文基于海床动力响应的观点分析了其成因机制。     

规则波作用下刚性挺水植物波生流试验研究

通过在波浪水槽中进行一系列物理模型试验研究规则波作用下刚性挺水植物波生流特性,测试了2种水深和不同入射波况条件下的波生时均流速并对比了挺水植物有茎与无茎时波生流的情况,基于假设检验分别对挺水植物的拖曳力长度尺度、植物特征参数建立模型。研究结果表明:在茎密度为60株/m2时,波生流仍以根部作用为主要因素,对沉水植物公式进行拓展,据此提出挺水植物时均流速最大值预测公式;茎部对时均流速均方根具有衰减作用,且茎部在根部及根部上方区域产生的相对衰减强度,与厄塞尔数呈抛物线函数关系;在簇状根系分布下,茎部对垂向流速结构起稳定作用,尤其是在时均流速最大值所对应的量纲一垂直位置。     

邯郸娄里金矿床原生晕垂向分带特征研究

娄里金矿位于邯郸洪山杂岩体北侧的火山凝灰岩中。该类型金矿在太行山南段为首次发现,其成矿规律和深部成矿潜力的研究对圈定成矿远景区及增加矿山资源量有重要的指导意义。据此系统采集了矿区ZK11-2和ZK16-1钻孔岩心样品,在分析区域成矿地质规律和矿床地质特征的基础上,采用格里戈良分带指数法和广义衬值法两种方法计算了元素垂向分带序列,结果显示不同方法虽在元素排列顺序上存在一定差异,但在反映热夜活动期次及元素空间叠加模式上基本相同。两个钻孔中均出现了明显的"前尾晕共存"及"反向分带"现象。在此基础上结合微量元素垂向变化规律,分析了区内金矿体形成的热夜活动过程。研究表明两个钻孔的垂向分带都由两个完整的小循环组成,说明区内金矿体的形成至少经历了两次热夜活动叠加。     

天津地区地热流体化学成分垂向演化特征研究

天津地区地热流体化学成分具有明显的垂向分层演化特征,由浅至深不仅有TDS增高、水化学类型由HCO3-Na型向SO4-Na或Cl-Na型演化的特征,而且具有分层演化特征,即不同热储层段的地热流体水化学特征及演化规律不同,说明其补给、径流、排泄自成系统。     

The vertical turbulence structure of experimental turbidity currents encountering basal obstructions: implications for vertical suspended sediment distribution in non‐equilibrium currents

Large roughness features, caused by erosion of the sea floor, are commonly observed on the modern sea floor and beneath turbidite sandstone beds in outcrop. This paper aims to investigate the effect of such roughness elements on the turbulent velocity field and its consequences for the sediment carrying capacity of the flows. Experimental turbidity currents were run through a rectangular channel, with a single roughness element fixed to the bottom in some runs. The effect of this roughness element on the turbulent velocity field was determined by measuring vertical profiles of the vertical velocity component in the region downstream of the basal obstruction with the Ultrasonic Doppler Velocity Profiling technique. The experiments were set up to answer two research questions. (i) How does a single roughness element alter the distribution of vertical turbulence intensity? (ii) How does the altered profile evolve in the downstream direction? The results for runs over a plane substrate are similar to data presented previously and show a lower turbulence maximum near the channel floor, a turbulence minimum associated with the velocity maximum, and a turbulence maximum associated with the upper flow interface. In the runs in which the flows were perturbed by the single roughness element, the intensity of the lower turbulence maximum was increased between 41% to 81%. This excess turbulence dissipated upwards in the flow while it travelled further downstream, but was still observable at the most distal measurement location (at a distance ca 39 times the roughness height downstream of the element). All results point towards a similarity between the near bed turbulence structure of turbidity currents and free surface shear flows that has been proposed by previous authors, and this proposition is supported further by the apparent success of a shear velocity estimation method that is based on this similarity. Theory of turbulent dispersal of suspended sediment is used to discuss how the observed turbulent effects of a single large roughness element may impact on the suspended sediment distribution in real world turbidity currents. It is concluded that this impact may consist of a non‐equilibrium net‐upwards transport of suspended sediment, counteracting density stratification. Thus, erosive substrate topography created by frontal parts of natural turbidity flows may super‐elevate sediment concentrations in upper regions above equilibrium values in following flow stages, delay depletion of the flow via sedimentation and increase their run‐out distance.     

Sediment suspension under waves and currents: time scales and vertical structure

Field measurements of the vertical structure of near-bed suspended sediment concentrations were obtained from arrays of fast response optical backscatter suspended solids sensors to examine the time-dependent response of sediment resuspension to waves and currents and the constraints imposed by bedforms. Data were recorded from both a nonbarred, marine shoreface and a barred lacustrine shoreface, under both shoaling and breaking waves (significant heights of 0·25–1·50m; peak periods of 3 and 8 s) and in water depths of 0·5–5·0 m. Sediment concentrations are positively correlated with increasing elevation above the bed, but lagged in time. The time lag varies directly with separation distance between measurement locations and inversely with the horizontal component of the near-bed oscillatory velocity. Both the presence of wave groups and the settling velocities of the sediment particules in suspension influence the temporal changes in concentration at a given elevation. Sediment concentrations appear to respond more slowly to the incident wind-wave forcing with distance away from the bed as a result of two factors: (1) the sequential increase in concentration induced by a succession of large waves in a group; and (ii) the relative increase in finer sediments with smaller settling velocities. Bedforms interact with the near-bed horizontal currents to impose a distinct constraint upon the timing of suspension events relative to the phase of the fluid motion, and, therefore, the vertical structure of the suspended sediment concentration at a range of time scales. The near-bed concentrations appear to be strongly dependent upon the vertical convection of sediment associated with the ejection from the wave boundary layer of separation vortices generated in the lee of ripple crests. Concentration gradients in the presence of vortex ripples are large, as are the correlation between concentrations measured at different elevations within the fluid.     

Flow structure of turbidity currents

A two‐dimensional numerical model is used to describe the flow structure of turbidity currents in a vertical plane. To test the accuracy of the model, it is applied to historical flows in Bute Inlet and the Grand Banks flow. The two‐dimensional spatial and temporal distributions of velocity and sediment concentration and non‐dimensionalized vertical profiles of velocity, turbulent kinetic energy and sediment concentration are discussed for several simple computational currents. The flows show a clear interaction between velocity, turbulence and sediment distribution. The results of the numerical tests show that flows with fine‐grained sediment have low vertical and high horizontal gradients of velocity and sediment concentration, show little increase in flow thickness and decelerate slowly. Steadiness and uniformity in these flows are comparable for velocity and concentration. In contrast, flows with coarse‐grained sediment have high vertical and low horizontal velocity gradients and high horizontal concentration gradients. These flows grow considerably in thickness and decelerate rapidly. Steadiness and uniformity in flows with coarse‐grained sediment are different for velocity and concentration. The results show the influence of spatial and temporal flow structure on flow duration and sediment transport.     

Fine vertical structure of the galactic gaseous disk

The three-dimensional evolution of an ensemble of N particles (N = 8 × 10⁵) in the external gravitational field of a galaxy perturbed by a spiral density wave is considered. The particles simulate clouds of interstellar gas, and inelastic two-body collisions between them are taken into account. The three-dimensional structure of the gaseous galactic layer and the vertical profile of the spiral arms are computed. It is shown that: (1) the local thickness of the gaseous galactic disk has a minimum where the volume gas density has a maximum (the maximum density of the interstellar medium is shifted downstream relative to the galactic shock front), (2) the configuration of the vertical profile of the spiral arms changes radically when the corotation region is crossed. Our first result explains the negative correlation between the thickness of the gas layer and the density derived from neutral-hydrogen observations. The second result can be used in the next generation of neutral-hydrogen observations to localize the corotation radius in the Galaxy.     

On the spatial structure of currents across the Chesapeake and Delaware Canal

The Chesapeake and Delaware (C&;D) Canal is a man-made waterway connecting two of the largest estuaries on the east coast of the United States: Chesapeake Bay and Delaware Bay. A set of current meter data collected during April–May 1975 along two cross-sections of the C&;D Canal was used to examine the spatial distributions of the currents at tidal and subtidal time scales. The different responses of the Chesapeake and Delaware Bays to tidal and wind forcing produce significant differences in sea level fluctuations between the two ends of the canal. These alongcanal surface slopes produce significant barotropic current fluctuations at both tidal (semidiurnal and diurnal) and subtidal (2-d to 3-d) time scales. Under the influence of bottom friction, the barotropic currents near the surface are stronger than those at depth, but these currents do not exhibit significant lateral variations across the canal. On the other hand, the long-term flow in the canal exhibits strong lateral variability with eastward flow off the south shore of the canal and westward flow off the north shore of the canal. The lateral structure of the long-term flow may carry significant implications for the long-term exchange of material between the two bays.     

波浪作用下沙坝剖面形成过程的数值模拟

建立了波浪、底部离岸流、泥沙运动和沙滩剖面演变耦合数学模型,模拟了实验室中波浪作用下沙坝剖面的形成过程,讨论了各个物理参数对剖面上的水动力和地形变化的影响。通过参数率定,较好地复演了不规则波作用下沙滩地形由均匀斜坡向沙坝剖面的演变,波高、底部离岸流、输沙率和剖面演变的计算结果与实测数据吻合良好。采用数值实验和误差分析方法,讨论了破碎波水滚倾角、泥沙扩散系数和床面休止角等物理参数对计算结果的影响。研究表明,所建立的模型能较好地描述波浪作用下沙坝剖面形成过程中的重要物理机制。     

Jet structure of electric currents in coronal magnetic loops

We analyze the properties of the electric-current distribution over the cross sections of fairly dense coronal magnetic flux tubes in which the plasma pressure exceeds the magnetic pressure, so that the equilibrium is maintained by the ambient magnetic field. If the plasma is fully ionized, the distributions of the longitudinal and azimuthal currents over the cross section of the loop have the same spatial scale as the pressure distribution. However, even a small number of neutral atoms in the corona (with a mass fraction of the order of 10^?5, taking into account the partial ionization of helium) substantially modifies the current distribution over the tube cross section: in this case, a considerable fraction of the full current flowing along the tube is concentrated in a thin region near the axis with a radius of the order of (10^?2–10^?3)r ₀ (where r ₀ is the characteristic scale of the plasma-pressure distribution over the tube), thus forming a sort of a jet current. This comes about because the pattern of the conductivity anisotropy is substantially modified in the presence of ion-atom collisions in the magnetoactive plasma of the tube, and the Cowling conductivity dominates over the Hall and Pedersen conductivities. The high current density near the axis of the tube can ensure heating of the plasma to coronal temperatures via Joule dissipation.     

A semi-analytical method for the wave-induced seabed response

A general semi-analytical method is presented for the analysis of seabed response under wave action during a storm. The seabed is idealized as a poro-elastic medium filled with a single compressible fluid with anisotropic flow. The coupled process of fluid flow and deformation of soil skeleton is formulated in the framework of Biot's theory. The analysis for the response of homogeneous seabed of finite thickness under a plane progressive wave is developed first, followed by an extension for the case of a layered seabed. A generalization for three-dimensional response of seabed is also developed for a general wave field which provides the analyses for seabed response under short-crested or standing waves in the vicinity of a structure. Some numerical examples illustrating the proposed analyses are also presented.