Toward a new conceptual model for groundwater flow in merokarst systems: Insights from multiple geophysical approaches

Merokarst aquifers — relatively thin (<1–2 m) karstified carbonate units interbedded between mudstone, shale, or sandstone — constitute a significant proportion of carbonate terrain and underlie a large portion of the west- and south-central USA, yet few advances have been made in our understanding of porosity development and flow-path generation in these complex systems in decades. Toward this end, we used a multi-geophysical approach at the well-studied Konza Prairie Biological Station (KPBS), a part of the larger Flint Hills (25,734 km²), underlain by thin limestone units (1–2 m thick) interbedded with mudstone/shale units (2–4 m thick), to elucidate hydrologic connectivity and potential controls on known groundwater flow directions. We combined electrical resistivity tomography (ERT), surface and borehole nuclear magnetic resonance (NMR), and ground penetrating radar (GPR) measurements across a low order catchment where over 25 boreholes and groundwater wells sampling perched aquifers could be used to constrain interpretation of lithology, potential flow paths, and permeability. Data revealed that groundwater export may be an unappreciated component of lateral-flow-dominated models used to represent merokarst in that: (a) potentiometric surfaces from two limestone units showed groundwater flows toward a hydrologic depression, opposite the direction of stream flow, in the upstream portion of the catchment, (b) long term measures of groundwater levels revealed a greater variance and overall water storage in this same upstream area compared to wells near the outlet, and (c) ERT and NMR results indicate pronounced lateral heterogeneity within a given unit, suggestive of a greater degree of vertical hydrological connectivity than usually considered for horizontally-layered merokarst. Our data suggest vertical connectivity can shunt water to depth in these “sandwiched” merokarst aquifers, yielding atypical groundwater flow directions and unrealized deep export of weathering solutes and carbon.     

地球物理联合探测在识别岩溶地面塌陷精细结构中的应用—以武汉市为例

岩溶地面塌陷是隐伏岩溶区典型的地质灾害类型之一,查明其精细地质结构对城市地下空间开发利用和地质灾害防治都具有重要意义。本文以武汉市为例,梳理了以覆盖层厚度与结构、岩溶发育程度和地下水为核心的岩溶塌陷探测地质要素特征,开展了三维高密度电阻率法、多源面波勘探和地面核磁共振方法等地球物理联合探测方法研究,结果表明：三维高密度电阻率法可用于圈定塌陷体边界;多源面波勘探对岩土体界面识别效果明显;地面核磁共振方法可量化地下含水量,划定含水层顶底板埋深,对岩溶发育程度和含水层孔隙度有一定指示作用;三种方法联合探测对提升岩溶塌陷体结构精细度具有重要意义。本文研究成果可为隐伏岩溶区地质结构精细探测方法选择提供参考。     

周期性四组元局域共振桩带隙特征及隔振性能研究

基于既有的三组元局域共振结构,考虑到岩土材料流塑性及弹塑性特征,提出一种由管桩、软体层、桩芯和基体组成的四组元局域共振桩结构,并基于周期性理论,利用频散曲线分析了管桩密度、弹性模量、厚度对四组元结构带隙特征影响。研究发现：相比于无管桩的三组元结构,带管桩的四组元结构更加容易出现完全带隙,且管桩只对带隙上边界产生影响;相比于管桩密度,管桩弹性模量可对四组元结构带隙产生质变影响,影响程度更大;通过四组元结构带隙边界点模态分析,确立了四组元结构带隙计算公式。通过时域分析发现,相比于三组元结构,四组元结构在带隙内振动衰减效果可提升近30%,其内在原因是四组元结构中管桩提高了基体整体刚度,从而有效地预防了基体“击穿”带及泄露波的出现。相关研究可为新型局域共振结构设计及隔振效果改善提供新思路。     

基于核磁共振的不同含水率黄土古土壤冻融循环试验研究

为明确不同含水率古土壤在冻融循环作用下其微观结构特性及古土壤损伤机制,采用核磁共振扫描仪对不同含水率冻融循环后的古土壤试样进行测试,研究冻融循环和含水率共同作用对古土壤微观结构的影响及土体内部损伤变化。结果表明:冻融循环下不同含水率使土体内部产生不同程度的损伤,损伤程度为含水率大的土体大于含水率小的土体。随着冻融循环次数增加,T₂谱曲线信号幅度增加,孔隙结构改变,大孔隙、最大孔隙含量增加,中孔隙含量减小; 同时含水率较大的土体孔隙体积增大幅度大于含水率较小的土体,说明在季节冻土区建设工程中含水率越大土体越易发生破坏,因此在工程中应注意防排水问题。依据损伤力学原理,得出土体颗粒连续性与孔隙率关系,进而得出有效应力与孔隙率关系; 根据核磁扫描结果,建立孔隙率与冻融循环次数关系,最终推导出古土壤有效应力与冻融循环次数关系表达式。研究成果为季节冻土区古土壤地层建设工程提供理论指导。     

碳封存砂岩储层的孔隙迂曲度分形特征研究

目前“双碳”减排的背景下,砂岩咸水层是最具潜力的二氧化碳封存介质,其孔隙结构特征决定了流体的运移及封存效率,对其进行研究具有重要的科学意义。分形维数常用于定量表征砂岩孔隙结构在三维空间中的分布规律,但以往对砂岩孔隙迂曲度分形特征进行研究时,多依赖于孔隙率、平均孔隙半径、平均迂曲度等特征值计算分形维数,并不能很好的反映岩石孔隙排列分布及孔隙连接的非均质性。文章将砂岩孔隙结构概化为具有迂曲度的毛细管,首先通过核磁共振成像（MRI）技术获取岩石片层图像,以盒维数法获得孔径分布分形维数;随后以孔径分布分形维数为目标值,结合不同迂曲度下的毛细管数量与孔径的分形标量关系,迭代计算迂曲度分形维数。与传统的分形维数计算方法相比,该研究确定的迂曲度分形维数更能体现岩石内部不同孔隙排列方式引起的迂曲度差异。     

Viscosity and nonlinearity effects on the forces and waves generated by a floating twin hull under heave oscillation

Nonlinear hydrodynamics of a twin rectangular hull under heave oscillation is analyzed using numerical methods. Two-dimensional nonlinear time-domain solutions to both inviscid and viscous problems are obtained and the results are compared with linear, inviscid frequency-domain results obtained in [26] to quantify nonlinear and viscous effects. Finite-difference methods based on boundary-fitted coordinates are used for solving the governing equations in the time domain [2]. A primitive-variables based projection method [6] is used for the viscous analysis and a mixed Eulerian–Lagrangian formulation [11] for inviscid analysis. The algorithms are validated and the order of accuracy determined by comparing the results obtained from the present algorithm with the experimental results of Vugt [22] for a heaving rectangle in the free surface. The present study on the twin-hull hydrodynamics shows that at large and non-resonant regular frequencies, and small amplitude of body oscillation, the fluid viscosity does not significantly affect the wave motion and the radiation forces. At low frequencies however the viscosity effect is found to be significant even for small amplitude of body oscillation. In particular, the hydrodynamic force obtained from the nonlinear viscous analysis is found to be closer to the linear inviscid force than the nonlinear inviscid force to the linear inviscid force, the reason for which is attributed to the wave dampening effect of viscosity. Since the wave lengths generated at smaller frequencies of oscillation are longer and therefore the waves could have a more significant effect on the dynamic pressure on the bottom of the hulls which contribute to the heave force, the correlation between the heave force and the wave elevation is found to be larger at smaller frequencies. Because of nonlinearity, the wave radiation and wave damping force remained nonzero even at and around the resonant frequencies – with the resonant frequencies as determined in [26] using linear potential flow theory. As to be expected, the nonlinear effect on the wave force is found to be significant at all frequencies for large amplitude of oscillation compared to the hull draft. The effect of viscosity on the force, by flow separation, is also found to be significant for large amplitude of body oscillation.     

太阳射电尖峰爆与电子回旋共振吸收

本文采用等离子体动力学方法，研究了日冕条件下磁化非相对论热等离子体对太阳射电辐射产生的电子回旋共振吸收，并在辐射频率等于电子回旋谐波频率时求得ｎ≥２谐波吸收率的近似表示式，以及其对等离子体温度，出射角度和谐波数的变化规律，在应用部分，讨论了电子回旋共振吸收对于太阳射电尖峰爆发的影响，认为目前观测到的尖峰爆，大多数高能电子束来自日冕的内层。     

基阶和高阶瑞雷波频散曲线反演对比分析

通过瑞雷波频散曲线反演,可以建立近地表横波速度模型。瑞雷波频散曲线具有多个模态,一般情况下,基阶模态是具有最大能量的模态。目前对瑞雷波频散曲线的反演大多数是基于基阶模态,对于高阶模态的反演比较少。通过对基阶和高阶瑞雷波频散曲线反演对比分析,认为高阶模态和基阶模态频散曲线反演有比较好的一致性趋势,相对而言高阶模态频散曲线反演的横波速度比基阶要高;高阶和基阶模态能量加权反演,不仅保留了频散曲线基阶模态能量强,勘探深度深的特点,而且对浅层低速层也具有更好的响应,可有效提高反演的精度。     

Loss-Cone Instability and Formation of Zebra Patterns in Type IV Solar Radio Bursts

The loss-cone instability of energetic electrons at double plasma resonance is considered. Conditions required for the formation of a zebra pattern in type IV solar radio bursts are determined. It is shown that electrons with a power-law energetic spectrum can effectively excite upper-hybrid waves at double plasma resonance. Stripes of a zebra pattern become more pronounced with an increase of the loss-cone opening angle and the power-law spectral index. The growth rate at the resonance frequencies decreases with an increase of the cyclotron harmonic number. Interpretation of observations and diagnostics of plasma for the April 21, 2002, event are performed. Conclusions about the impulsive mode of injection of energetic electrons into a coronal arc are made.