Multi-resolution adaptive modeling of groundwater flow and transport problems

Many groundwater flow and transport problems, especially those with sharp fronts, narrow transition zones, layers and fingers, require extensive computational resources. In this paper, we present a novel multi-resolution adaptive Fup approach to solve the above mentioned problems. Our numerical procedure is the Adaptive Fup Collocation Method (AFCM), based on Fup basis functions and designed through a method of lines (MOL). Fup basis functions are localized and infinitely differentiable functions with compact support and are related to more standard choices such as splines or wavelets. This method enables the adaptive multi-resolution approach to solve problems with different spatial and temporal scales with a desired level of accuracy using the entire family of Fup basis functions. In addition, the utilized collocation algorithm enables the mesh free approach with consistent velocity approximation and flux continuity due to properties of the Fup basis functions. The introduced numerical procedure was tested and verified by a few characteristic groundwater flow and transport problems, the Buckley–Leverett multiphase flow problem, the 1-D vertical density driven problem and the standard 2-D seawater intrusion benchmark–Henry problem. The results demonstrate that the method is robust and efficient particularly when describing sharp fronts and narrow transition zones changing in space and time.     

利用重复测线校正海空重力仪格值及试验验证

通过分析海空重力测量系统误差的形成机理,我们发现海空重力仪格值标定误差是引起系统性测量偏差的主要因素之一.本文简要介绍了重力仪格值的标定方法,分析论证了格值标定的精度要求,提出了利用东西正反向重复测线检测校正海空重力仪格值的计算模型和补偿方法,分析讨论了该方法的校正精度及其适用条件,利用航空重力实际观测网数据对该方法的合理性和有效性进行了数值验证,证明该方法对消除海空重力测量系统性偏差具有显著作用.

     

利用模板匹配技术检测米尺度岩石断层黏滑实验中的声发射事件

声发射观测是岩石摩擦实验中研究强震孕育机制的重要观测手段之一.传统的声发射观测以触发记录为主,难以捕捉并提取微弱的声发射信号.近年来天然地震学中发展的模板匹配技术可以在连续波形记录中识别出微弱的事件,大大提升事件检测能力.本文发展了适用于多通道连续记录的声发射信号的模板匹配技术,并将该技术应用于分析米尺度岩石断层黏滑失稳过程.研究结果显示：模板匹配技术识别出的声发射数量约为传统方法识别的5倍,相对完备震级降低约0.3.临近失稳,检测目录完备震级以上声发射率表现出幂律增加的特征,且声发射事件聚集在失稳破裂起始点附近.失稳后,检测目录完备震级以上声发射率表现出幂律衰减,其中失稳后早期表现出较慢的衰减速度.另外,失稳后早期声发射事件表现出随对数时间沿断层迁移的趋势.最后,通过与天然地震学观测的结果对比表明,声发射模板匹配识别技术是研究大尺度岩石摩擦滑动失稳过程中前震和余震时空演化规律的有效手段,该技术助于在方法上为室内实验和野外观测进行融合研究提供思路.

     

Horizontal rough fracture type curves for groundwater movement

Abstract

Type curves are derived analytically for radial flow in rough horizontal fractures toward a well. The basic assumptions are that there is no turbulent flow near the borehole and the well storage is ignored. The basis of the methodology is to write explicit expressions for the continuity and cubic law flow equations, which are combined using a Boltzmann transformation leading to a simple ordinary differential equation for groundwater movement. Solutions are presented as a set of type curves for different fracture apertures. It is observed that the solutions provide a method of uniquely identifying fracture hydraulic parameters when the fracture is smooth, but pose ambiguity for rough fracture parameter estimations. However, large time portions of these type curves appear as straight lines on semi-logarithmic paper, which provides a unique way for rough fracture parameter determination. Identification of the fracture parameters, namely, the aperture and relative roughness, is possible in a unique manner with the use of these lines and the dimensionless time drawdown concept. The cubic law is the asymptotic behaviour, either for large times or large fracture apertures. Prior to this asymptotic part, there is a non-cubic portion which gives rise to systematic deviations from the cubic law. The technique presented is useful, especially for evaluating pumping tests from a single major fracture isolated by packers.     

日侧极光弧的发光强度与沉降电子能谱的相关关系

本文利用中国北极黄河站多波段全天空极光观测数据,选取稳定的日侧极光弧,统计研究了极光强度比I_557.7/I_630.0与极光发光强度I_557.7的相关关系.发现I_557.7在午前暖点和午后热点区附近出现极大值,分别为2.2 kR和2.9 kR;而I_630.0在磁正午出现极大值,为1.5 kR.当I_557.7从0.1 kR增加到10 kR时,极光强度比I_557.7/I_630.0也由0.2增加到9.结合DMSP卫星探测的沉降粒子能谱数据,找到17个DMSP卫星穿越黄河站上空极光弧的事件,共穿越40条极光弧.得到了沉降电子的平均能量正比于极光强度比I_557.7/I_630.0,沉降电子的总能通量正相关于极光强度I_557.7的关系式.利用该关系式反演所有极光弧的电子能谱,发现在午前和午后扇区,产生极光弧的沉降电子主要来源于等离子体片边界层;在高纬出现强度较弱的弧,对应等离子体幔区域.在磁正午附近,沉降电子的平均能量较低,极光弧处于低纬一侧,粒子源区主要是低纬边界层.

     

Fracture in Westerly Granite under AE Feedback and Constant Strain Rate Loading: Nucleation, Quasi-static Propagation, and the Transition to Unstable Fracture Propagation

New observations of fracture nucleation are presented from three triaxial compression experiments on intact samples of Westerly granite, using Acoustic Emission (AE) monitoring. By conducting the tests under different loading conditions, the fracture process is demonstrated for quasi-static fracture (under AE Feedback load), a slowly developing unstable fracture (loaded at a `slow' constant strain rate of 2.5 × 10<sup>−6</sup> /s) and an unstable fracture that develops near instantaneously (loaded at a `fast' constant strain rate of 5 × 10⁻⁵ /s). By recording a continuous ultrasonic waveform during the critical period of fracture, the entire AE catalogue can be captured and the exact time of fracture defined. Under constant strain loading, three stages are observed: (1) An initial nucleation or stable growth phase at a rate of ~ 1.3 mm/s, (2) a sudden increase to a constant or slowly accelerating propagation speed of ~ 18 mm/s, and (3) unstable, accelerating propagation. In the ~ 100 ms before rupture, the high level of AE activity (as seen on the continuous record) prevented the location of discrete AE events. A lower bound estimate of the average propagation velocity (using the time-to-rupture and the existing fracture length) suggests values of a few m/s. However from a low gain acoustic record, we infer that in the final few ms, the fracture propagation speed increased to 175 m/s. These results demonstrate similarities between fracture nucleation in intact rock and the nucleation of dynamic instabilities in stick slip experiments. It is suggested that the ability to constrain the size of an evolving fracture provides a crucial tool in further understanding the controls on fracture nucleation.     

断裂系扩展过程与地球物理场变化特征的实验研究

给出了判断断裂之间有无相互作用的力学判据.为了研究断裂间相互作用对地球物理场的影响,对若干类常见断裂系的扩展过程及扩展过程中的声发射b 值和波速变化进行了模拟实验研究.结果表明,各种断裂系由于相互作用的不同,扩展过程,声发射b 值和波速度化均有明显的不同特征.     

岩石摩擦-滑动中的声发射b值动态特征及其地震学意义

双剪试验的结果表明,岩石摩擦滑动过程中的AE b 值动态曲线可归纳为锯齿状和台柱状两类,分别对应于不规则粘滑和规则粘滑.滑动瞬间滑面上摩擦系数变化△μ的分布决定岩石的滑动方式,不均匀分布易发生不规则粘滑.均匀分布易引起规则粘滑;破裂和滑动是摩擦滑动过程中的两类AE 机制,两类AE 的能量分别与破裂面面积及△μ成正比.研究结果提示,平直光滑的地震断层可能显示出特征地震的发震特点.     

破裂岩石的电子发射与压缩原子模型

本文根据压缩原子模型提出岩石在高压及破裂过程中产生自由电子的一种可能机制。根据这种机制可在半定量水平上较为满意地解释岩石破裂中的电子发射、光辐射和岩石的爆炸式破裂。