An analytical solution to multi-component NAPL dissolution equations

This note presents a novel method for determining the changing composition of a multi-component NAPL body dissolving into moving groundwater, and the consequent changes in the aqueous phase solute concentrations in the surrounding pore water. A canonical system of coupled non-linear governing equations is derived which is suitable for representation of both pooled and residual configurations, and this is solved. Whereas previous authors have handled such problems numerically, it is shown that these governing equations succumb to analytical solution. By a suitable substitution, the equations become decoupled, and the problem collapses to a single first-order equation. The final result is expressed implicitly, with time as a function of the number of moles of the least soluble component, m₁. The number of moles of each other component is expressed explicitly in terms of m₁. It is shown that the time-m₁ relationship has a well behaved inverse. An example is given in which the analytic solution is verified against traditional finite difference analysis, and its computational efficiency is shown.     

Analytical models of contaminant transport in coastal aquifers

The Henry formulation, which couples subsurface flow and salt transport via a variable-density flow formulation, can be used to evaluate the extent of sea water intrusion into coastal aquifers. The coupling gives rise to nontrivial flow patterns that are very different from those observed in inland aquifers. We investigate the influence of these flow patterns on the transport of conservative contaminants in a coastal aquifer. The flow is characterized by two dimensionless parameters: the Péclet number, which compares the relative effects of advective and dispersive transport mechanisms, and a coupling parameter, which describes the importance of the salt water boundary on the flow. We focus our attention on two regimes – low and intermediate Péclet number flows. Two transport scenarios are solved analytically by means of a perturbation analysis. The first, a natural attenuation scenario, describes the flushing of a contaminant from a coastal aquifer by clean fresh water, while the second, a contaminant spill scenario, considers an isolated point source.     

Effective macroscopic transport parameters between parallel plates with constant concentration boundaries

A macroscopic transport model is developed, following the Taylor shear dispersion analysis procedure, for a 2D laminar shear flow between parallel plates possessing a constant specified concentration. This idealized geometry models flow with contaminant dissolution at pore-scale in a contaminant source zone and flow in a rock fracture with dissolving walls. We upscale a macroscopic transient transport model with effective transport coefficients of mean velocity, macroscopic dispersion, and first-order mass transfer rate. To validate the macroscopic model the mean concentration, covariance, and wall concentration gradient are compared to the results of numerical simulations of the advection–diffusion equation and the Graetz solution. Results indicate that in the presence of local-scale variations and constant concentration boundaries, the upscaled mean velocity and macrodispersion coefficient differ from those of the Taylor–Aris dispersion, and the mass transfer flux described by the first-order mass transfer model is larger than the diffusive mass flux from the constant wall. In addition, the upscaled first-order mass transfer coefficient in the macroscopic model depends only on the plate gap and diffusion coefficient. Therefore, the upscaled first-order mass transfer coefficient is independent of the mean velocity and travel distance, leading to a constant pore-scale Sherwood number of 12. By contrast, the effective Sherwood number determined by the diffusive mass flux is a function of the Peclet number for small Peclet number, and approaches a constant of 10.3 for large Peclet number.     

p—范分布的近似表示

p-范分布是一个包含拉普拉斯分布、正态分析、均匀分布等常见分布的分布族。用p-范分布描述观测误差的统计特性，只需假定误差的分布为单峰、对称，因此、p-范分布似然平差可以避免事先假定误差的具体分布模式，而在平差过程中确定未知参数及误差的分布具有自适应的特点。但是p-范分布的密度函数比较复杂，不利于理论分析和实际应用。 的研究表明，p-范分布可以近似地表示为拉普拉斯分布与正态分析或正态分布均均匀分布的线性组全。p-范分布与本文给出的近似分布具有相的前四阶矩。由于拉普拉斯分布。正态分布。均匀分布的密度函数都比较简单，用近似分布代替p-范分布会使相关的问题得到简化。     

拟合推估新用——新逼近模式

作为拟合推估的新用之一，提出一种带有随机量的新逼近模式，将拟合推估的通行算法，去短取长，藉以推得良好的插值。     

圆弧形凹陷地形表面覆盖层对入射平面SV波的影响

利用Fourier Bessel级数展开法，给出了表面具有覆盖层的圆弧形凹陷地形对入射平面SV波散射问题的一个解析解，并利用该解分析了不同深宽比凹陷地形表面覆盖层刚度和厚度对入射SV波的影响．结果表明，凹陷地形表面覆盖层的存在，即使厚度很薄，也会显著增加对入射SV波的放大作用，该放大作用可达到单一凹陷地形的2.5倍以上．覆盖层刚度和厚度的变化对入射平面SV波也具有很大影响．      

宁夏南部地区地震活动特征

研究了宁夏南部地区（35°20＇～37°40＇N，104°30＇～106°30＇E）小震综合机制解的变化特征．发现在中强地震前，该地区小震综合机制解的P轴方位偏离平均值，向E或SE偏转，其仰角有增大的趋势，同时小震纵波（Pn和Pg）初动符号一致．6次中强地震的震源机制解表明，该地区P轴方位与主压应力方向基本一致，因此该地区发生破坏性地震的可能性不大．近几年来的地震活动图像显示，该地区的地震活动有向中卫、中宁和同心地区集中的迹象．     

天津及邻区1999～2000年地震震源机制研究

文中使用的是天津市、河北省地震台的 P,Pn波初动方向资料 ,研究了 1 999～ 2 0 0 0年发生在天津及邻区的 4个地震的震源机制。其中 1 999年 8月 1 1日河北省大城县 ML 3.1地震发生在大城断裂上 ,震源机制解的节面 (走向 2 0 8°,倾角 54°,滑动角 1 2°)与大城断裂的走向一致。2 0 0 0年 1月 2 6日蓟县 ML 3.3地震发生在蓟县 -邦均断裂上 ,震源机制解的节面 (走向 2 85°,倾角 77°,滑动角 1 38°)与蓟县 -邦均断裂的走向一致。表明这两个断层在活动 ,这两个地震发生在多年来 (有感 )地震不多的地区 ,并且是发生在首都圈中部 ,应该引起注意     

盆地超压层段非幕式突破期的地热场模型数值解法

盆地超压层段的地热场模拟对超压成因研究及油气生成，排放机理分析有重要意义，但求解超压层段的热传递方程的数值解法至今仍没有被很好地解决，问题的关键在于确定超压流体速度场，由于超压流体发生幕式突破之前的排出速率极低，其速度场v=vxi vzk近似满足条件Эux/Эx Эuz/Эz＝0，可将其近似地视为稳定的不可压缩的无源流体，利用这一条件及相应的边界压力条件，可使整个计算过程得到简化，在此基础上所建立的地热场模型有限元数值解法和模拟软件，能够实现对含油气盆地在幕式突破之前的超压层段地热场进行动态模拟。     

边界条件对土层粘弹性地震反应的影响

边界条件是土层粘弹性地震反应必须考虑的问题。对于地表面，一般作为自由边界，而对于地下深处基岩面，边界条件有两种处理方法：一是静止边界，二是非静止边界的。首先给出了在时间域内单层地基一维土层粘弹性地震反应封闭解析解，然后从理论上分析边界条件对土层粘弹性地震反应的影响，认为采用静止边界进行土层地震反应计算得到的地震动加速度时程符合天然地震记录形式，而采用非静止边界进行土层地震反应计算得到的地震动加速度时程初始阶段出现较大幅度的振荡现象，随着土层厚度增加，这种振荡现象逐渐减弱直至消失。当土层阻尼较小或剪切波速增长时，振荡现象加剧；当土层阻尼比较大或剪切波速减小时，振荡现象减弱。此外，采用非静止边界进行土层地震反应计算，表现出随着土层厚度的增加，土层对基岩输入地震动的放大作用逐步转化为吸收作用。