A generalized transformation approach for simulating steady-state variably-saturated subsurface flow

A numerical method is proposed to accurately and efficiently compute a direct steady-state solution of the nonlinear Richards equation. In the proposed method, the Kirchhoff integral transformation and a complementary transformation are applied to the governing equation in order to separate the nonlinear hyperbolic characteristic from the linear parabolic part. The separation allows the transformed governing equation to be applied to partially- to fully-saturated systems with arbitrary constitutive relations between primary (pressure head) and secondary variables (relative permeability). The transformed governing equation is then discretized with control volume finite difference/finite element approximations, followed by inverse transformation. The approach is compared to analytical and other numerical approaches for variably-saturated flow in 1-D and 3-D domains. The results clearly demonstrate that the approach is not only more computationally efficient but also more accurate than traditional numerical solutions. The approach is also applied to an example flow problem involving a regional-scale variably-saturated heterogeneous system, where the vadose zone is up to 1 km thick. The performance, stability, and effectiveness of the transform approach is exemplified for this complex heterogeneous example, which is typical of many problems encountered in the field. It is shown that computational performance can be enhanced by several orders of magnitude with the described integral transformation approach.     

承压含水层井流-盖层弯曲效应的解析理论

指出承压含水层盖层的弯曲变形与开采井周围的径向地下水运动存在相互作用, 而这一效应在传统的井流理论中没有被认识到.通过引入弹性薄板理论, 建立了无越流的承压含水层井流-顶板弯曲效应的解析模型, 同时考虑了含水层和水的压缩性, 结果表明Theis井流方程给出的抽水降深偏小.在此基础上推导了有越流承压含水层井流-盖层弯曲效应的偏微分方程, 求出了解析解, 并与传统理论的结果进行了对比, 表明Hantush-Jacob公式计算的降深也是偏小的.在抽水井附近和抽水初期, 传统理论可能导致显著的相对误差.      

岩石单轴压缩下能量与损伤演化规律研究

能量变化是材料各参数变化的本质特征。通过颗粒流模拟岩石单轴压缩试验,从细观力学角度分析能量转化、裂纹扩展及损伤演化规律。颗粒流利用细观颗粒的运动克服了宏观力学理论不易实现试件多裂纹破坏形态的缺点,以此研究能量变化更加合理。应力-应变各阶段能量与微裂纹及损伤存在着相互对应的发展关系。通过3种工况论证了宏-细观力学参数对应关系,采用最小二乘法拟合得出微裂纹与轴向应变呈幂次函数关系。采用割线模量定义损伤变量,取模量加速下降处（对应裂纹加速扩展）为损伤门槛,对应损伤阈值为0.158。3种工况下微裂纹数量与损伤呈线性发展关系,为损伤演化的深入研究提供参考。     

A pleistocene mummified ostracod from the Wairarapa district

A mummified specimen of Scottia insularis Chapman is recorded from beds of Upper Nukumaruan age (Pleistocene) at Gladstone in the Wairarapa Valley. This is the first record of a mummified species from New Zealand.     

空心块体水沙动力及泥沙淤积特性研究

空心块体具有良好的阻水和促淤功能,近年来被广泛用于生态修复工程。本文结合水槽试验及Flow-3D数值模拟,分析了开敞型和半封闭型空心块体的阻水效应和泥沙淤积特性。结果表明：空心块体的开孔率对内部水流流速、紊动强度起主导作用,开孔率较小的半开敞型空心块体减速、制紊效果更强;开敞型和半封闭型空心块体近底层悬沙浓度分别增大56%和75%,两者均有利于促进泥沙在块体内部淤积,近底层水流紊动越强,泥沙淤积形态差异越大;空心块体所营造的低流速、泥沙微淤、内外连通的水沙环境是大型底栖生物的生境需求,半封闭型空心块体内部的低紊动水流结构更有利于大型底栖生物的栖息、繁衍。     

Distribution and expression of gas seeps in a gas hydrate province of the northeastern Sakhalin continental slope, Sea of Okhotsk

Multidisciplinary surveys were conducted to investigate gas seepage and gas hydrate accumulation on the northeastern Sakhalin continental slope (NESS), Sea of Okhotsk, during joint Korean–Russian–Japanese expeditions conducted from 2003 to 2007 (CHAOS and SSGH projects). One hundred sixty-one gas seeps were detected in a 2000 km² area of the NESS (between 53°45′N and 54°45′N). Active gas seeps in a gas hydrate province on the NESS were evident from features in the water column, on the seafloor, and in the subsurface: well-defined hydroacoustic anomalies (gas flares), side-scan sonar structures with high backscatter intensity (seepage structures), bathymetric structures (pockmarks and mounds), gas- and gas-hydrate-related seismic features (bottom-simulating reflectors, gas chimneys, high-amplitude reflectors, and acoustic blanking), high methane concentrations in seawater, and gas hydrates in sediment near the seafloor. These expressions were generally spatially related; a gas flare would be associated with a seepage structure (mound), below which a gas chimney was present. The spatial distribution of gas seeps on the NESS is controlled by four types of geological structures: faults, the shelf break, seafloor canyons, and submarine slides. Gas chimneys that produced enhanced reflection on high-resolution seismic profiles are interpreted as active pathways for upward gas migration to the seafloor. The chimneys and gas flares are good indicators of active seepage.     

Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.     

双平板式透空堤消能效果评价方法

双平板式透空堤具有较为优越的消浪性能备受专家学者关注.目前关于其消浪性能的评价多采用透射波高法开展,仅考虑波高一个参数.本文采用透射波高法、波浪能量法和波能流法分别对平板式透空堤的消能效果进行评价,结果表明,综合考虑波高、水深和周期三个参数的波能流法更加全面与深入.探讨了双平板式透空堤迎浪向与背浪向处波能流的主要影响因...     

规则波作用下台阶地形上墩柱绕流三维数值研究

利用Flow-3D建立三维数值波浪水槽,模拟波浪在不对称台阶地形上的传播。系统研究规则波作用下墩柱周围水流的流动特性,分析墩柱周围的瞬时速度场、涡量场以及KC值变化,不同相位时墩柱前、后水平流速分布情况。结果表明:波浪在台阶地形传播的过程中,墩柱迎水面的涡动结构不够明显;高涡量呈对称状聚集在墩柱的背水面,并形成一对旋转方向相反的涡结构;周期对KC值的影响比波高的影响要明显;墩柱迎水面水平方向流速变化较大,侧面水平流速变化最为剧烈,背面由于受到墩柱的掩护作用水平方向流速变化不大,在墩柱的正面和侧面竖向环流明显。     

Responses of a Mediterranean soft bottom community to short-term (1993-1996) hydrological changes in the Rhone river

The polychaete fauna of muddy bottoms off the Rhone delta (NW Mediterranean) was seasonally sampled at two sites at 70 m depth, from 1993 to 1996. During this period, five severe flood events occurred. A clustering analysis (distance coefficient of Whittaker) and three way fixed factor ANOVAs (site x season x year) showed strong year-to-year changes in species density and community structure, changes that masked seasonal variations. Total density increased by a factor of 3 and density of most species significantly increased following the first flooding event. Changes in the community structure were due to the successive dominance of a few species. The opportunistic species, such as Cossura sp., Mediomastus sp. and Polycirrus sp., exhibited peaks in density 1-3 months after the flood. These peaks were followed by a drastic decline. For species with a long life span, such as Sternaspis scutata, a slower but continuous increase in density was observed which was maintained for several months. The density fluctuations of these species exhibited a good correlation with river flow with time lags of 1-2 years. The successional dynamics observed are explained according to the geographical origin of the floods and the biology and feeding ecology of species.