真空负载方式对疏浚淤泥脱水过程中脱水规律的影响

采用自主改装的变压头黏土渗透仪,模拟疏浚淤泥脱水过程中水分在淤泥颗粒多孔介质中的迁移运动过程,通过测定疏浚淤泥脱水过程中渗透系数与过滤常数、渗水量与滤液量、泥饼含水率与孔隙率,探究了负载压力为100 kPa时,不同负载方式（负载时间及负载压力梯度变化）对疏浚淤泥脱水过程中过滤和渗流规律的影响。研究发现：负载下的过滤和渗流同时进行,负载1 h期间,前期0～40 min内过滤占主导地位,后期40～60 min内渗流起关键作用。负载压力梯度由原来1 h内间隔变化1次,变为以20 kPa为梯度递增变化5次,随着压力变化梯度减小,渗透系数、过滤常数增加,淤泥脱水性能得到较明显的改善,并且前期负载压力越小,脱水效果越好。     

基于时空守恒元和解元(CE/SE)方法的孔隙介质多相流动计算

将时空守恒元/解元(CE/SE)方法推广到二维孔隙介质多相流问题的数值计算中,采用人工压缩法耦合速度和压力,同时结合杂交粒子水平集方法捕捉物质界面.提出一套完整的二维欧拉型孔隙介质非稳态多相不可压缩黏性流动计算方案.通过对溃坝和液滴在重力作用下的运动和变形问题的数值模拟,验证了方法的精度和有效性.在此基础上,提出了一个新的孔隙介质两相流物理模型——双层流体顶盖驱动方腔流.     

任意空间取向TI介质中长排列非双曲动校正公式优化(英文)

常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI(ATI)条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的"精确走时"结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。     

Diurnal evolution of the barrier layer and its local feedback in the central Taiwan Strait

Diurnal evolution of the barrier layer (BL) and its local feedback features in the central Taiwan Strait (119.2oE,24.3oN) during summertime monsoon are investigated using in situ moored observations conducted by the Yan-Ping 2 research vessel in late June 2005.During the initiation phase,for the non-solar radiation tends to be trapped in the upper mixed layer,whereas the solar radiation can penetrate deeply through the mixed layer approaching the thermocline,most heat is accumulated inside the BL inducing a...     

非饱和土一维固结混合流体方法的解析分析

对采用混合可压缩流体方法分析非饱和土一维固结问题的固结方程进行了求解,在得到的解析解的基础上,对影响非饱和土一维固结的因素进行了分析。分析结果表明,在采用混合流体方法计算非饱和土一维固结的孔隙水压力时,所用公式与计算饱和土一维固结的太沙基理论公式基本相同,不同之处在于引入Bishop有效应力系数来体现孔隙气对孔隙水的影响。而在非饱和土孔隙气压的计算公式中除了体现孔隙水对孔隙气的影响参数以外,还有体现孔隙气体的可压缩性对固结影响的参数。在所有影响因素中,影响非饱和土一维固结最重要的因素是孔隙流体的渗流路径。     

城市温室气体清单研究

介绍了国际上城市温室气体清单研究进展。分析了城市清单主流方法体系、模式以及编制原则、边界、范围,并且比较了城市清单和国家清单在方法体系及模式上的差异和其自身特点。重点分析了城市清单编制的“混合模式”和3个尺度范围。最后提出国内城市清单研究面临的困难和建议。     

An analytical solution for one-dimensional contaminant diffusion through multi-layered system and its applications

An analytical solution for one-dimensional contaminant diffusion through multi-layered media is derived regarding the change of the concentration of contaminants at the top boundary with time. The model accounts for the arbitrary initial conditions and the conditions of zero concentration and zero mass flux on the bottom boundary. The average degree of diffusion of the layered system is introduced on the basis of the solution. The results obtained by the presented analytical solutions agree well with those obtained by the numerical methods presented in the literature papers. The application of the analytical solution to the problem of landfill liner design is illustrated by considering a composite liner consisting of geomembrane and compacted clay liner. The results show that the 100-year mass flux of benzene at the bottom of the composite liner is 45 times higher than that of acetone for the same composite liner. The half-life of the contaminant has a great influence on the solute flux of benzene diffused into the underlying aquifer. Results also indicates that an additional 2.9–5.0 m of the conventional (untreated) compacted clay liner under the geomembrane is required to achieve the same level of protection as provided by 0.60 m of the Hexadecyltrimethylammonium (HDTMA)-treated compacted clay liners in conjunction with the geomembrane. Applications of the solution are also presented in the context of a contaminated two-layered media to demonstrate that different boundary and initial conditions can greatly affect the decontamination rate of the problem. The method is relatively simple to apply and can be used for performing equivalency analysis of landfill liners, preliminary design of groundwater remediation system, evaluating experimental results, and verifying more complex numerical models.     

Towards a thermodynamics of immiscible two-phase steady-state flow in porous media

We propose that steady-state two-phase flow in porous media may be described through a formalism closely resembling equilibrium thermodynamics. This leads to a Monte Carlo method that will be highly efficient in studying two-phase flow under steady-state conditions numerically. This work was partially supported by the Norwegian Research Council through grants nos. 154535/432 and 180296/S30.     

Steady flows driven by sources of random strength in heterogeneous aquifers with application to partially penetrating wells

Average steady source flow in heterogeneous porous formations is modelled by regarding the hydraulic conductivity K(x) as a stationary random space function (RSF). As a consequence, the flow variables become RSFs as well, and we are interested into calculating their moments. This problem has been intensively studied in the case of a Neumann type boundary condition at the source. However, there are many applications (such as well-type flows) for which the required boundary condition is that of Dirichlet. In order to fulfill such a requirement the strength of the source must be proportional to K(x), and therefore the source itself results a RSF. To solve flows driven by sources whose strength is spatially variable, we have used a perturbation procedure similar to that developed by Indelman and Abramovich (Water Resour Res 30:3385–3393, 1994) to analyze flows generated by sources of deterministic strength. Due to the linearity of the mathematical problem, we have focused on the explicit derivation of the mean head distribution G _d (x) generated by a unit pulse. Such a distribution represents the fundamental solution to the average flow equations, and it is termed as mean Green function. The function G _d (x) is derived here at the second order of approximation in the variance σ² of the fluctuation (where K _A is the mean value of K(x)), for arbitrary correlation function ρ(x), and any dimensionality d of the flow domain. We represent G _d (x) as product between the homogeneous Green function G _d ⁽⁰⁾(x) valid in a domain with constant K _A , and a distortion term Ψ _d (x) = 1 + σ²ψ _d (x) which modifies G _d ⁽⁰⁾(x) to account for the medium heterogeneity. In the case of isotropic formations ψ _d (x) is expressed via one quadrature. This quadrature can be analytically calculated after adopting specific (e.g.. exponential and Gaussian) shape for ρ(x). These general results are subsequently used to investigate flow toward a partially-penetrating well in a semi-infinite domain. Indeed, we construct a σ²-order approximation to the mean as well as variance of the head by replacing the well with a singular segment. It is shown how the well-length combined with the medium heterogeneity affects the head distribution. We have introduced the concept of equivalent conductivity K ^eq(r,z). The main result is the relationship where the characteristic function ψ^(w)(r,z) adjusts the homogeneous conductivity K _A to account for the impact of the heterogeneity. In this way, a procedure can be developed to identify the aquifer hydraulic properties by means of field-scale head measurements. Finally, in the case of a fully penetrating well we have expressed the equivalent conductivity in analytical form, and we have shown that (being the effective conductivity for mean uniform flow), in agreement with the numerical simulations of Firmani et al. (Water Resour Res 42:W03422, 2006).     

天然气水合物似海底反射层(BSR)AVA特征:双相介质模型

The bottom simulating reflector （BSR） in gas hydrate-bearing sediments is a physical interface which is composed of solid, gas, and liquid and is influenced by temperature and pressure. Deep sea floor sediment is a porous, unconsolidated, fluid saturated media. Therefore, the reflection and transmission coefficients computed by the Zoeppritz equation based on elastic media do not match reality. In this paper, a two-phase media model is applied to study the reflection and transmission at the bottom simulating reflector in order to find an accurate wave propagation energy distribution and the relationship between reflection and transmission and fluid saturation on the BSR. The numerical experiments show that the type I compressional （fast） and shear waves are not sensitive to frequency variation and the velocities change slowly over the whole frequency range. However, type II compressional （slow） waves are more sensitive to frequency variation and the velocities change over a large range. We find that reflection and transmission coefficients change with the amount of hydrate and free gas. Frequency, pore fluid saturation, and incident angle have different impacts on the reflection and transmission coefficients. We can use these characteristics to estimate gas hydrate saturation or detect lithological variations in the gas hydrate-bearing sediments.