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).     

约束等效源稳定向下延拓方法研究

在等效源向下延拓方法的理论基础上,研究提出加入观测平面以下空间实测数据作为约束的等效源约束向下延拓法,用于二维位场数据向下延拓.与传统的向下延拓方法和等效源法延拓相比,等效源约束向下延拓法抗噪声性强,能延拓出由于延拓深度过大而衰减至难以探测的高频信号,使延拓结果更接近真实值.理论重力模型数据的延拓结果显示该方法有效,延...     

The Influence of Thrust Fault Structure on Cross-fault Short-leveling Survey

Aiming at different variation patterns of cross-fault short-leveling before earthquakes, the paper establishes the 2-D finite-element models with different crustal stratification and fault occurrence perpendicular to Longmenshan fault zone. By contact analysis and viscoelastic finite element method, the influence of fault structure on cross-fault short-leveling is obtained under the same constraint conditions, the results show that:with the increase of the horizontal projection distance of fault, the cumulative displacements of surface increase gradually in the models with fixed dip angles of the fault plane (model 1). However, when the horizontal projection distance exceeds 20km, the influence of fault''s dip angle on the cumulative displacements of surface short-leveling will decrease significantly, and the cumulative displacements are maintained at about 1.5m. However, in the listric fault models (model 2),when the projection distance is less than 20km, the listric fault structure impedes the sliding of fault. The short-leveling variation rate is only half of model 1; as a result, the ability to reflect the regional stress enhancement by cross-fault short-leveling is further weakened. But when the horizontal projection distance exceeds 25km, the cumulative displacements significantly increase, with the maximum displacement reaching 1.75m. The results of equivalent stress show that the listric fault structure causes a sudden increasement in equivalent stress when the horizontal projection distance is 10km, higher equivalent stress values are accumulated between the projection distance of 5-20km, and then high-low stress difference zones are formed at the bottom of the fault plane and near the transition zone of low-high dip angle.     

利用电磁台站观测数据推算汶川MS8.0地震等效电偶矩

本文利用地下电偶极子产生的电磁场模型, 尝试模拟和推算地震前地面电磁台站观测数据所需震源深度处等效电偶矩。 结果显示, 在观测频率为1 Hz, 地下平均电导率取7.0×10^-4 S/m 时, 2008年汶川M_S8.0地震前电磁异常高潮期, 距离震中1440 km 的高碑店台NS测向观测到1.3 mV/m的电场强度所需震源深度处等效电偶矩量级为10¹³A·m; 若考虑地震主破裂长度150 km, 则产生的地电流量级为10⁸ A。 同时, 这一等效电偶源在地面产生电场的2D分布图显示出很强的方向性(或不均匀性); 信号强度和电流源方向存在着密切的关系, 且电磁场值随观测距离增大而衰减迅速。 这部分地解释了地震相关电磁异常信号的方向性(或不均匀性)的观测事实; 如果观测点位于非强信号区, 即使在有效的探测范围内也可能记录不到异常信息, 这对地震监测台网仪器选址和布置提出了要求, 具有现实的指导意义。     

基于等效线性化方法的一维土层地震反应通用计算程序对比

基于等效线性化的一维土层地震反应计算是目前国内外普遍采用的方法,国外的SHAKE91、DEEPSOIL和我国的LSSRLI-1即是根据这一方法编制的通用计算程序。本文采用这3个程序进行了不同地震波、不同输入地震动幅值下不同场地类型的土层地震反应计算,并对三者的结果进行了全面的比较分析。结果表明:①SHAKE91和DEEPSOIL程序的计算结果完全相同;②当土层最大剪应变均采用时域计算时,LSSRLI-1程序的计算结果与SHAKE91和DEEPSOIL程序基本相同,但有微小差别,其原因是:在基于等效剪应变通过离散形式的剪切模量和阻尼比随等效剪应变变化的关系曲线确定等效剪切模量和阻尼比时,DEEPSOIL和SHAKE91采用的插值方法与LSSRLI-1不同;③当LSSRLI-1程序采用频域经验关系计算土层最大剪应变时,特别是在强地震动输入下得到的土层地表加速度峰值和加速度反应谱与另外两个程序的计算结果有差别,且土层最大剪应变随着输入加速度的增大出现较大的差别。因此,本文建议:当采用LSSRLI-1程序计算土层地震响应时,应使用程序中的时域解方法代替以往默认的频域经验关系方法。     

正则化等效层重力向下延拓方法

近年来,利用时移微重力技术进行储层开发监测受到国内外学者广泛关注.时移微重力观测数据存在信噪比低,信号弱的问题,难以实现储层内物质运移的定量解释.为压制数据噪声,增强有效弱信号,本文研究了利用Tikhonov正则化方法反演等效层（源）,并由等效源实现重力场向下延拓的方法;在此基础上,本文推导了波数域正则化等效源向下延拓算子.针对向下延拓场幅值衰减问题,提出了正则化等效源迭代补偿算法.通过模拟数据实验研究了不同深度正则化等效源滤波算子及向下延拓算子的波数响应;与波数域Tikhonov正则化向下延拓方法相比,正则化等效源向下延拓方法的延拓精度更高、更稳定.最后,将基于迭代补偿的正则化等效源向下延拓技术应用于实测时移微重力数据证实了该方法能够有效增强局部异常,实现时移微重力数据大深度稳定向下延拓.     

多层等效源曲面磁异常转换方法

在磁异常数据处理中,利用等效源技术重构磁异常场具有较好的稳定性和较高的计算精度,因而被广泛应用.传统方法是采用设置在近地表的单层等效源拟合实测磁异常数据,尽管拟合精度很高,但向上延拓之后往往会出现较大的拟合误差,即可能存在磁异常信号的"泄漏",尤其在原始数据中存在背景场时更容易出现此种误差.本文提出一种多层等效源技术方案,应用分布于不同深度范围内的等效源模拟实测数据,减少了等效源参数设置的盲目性.理论模型试验表明,采用多层等效源方法重构的磁异常及其梯度与分量,较单层等效源方法具有更高精度,可以吸收更完整的实测磁异常信息.论文详细地讨论了如何优化多层等效源设置、等效源参数选择以及计算方法,通过二维和三维理论模型试验,验证了在复杂条件下多层等效源方法的可行性和适应性,并且将该方法应用于广西某地的实测磁异常数据转换之中,取得了较好的应用效果.     

水平分层土层系统等效阻尼比的简化计算方法

在实际工程场地中,很多土层可视为水平分层,各层土的物理和力学性质存在差异,其中包括土的振动阻尼比。本文讨论水平分层土层系统的等效阻尼比的近似计算方法,基于5个不同的加权函数推导了10种等效阻尼比的计算公式。通过2个算例,分别以等效阻尼比为参数计算水平分层土层的地震反应,并与准确解相比较,分析了不同等效阻尼比近似计算方法的计算精度。数值结果表明,若等效阻尼比计算方法选择不恰当,会导致土层地震反应的计算结果出现较大误差。针对2种不同类型的水平分层土层,建议采用基于三角形分布的加权函数来计算土层系统的等效阻尼比。     

调频液体阻尼器(TLD)的等效力学模型研究

本文根据流体力学理论,建立了圆柱形ＴＬＤ中晃动液体的势流场,并从液体晃动的动力效应等效原则出发,导出了晃动液体的等效力学模型,并阐明其适用范围。文量计算的结果表明,在工程应用中,仅用ＴＬＤ的一阶晃动等效力学模型就能满足工程需要。     

Atmospheric driving mechanisms of extreme precipitation events in July of 2017 and 2018 in western Japan

The occurrence of extreme precipitation events is now a serious concern in recent years in Japan. This study explores the atmospheric driving mechanisms of two extreme precipitation events occurred during 5–6 July 2017 and 5–8 July 2018 over western Japan. We identified that the atmospheric transport of large amounts of moisture and wind streams with wind speed of minimum 15 m s⁻¹ from south of Japan towards north on the days before these torrential precipitation events are mainly responsible for the July 2017 and July 2018 floods over western Japan. However, the contributions from the moisture advections (both vertical and horizontal) to the atmospheric water budget plays key roles to intensify the precipitations during the said torrential events. We also find that the prominent moisture flux convergence and well-developed moist conditions mainly maintain these heavy precipitation events over the downpour affected area. Our overall analysis suggests that the atmospheric factors driving to all these two heavy precipitation events are qualitatively robust and vital to explain the mechanism of extreme precipitations.