Using a thermal tracer to estimate overland and rill flow velocities

Flow velocity is a basic hydraulic property of surface flows and its precise calculation is necessary for process based hydrological models, such as soil erosion and rill development models, as well as for modelling sediment and solute transport by runoff. This study presents a technique based on infrared thermography to visualize very shallow flows and allow a quantitative measurement of overland flow and rill flow velocities. Laboratory experiments were conducted to compare the traditional dye tracer technique with this new thermal tracer technique by injecting a combined tracer (heated dye) into shallow flowing surface water. The leading edge tracer velocities estimated by means of infrared video and by the usual real imaging video were compared. The results show that thermal tracers can be used to estimate both overland and rill flow velocities, since measurements are similar to those resulting from using dye tracers. The main advantage of using thermography was the higher visibility of the leading edge of the injected tracer compared with the real image videos. Copyright © 2013 John Wiley & Sons, Ltd.     

利用多种地震数据联合反演剪切波速度结构的可靠性检测

本文模拟使用青藏高原东南缘区域台网及国家台网的170个宽频台站基于背景噪声、天然地震面波、P波接收函数反演时的实际数据,对青藏高原东南缘假定的初始模型进行恢复,通过计算初始模型台站下方纯路径频散、提取各台站对间的瑞雷波频散曲线、计算理论接收函数以及反演剪切波速度结构来测试使用不同单项数据与联合使用多种数据反演对初始模型的恢复程度。结果表明,同时使用接收函数、基于噪声经验格林函数的群速度、相速度频散以及基于天然地震面波的相速度频散联合反演的剪切波速度结构,充分利用了几种数据的分辨率优势,清晰地分辨出中下地壳及上地幔顶部的低速层。此外,本文也分析了实际数据处理中出现的计算误差、随机噪声干扰对计算结果稳定性的影响。结果显示:对于面波频散,加入1%的误差后,联合反演的结果仍可很好地反映低速层的形态,但是当误差提升至5%后,对最终结果则产生了一定程度的影响;而在接收函数中加入4%的随机噪声时,虽然地幔低速层的上界面和下界面会略微受到随机噪声的影响,但是低速层的深度范围和速度值均得到了较好的恢复。     

Anisotropy parameter inversion in vertical axis of symmetry media using diffractions

Diffractions play a vital role in seismic processing as they can be utilized for high‐resolution imaging applications and analysis of subsurface medium properties like velocity. They are particularly valuable for anisotropic media as they inherently possess a wide range of dips necessary to resolve the angular dependence of velocity. However, until recently, the focus of diffraction imaging or inversion algorithms have been only on the isotropic approximation of the subsurface. Using diffracted waves, we develop a framework to invert for the effective η model. This effective model is obtained through scanning over possible effective η values and selecting the one that best fits the observed moveout curve for each diffractor location. The obtained effective η model is then converted to an interval η model using a Dix‐type inversion formula. The inversion methodology holds the potential to reconstruct the true η model with sufficiently high accuracy and resolution properties. However, it relies on an accurate estimation of diffractor locations, which in turn requires good knowledge of the background velocity model. We test the effectiveness and applicability of our method on the vertical transverse isotropic Marmousi model. The inversion results yield a reasonable match even for the complex Marmousi model.     

LIMITS OF STRUCTURAL INVERSION OF SEISMIC HORIZONS*

Time horizons can be depth-migrated when interval velocities are known; on the other hand, the velocity distribution can be found when traveltimes and NMO velocities at zero offset are known (wavefront curvatures; Shah 1973). Using these concepts, exact recursive inversion formulae for the calculation of interval velocities are given. The assumption of rectilinear raypath propagation within each layer is made; interval velocities and curvatures of the interfaces between layers can be found if traveltimes together with their gradients and curvatures and very precise V_NMO velocities at zero offset are known. However, the available stacking velocity is a numerical quantity which has no direct physical significance; its deviation from zero offset NMO velocity is examined in terms of horizon curvatures, cable length and lateral velocity inhomogeneities. A method has been derived to estimate the geological depth model by searching, iteratively, for the best solution that minimizes the difference between stacking velocities from the real data and from the structural model. Results show the limits and capabilities of the approach; perhaps, owing to the low resolution of conventional velocity analyses, a simplified version of the given formulae would be more robust.     

Seismic inversion for underground fractures detection based on effective anisotropy and fluid substitution

Underground fractures play an important role in the storage and movement of hydrocarbon fluid. Fracture rock physics has been the useful bridge between fracture parameters and seismic response. In this paper, we aim to use seismic data to predict subsurface fractures based on rock physics. We begin with the construction of fracture rock physics model. Using the model, we may estimate P-wave velocity, S-wave velocity and fracture rock physics parameters. Then we derive a new approximate formula for the analysis of the relationship between fracture rock physics parameters and seismic response, and we also propose the method which uses seismic data to invert the elastic and rock physics parameters of fractured rock. We end with the method verification, which includes using well-logging data to confirm the reliability of fracture rock physics effective model and utilizing real seismic data to validate the applicability of the inversion method. Tests show that the fracture rock physics effective model may be used to estimate velocities and fracture rock physics parameters reliably, and the inversion method is resultful even when the seismic data is added with random noise. Real data test also indicates the inversion method can be applied into the estimation of the elastic and fracture weaknesses parameters in the target area.     

区域面波群速度反演的球谐函数法

一个定义在球面局部区域的复杂的面波速度函数如果直接利用球谐函数拟合可能需要展开到很高阶的球谐系数．通过保角变换，把一个球面局部区域扩展到球面上更大的区域上，变换过程中面波速度保持不变，在变换后的球面域上用球谐函数来拟合速度函数，达到降低球谐系数阶数的目的，使面波群速度的反演变成了球谐系数的线性化反演.通过球谐系数分析，可得到反演的分辨率.该方法不仅适用于面波群速度反演，同样适用于各种球面区域场的分析.     

汶川地震断裂带东北端浅部结构的人工地震探测

结合汶川地震断裂带动态监测,利用快速响应探测系统,开展了断层带浅部结构人工地震探测.针对地震断裂带动态监测条件下的复杂波场和低信噪比的情况,在f-k波场分离的基础上,分别利用了折射波共中心点成像、面波速度反演、反射波叠加成像方法,进行了浅层断层和构造成像处理,并对处理结果进行了综合解释,给出了断裂带浅部断层分布和速度特征.为汶川地震龙门山断裂带东北端动态监测提供了基础结构信息,所发展的断裂带快速响应探测技术对于地震应急动态监测具有重要意义.     

保角变换在区域面波群速度反演中的应用

通过保角变换，把一个球面局部区域扩展到球面上更大的球台域，再由球台域解析延拓到整个球面域.作为约束条件，变换过程中面波速度保持不变.在变换后的球面域上用球谐函数来拟合速度函数，达到降低球谐系数阶数的目的.面波群速度的反演变成了球谐系数的线性化反演，由球谐系数计算反演分辨核以及方差.该方法的优点是计算速度快，等值线光滑，构造界限清晰，该方法同样适用于各种局部球面区域的分析.     

高频面波方法的若干新进展

面波多道分析方法(MASW)通过分析高频瑞雷波确定浅地表剪切波速度.在过去的20年中,由于该方法具有非侵入性、无损、高效及价格低的特点,越来越受到浅地表地球物理和地质工程学界的重视,视为未来最有希望的技术之一.这篇综述论文将介绍中国地质大学(武汉)浅地表地球物理团队近年来在研究高频面波的传播理论和应用中取得的部分成果.非几何波是一种仅存在于浅地表介质,尤其是未固结的沉积物中的独特的地震波.它的存在对快速而准确地获得表层S波速度有一定价值.我们的研究表明非几何波是一种具有频散特性的泄漏波.泄漏波的存在可能导致将其误认为瑞雷波的基阶或高阶能量,从而造成模式误判.这种模式误判会导致错误的反演结果.我们通过求取高基阶分离后的瑞雷波格林函数证明虚震源法瑞雷波勘探的可行性.这个结果将极大地降低野外瑞雷波勘探成本.勒夫波多道分析方法(MALW)中未知参数比瑞雷波的少,这使得勒夫波的频散曲线比瑞雷波的简单.因此,勒夫波反演更稳定,非唯一性更低.勒夫波数据生成的能量图像通常比瑞雷波的清晰,并具有更高的分辨率,从而可以更容易地拾取精确的勒夫波的相速度.利用雅克比矩阵分析波长与探测深度的关系表明对相同波长的基阶模式而言,瑞雷波的探测深度是勒夫波的1.3~1.4倍;而两种波的相同波长的高阶模式波的探测深度相同.我们也尝试了时间域勒夫波反演.按照勒夫波分辨率将地球模型剖分成了不同尺寸的块体,利用反卷积消除了地震子波对勒夫波波形的影响,通过更新每个块体的S波速度来拟合勒夫波波形,从而获得地下S波速度模型.该方法不基于水平层状模型假设,适用于任意二维介质模型.     

High-resolution Electrical Resistivity Tomography monitoring of a tracer test in a confined aquifer

A permanent geoelectrical subsurface imaging system has been installed at a contaminated land site to monitor changes in groundwater quality after the completion of a remediation programme. Since the resistivities of earth materials are sensitive to the presence of contaminants and their break-down products, 4-dimensional resistivity imaging can act as a surrogate monitoring technology for tracking and visualising changes in contaminant concentrations at much higher spatial and temporal resolution than manual intrusive investigations. The test site, a municipal car park built on a former gasworks, had been polluted by a range of polycyclic aromatic hydrocarbons and dissolved phase contaminants. It was designated statutory contaminated land under Part IIA of the UK Environmental Protection Act due to the risk of polluting an underlying minor aquifer. Resistivity monitoring zones were established on the boundaries of the site by installing vertical electrode arrays in purpose-drilled boreholes. After a year of monitoring data had been collected, a tracer test was performed to investigate groundwater flow velocity and to demonstrate rapid volumetric monitoring of natural attenuation processes. A saline tracer was injected into the confined aquifer, and its motion and evolution were visualised directly in high-resolution tomographic images in near real-time. Breakthrough curves were calculated from independent resistivity measurements, and the estimated seepage velocities from the monitoring images and the breakthrough curves were found to be in good agreement with each other and with estimates based on the piezometric gradient and assumed material parameters.     

三维复杂地形近地表速度估算及地震层析静校正

在地表一致性模型的基础上提出一种可适用于宽线剖面、弯曲测线、传统的二维和目前广泛使用的三维地震观测.在地形及近地表低降速带地质结构复杂的探区，低降速带厚度及速度估算的精度是静校正处理的关键.本研究根据三维地震观测的初至走时数据，利用最小平方与QR分解相结合的算法，在三维空间重建近地表低降速带速度模型，根据重建速度模型实现了静校正长波长分量与短波长分量的同步计算.分析了复杂的近地表低降速带模型初至波的性质，在观测值的自动拾取以及理论值的计算中充分考虑了可能成为初至波的直达波、折射波和反射波的利用，提高了低降速带速度模型反演的精度.在初至走时观测数据的拾取中，本研究采用分形算法克服了初至波波形差异以及折射波相位反转导致的拾取误差，实现了三维初至拾取的大规模全自动化运算.在射线路径与初至波理论走时的计算中，本研究采用一种计算量与模型复杂程度无关的三维射线追踪方法，该方法以最小走时射线路径保证了与观测数据有同等意义的初至波的射线追踪及理论走时的计算.野外实际资料的处理结果表明了方法的有效性.     

Geostatistical integration of near-surface geophysical data

Accurate statics calculation and near‐surface related noise removal require a detailed knowledge of the near‐surface velocity field. Conventional seismic surveys currently are not designed to provide this information, and 3D high‐resolution reflection/refraction acquisition is not feasible for large survey areas. Satellite images and vibrator plate attributes are dense low‐cost data, which can be used in spatially extrapolating velocities from sparse uphole data by geostatistics. We tested this approach in two different areas of Saudi Arabia and found that the optimal recipe depends on the local geology.     

基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究

伴随层析成像（Adjoint Tomography）通过求解全波方程来准确模拟地震波在复杂介质中的传播,并利用波形信息来反演地下结构,是新一代的高分辨率成像方法.其中3-D伴随层析成像需要庞大的计算资源,而2-D反演相对更具计算效率.面波和远震体波是研究地壳上地幔速度结构的重要方法,它们对S波速度及Moho面的敏感度不同,通过联合反演,可以得到更为准确的S波速度结构及Moho面.通过两种数据的高度互补性,本文提出基于伴随方法的线性台阵背景噪声面波和远震体波联合成像方法,同时约束台阵下方S波速度结构及Moho面形态.我们将该方法应用到符合华北克拉通岩石圈典型结构特征的理论模型上,测试结果表明联合反演方法优势明显,相比于面波伴随层析成像,能获得更高分辨率的S波速度结构,同时能精准约束Moho面形态.相比于体波伴随层析成像,联合反演能有效压制高频假象,降低波形反演过程中的非线性化程度.本研究有望提供一种更为高效精准的线性台阵成像方法,搭建联合伴随层析成像理论框架,提升岩石圈成像分辨率,并为后续其他类型波形数据的引入提供思路和方法.     

基于卫星测高交叉点的海洋表面地转流速度

在流体静力平衡状态下,海洋Coriolis力和压力梯度平衡就形成地转流,世界上大多数海流都近似为地转流.本文利用卫星测高交叉点方法计算海洋表面地转流速度,分析了利用测高交叉点计算地转流速度的不确定性,上升和下降弧段的海面倾斜在分辨率50 km上可以达到10^-7量级,才可能获得优于10 cm/s的地转流速度.在低纬度或者纬度接近卫星轨道倾角的地区,由交叉点方法计算的地转流速度精度低于中纬度地区.以中国台湾东部黑潮为试验区,利用最新的中国台湾周边海域大地水准面模型参考场计算高精度的大地水准面高,利用TOPEX/Poseidon和Jason-1的GDR数据(2002~2005年)计算海面高,然后计算交叉点的动力高,确定交叉点的地转流速度,结果与中国台湾NCOR(National Center for Ocean Research)的流速基本一致.     

Assessment of Density‐Induced Tracer Movement in Groundwater Velocity Measurements with Point Velocity Probes (PVPs)

The concept of equivalent freshwater head was adapted to predict the conditions under which density‐driven flow would adversely impact measured groundwater velocities using point velocity probes (PVPs). Theoretically, vertical flow will result from any density contrast between the PVP tracer and the groundwater. However, laboratory testing of tracers with salinities ranging from 0 to 2000 mg NaCl/L showed that horizontal velocities could be determined with good accuracy with up to 60% of the total flow being vertical due to density effects in a gravel medium. The available data suggest that density effects are less likely to be pronounced in sandy sediments. The relative amount of vertical flow due to tracer density can be estimated from vertical and horizontal velocities measured with PVPs, or from the ratio of vertical to horizontal hydraulic gradients. The equivalent freshwater gradient produced from a given tracer salinity at 10 °C (a typical groundwater temperature at moderate latitudes) can be estimated from 7.80 × 10^?7 × (M_NaCl), where M_NaCl is the mass of NaCl added, in mg, to 1 L of site groundwater in the mixing of the tracer. Equations for other temperatures were also determined.     

基于谱元法的全波形反演及其在海底地震数据中的应用

全波形反演已被广泛应用于获取地下速度结构.而反演问题与正演方法密切联系,针对特定反演问题,合适的正演方法能极大提高反演效率和精度.本文首先验证谱元法在含起伏界面模型数值模拟方面的优势,在此基础上将谱元法作为正演引擎应用于全波形反演,并为克服未知子波的影响,采用一种归一化的频率域目标泛函.结果表明,起伏地表情况下,基于谱元法的全波形反演相比于基于传统有限差分法反演,具有更高的反演精度.进一步,本文将基于谱元法的波形反演方法应用于OBS观测系统的理论合成数据和野外采集数据.谱元法非结构化网格剖分自然满足自由边界条件,能很好地适应不规则海床并模拟多次波.理论实验表明,即使在OBS观测系统很稀疏的情况下,基于谱元法的全波形反演仍能获得海底以下正确的高波数速度结构.在处理实际OBS数据时,本文采用分频策略以减少反演非线性,初始模型成功更新,其结果揭示了西沙海槽海底以下更多的细节信息.     

Interstitial flow through preferential flow paths in the hyporheic zone of the Oberer Seebach,Austria

We investigated interstitial flow velocities in the Oberer Seebach, Austria, with NaCl tracer injections at a sediment depth of 30 cm to estimate the hydraulic conditions experienced by invertebrates inhabiting the hyporheic zone. Flow velocity measured with tracers is taken as travel time of the water along a straight line between injection and sampling points, although the water flows around sediment particles, and thus travels a somewhat longer distance. From sections of stream sediment in which the interstitial spaces were replaced by concrete, we estimated that this difference amounts, on average, to 27% and used this factor to correct the results of our velocity measurements. Corrected interstitial water velocities ranged from 0.01 to 1.32 cm s^-1 and were independent of surface discharge. We also studied spatial flow patterns in the bed sediments with long-term tracer injections. The three-dimensional distribution of tracer concentrations 24 hours after the start of the injection indicated that interstitial water preferentially flows in a complex network of areas of high hydraulic connectivity. Reynolds numbers for flow in the hyporheic pore space ranged from 0.1 to 489, implying that the flow environment varies from laminar up to the zone of transition to turbulent flow. Therefore, invertebrates may have a size-related active choice of areas where either friction drag or pressure drag predominates. The consequence of flow patterns, such as those observed in our study, is that small-scale variability of hydraulic conditions may be an important determinant of the patchy invertebrate distribution in bed sediments.     

Recursive impedance inversion of ground-penetrating radar data in stochastic media

The travel time and amplitude of ground-penetrating radar (GPR) waves are closely related to medium parameters such as water content, porosity, and dielectric permittivity. However, conventional estimation methods, which are mostly based on wave velocity, are not suitable for real complex media because of limited resolution. Impedance inversion uses the reflection coefficient of radar waves to directly calculate GPR impedance and other parameters of subsurface media. We construct a 3D multiscale stochastic medium model and use the mixed Gaussian and exponential autocorrelation function to describe the distribution of parameters in real subsurface media. We introduce an elliptical Gaussian function to describe local random anomalies. The tapering function is also introduced to reduce calculation errors caused by the numerical simulation of discrete grids. We derive the impedance inversion workflow and test the calculation precision in complex media. Finally, we use impedance inversion to process GPR field data in a polluted site in Mongolia. The inversion results were constrained using borehole data and validated by resistivity data.     

HYDROLOGICALLY DRIVEN MECHANISMS OF HEADCUT DEVELOPMENT

An analytic investigation of the effect of surface seal mechanical properties, overland flow, and subsurface hydrology was performed on headcut development. Headcut growth rates on upland areas have been observed to be quite small （less than 0.00015 meter per second） and that they occur in increments in which chips break off at points where cracks have developed in surface seals. The substrate soil under the seal collapses and is removed by the flow. This mode of headcut development is the result of a strong interaction between the surface and the subsurface processes. The surface process is energetically controlled by the mechanical features of the seal whereas the subsurface process is hydrologically controlled. The analysis yields estimates of the temporal scale of headcut velocities. In cases of infiltration from the vertical gully wall into the substrate, the flexural wave velocity （seismic sound velocity） was found to inversely affect headcut velocity.     

SOTEM数据一维OCCAM反演及其应用于三维模型的效果

本文基于垂直磁场分量研究了SOTEM数据的一维OCCAM反演方法,并将其应用于理论三维数据及野外实测数据的反演.对于大部分一维模型,OCCAM反演可取得较好的反演效果,且反演结果不依赖于偏移距;噪声对SOTEM数据的OCCAM反演具有较大影响,但当信号含噪水平不超过5%时,反演结果仍具有较好的准确性;若浅层存在较厚的低阻层,OCCAM反演结果对下部地层的分辨能力下降,仅能获得具有平均效应的电阻率.将一维算法应用于SOTEM三维数据的反演,会产生较大的误差,尤其是在异常体边缘地带影响最为严重.该影响程度与异常体和背景电阻率之间的差异有关,对于大多数电性近似呈连续变化的真实大地而言,一维OCCAM反演算法仍可获得较好的效果.最后通过陕西某煤田深部富水性调查的实测SOTEM数据反演验证了本文的研究成果.