基于HAFMM的无射线追踪跨孔雷达走时层析成像

本文使用最小二乘线性迭代反演方法对跨孔雷达直达波初至时数据进行反演,每次迭代过程中,用有限差分法求解走时程函方程,并用高精度快速推进方法(HAFMM)进行波前扩展,通过追踪波前避免了进行射线追踪.为了验证该方案,我们对三组合成数据进行了测试,分析了单位矩阵算子、一阶差分算子和拉普拉斯算子等三种不同模型参数加权算子对模型的约束和平滑效果;讨论了FMM和HAFMM对反演精度的影响;测试了LSQR,GMRES和BICGSTAB等三种矩阵反演算法的反演效果.此外,我们还对一组野外实测数据进行了反演,对比了基于本方案以及基于平直射线追踪和弯曲射线追踪的走时层析成像反演效果.对比分析结果表明,使用拉普拉斯算子和HAFMM进行反演能较好地进行目标体重建,而三种矩阵反演方法对反演效果的影响差别不大;并且通过对波前等时线图的分析可以定性地判断异常体的性质和位置;而在对实测数据目标体的重建上,本方案能达到甚至优于弯曲射线算法的重建效果.     

Ray-based amplitude tomography for crosshole georadar data: a numerical assessment

Analyses of travel times and amplitudes of crosshole georadar data provide estimates of the electromagnetic velocity and attenuation of the probed media. Whereas inversions of travel times are well established and robust, ray-based inversions of amplitudes depend critically on the complex directive properties of the georadar antennae. We investigate the variations of radiation patterns in the presence of water-filled boreholes and/or changes of electrical material properties in the vicinity of the transmitters or receivers. To assess the implications of such complicating factors for ray-based georadar amplitude tomography, we generate crosshole georadar data for a suite of canonical models using a finite difference time domain (FDTD) solution of Maxwell's equations in cylindrical coordinates. The emitting dipole-type antenna is approximated by an infinitesimal vertical electric dipole, whereas a corresponding receiving antenna is emulated by recording the vertical component of the transmitted electric field. Inversions of the amplitudes of these synthetic data demonstrate that the presence of water-filled boreholes as well as changes in the material properties along the boreholes may cause substantial artifacts in the estimated attenuation structure. Furthermore, our results indicate that ray-based amplitude tomography of crosshole georadar data is unable to constrain absolute values of attenuation. Despite these inherent limitations, the method is surprisingly robust at detecting and constraining relative changes in attenuation. In particular, we find the method to be highly effective for locating conductivity contrasts that are not associated with corresponding changes in dielectric permittivity, and hence, cannot be located by travel time tomography alone.     

Inversion of multiple intersecting high-resolution crosshole GPR profiles for hydrological characterization at the Boise Hydrogeophysical Research Site

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.     

各向异性TI介质qP反射波走时层析成像

地震走时层析成像是反演地层各向异性参数分布的有效方法, 但是关于地震各向异性介质走时层析成像的研究并不多, 其技术远远没有达到成熟的阶段.在野外数据采集时, 地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低, 利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法, 它可以对任意强度的TI介质模型进行反演, 文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演, 由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同, 所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数, 不同模型的走时迭代反演达到了较好的收敛效果, 与各向同性介质走时反演结果相比较, 各向异性介质走时反演结果具有较好的识别能力.     

基于一种反余切紧凑约束的电磁波走时层析研究

电磁波走时层析可用于空洞和孤石等离散块体异常的探测中.在传统电磁波走时层析中,通常采用一个低阶差分泛函来稳定反演过程,而这类泛函都具有光滑约束的性质,在成像结果中难以对此类离散块体异常进行准确的解释.本文基于Tikhonov正则化思想,提出了一种反余切泛函,该泛函具有紧凑约束的性质且无需选择一个聚焦因子.结合弯曲胖射线层析理论和重加权正则化共轭梯度反演算法,对两个理论模型进行了成像效果研究.与低阶差分泛函和最小支撑泛函的成像结果相比,该方法能够更好地突出异常的边界,成像结果中的虚假异常也要更少.同时,还分析讨论了激发电磁波的中心频率以及数据噪音对成像结果的影响.此外,针对一个地下连续墙墙体内部缺陷的现场模型,利用该方法取得了理想的成像结果,良好地揭示了缺陷的形态和位置.

     

多步法快速射线追踪与圆台型高精度走时外插计算

为更好地适应复杂构造的地震偏移成像,本文提出了一套快速射线追踪算法和一种高精度的走时外插计算方法.采用线性多步法的预测-校正公式求解射线追踪方程组,与传统的四阶Runge-Kutta法相比,提高了计算效率.在网格节点上的走时计算中,应用一种基于圆台的外插方法,该方法以射线的方向为轴确定圆台,将轴上的走时外插到圆台内的网格节点上.与传统的矩形体外插方法相比,圆台走时外插方法提高了计算精度,且具有更好的稳定性.另外,该方法利用稀疏分布的射线即可获得高精度的走时表,节省计算量,对复杂构造的偏移成像非常有利,尤其是三维偏移.最后通过逆散射偏移成像算例,验证了算法的有效性和适用性.

     

Time-lapse borehole radar monitoring of an infiltration experiment in the vadose zone

Ground-penetrating radar (GPR) is an effective tool for imaging the spatial distribution of water content. An artificial groundwater recharge test was conducted in Nagaoka City in Japan, and time-lapse crosshole GPR data were collected to monitor the infiltration process in the vadose zone. Since electromagnetic wave velocities in the vadose zone are largely controlled by variations in water content, an increase in traveltime is interpreted as an increase in saturation. In the test zone, the infiltrated water penetrated downward with an average velocity of about 2.7 m/h. A finite-difference time-domain method using two-dimensional cylindrical coordinates is applied to simulate radargrams associated with the advancing wetting front and to quantify the effects of critical refraction. Standard zero-offset profiling for which all first-arrivals are assumed to be direct waves results in an underestimation of water content in the transition zone above the wetting front. As a result, correct velocity analysis requires identification of first-arriving critically refracted waves from the traveltime profile to accurately determine a water content profile.     

基于声波方程的井间地震数据快速WTW反演方法

WTW(Wave equation traveltime+Waveform inversion)反演是基于波动方程的走时反演（WT反演）和波形反演的联合反演方法.WT反演利用波动方程计算走时和走时关于速度的导数,和传统以射线为基础的走时反演相比,具有不必射线追踪、不必拾取初至、不必高频假设以及初始模型和实际模型差别较大时也能较好收敛等优点,但WT反演与波形反演相比其结果分辨率低.与之互补的是,波形反演的反演结果分辨率高,但是当所给初始模型和实际模型相差太大时,波形反演迭代算法容易陷入局部极小点.可见结合两种方法的WTW反演是一种比较好的联合反演方法.常规WTW迭代算法是首先以WT反演为主反演得到地质模型的整体特征,然后再以波形反演为主反演模型细节,该算法耗时和占用计算机存储空间接近WT反演或波形反演的两倍.为了节省运算耗时和计算机存储空间,往往采取首先单独利用WT反演然后再单独利用波形反演的算法.这样做的缺点是不能紧密结合两种反演方法,使得它们的优缺点在每一次迭代中无法得到互补,从而影响了最终的反演结果.针对以上事实,本文提出一种新的方法实现WTW,使得WTW运算速度和存储空间在任何情况下等同于WT反演或波形反演.模型计算表明新的算法具有更好的收敛性.     

基于声波测井资料的波动方程走时层析反演速度建模

地震勘探目标逐渐由浅层（< 2000 m）转向中深层（2000~3500 m）,高精度速度建模是实现精确地震成像的关键手段之一.目前,初至波或早至波速度层析建模方法已较为成熟,然而对深层速度建模更为重要的反射波层析建模方法依然面临走时拾取困难及计算量大等难题.针对当前中深层走时层析速度建模方法面临的难题,本文提出了一种基于声波测井资料的波动方程全波走时层析速度建模方法及策略.该方法首先通过早至波层析反演获得浅层及大尺度速度信息;其次,通过提取初次成像剖面的构造倾角信息,利用声波测井资料进行稀疏反演,进一步获得更高精度的速度模型;最后,以此为初始模型,进行全波层析反演获得最终偏移速度模型.数值结果表明,本文所提出的速度建模方法及策略不仅可以避免传统反射波层析建模方法因采用Born近似反偏移所带来的大计算量问题,且大幅提高了中深层速度建模的精度,进而能有效地提高中深层复杂构造区域的地震成像质量.

     

不需提取激发脉冲的探地雷达波形反演方法

利用麦克斯韦方程的波形反演方法,能充分利用雷达数据中反映地下介质特性的波形、走时、相位等信息.它要事先已知激发脉冲信号,需从实际资料中提取,但受到诸多因素的影响.本文尝试避免从实际雷达资料提取激发脉冲的方法.首先,借助于一个已知地质模型上的探地雷达系统采集的记录以及给定激发脉冲、给定激发接收天线的正演模拟信号,计算出转换算子.其次,将转换算子作用于实测雷达信号,将其转换成给定激发脉冲激发的、给定激发接收天线采集的信号.最后,用转换后的数据进行波形反演.二维模型数据的转换表明,转换后的数据与正演模拟数据一致,反演结果差异甚小.实验室实测数据的反演结果与实际情况相符.     

多震相走时和波形三维地震成象方法及在祁连山地壳结构研究中的应用

本文的理论方法是以几何射线理论为基础发展起来的、天然地震走时反演技术及天然地震层析成像技术。它存在着震源函数与介质参数的解耦问题。本研究采取了五种方法来改善反演结果。包括,利用,Pg,Pn等震相增大约束条件;用已有精度较高的人工地震测深结果作速度约柬：用波形反演来修改模型,把诸多物理量开发出来互为约束,以修改后的模型再作反演,使解的稳定性大大提高：采用最优化过程,选择遗传算法。可以进行震源定位,走时反演,波形反演：得到任意深度的速度分布及从地表到Moho面的速度剖面。用于在祁连山地区的结果表明,这些层析剖面对认识大地构造、重大深部事件动力学是很有益的。     

弯曲射线追踪中Dijkstra算法的改进与实现

文章针对图论中寻找最短路径的Dijkstra算法内存占用量大,效率低的缺点,对该算法进行了改进,修改后的算法计算效率是原来的四倍,内存使用量和图中节点数呈线性关系.在此基础上,用新算法求出了激发点和接收点的最短走时路径,并由激发接收点的旅行时结合联合迭代法对理论模型和实际场地进行了反演.结果表明：和直射线追踪相比,弯曲射线路径能更好地反演出地质体内部的速度场分布.理论模型和实际探测结果证实改进后的算法是有效的.     

跨孔雷达全波形反演成像方法的研究

跨孔雷达全波形反演是一种使用全波形信息反演两钻孔之间地下信息的层析成像技术.常规的层析成像反演大部分采用射线追踪方法,其中基于初至时的射线追踪方法可以反演出速度剖面（介电常数剖面）,基于最大振幅的层析成像可以反演出衰减剖面（电导率剖面）.常规射线追踪方法有许多不足,究其原因是该方法仅使用了小部分的信号信息.为了进一步提高成像分辨率,本文全面推导了全波形跨孔雷达层析成像反演方法,该方法利用雷达波全幅度相位信息能够反演出地下高分辨率的介电常数和电导率图像.本文通过基于局域网的分布式并行算法,有效地解决了巨量数据正演计算问题.文中首先建立了基于单轴各向异性介质完全匹配层的时间域有限差分二维正演算法,进而通过应用包括时间维度在内的全波场信息与残场逆向传播的全波场信息乘积来计算梯度方向,通过求取以步长为自变量的目标函数的极值确定步长公式,并提出以第一次介电常数反演作为同步反演的初始模型,能够有效提高收敛速度.本文对多组模型进行成像实验,取得了较好的反演效果.     

弹性矢量波成像域联合层析反演

在地震弹性矢量波场框架下,推导了多波联合层析速度反演方程以及走时残差与角道集剩余曲率的转换关系式,提出了一种利用成像域角道集更新P波、S波速度的走时层析反演方法.其实现过程可以概括为:将弹性波多分量数据作为输入,基于高斯束实现矢量波场成像并提取角道集,利用层析反演方程求解慢度更新量,最终获得多波联合反演结果.模型试算和实际资料处理验证了该方法的反演效果,能够为弹性矢量波联合深度偏移提供高质量的叠前速度场.     

Bayesian frequency-domain blind deconvolution of ground-penetrating radar data

Enhancing the resolution and accuracy of surface ground-penetrating radar (GPR) reflection data by inverse filtering to recover a zero-phased band-limited reflectivity image requires a deconvolution technique that takes the mixed-phase character of the embedded wavelet into account. In contrast, standard stochastic deconvolution techniques assume that the wavelet is minimum phase and, hence, often meet with limited success when applied to GPR data. We present a new general-purpose blind deconvolution algorithm for mixed-phase wavelet estimation and deconvolution that (1) uses the parametrization of a mixed-phase wavelet as the convolution of the wavelet's minimum-phase equivalent with a dispersive all-pass filter, (2) includes prior information about the wavelet to be estimated in a Bayesian framework, and (3) relies on the assumption of a sparse reflectivity. Solving the normal equations using the data autocorrelation function provides an inverse filter that optimally removes the minimum-phase equivalent of the wavelet from the data, which leaves traces with a balanced amplitude spectrum but distorted phase. To compensate for the remaining phase errors, we invert in the frequency domain for an all-pass filter thereby taking advantage of the fact that the action of the all-pass filter is exclusively contained in its phase spectrum. A key element of our algorithm and a novelty in blind deconvolution is the inclusion of prior information that allows resolving ambiguities in polarity and timing that cannot be resolved using the sparseness measure alone. We employ a global inversion approach for non-linear optimization to find the all-pass filter phase values for each signal frequency. We tested the robustness and reliability of our algorithm on synthetic data with different wavelets, 1-D reflectivity models of different complexity, varying levels of added noise, and different types of prior information. When applied to realistic synthetic 2-D data and 2-D field data, we obtain images with increased temporal resolution compared to the results of standard processing.     

Stochastic estimation of facies using ground penetrating radar data

Explicitly defining large-scale heterogeneity is a necessary step of groundwater model calibration if accurate estimates of flow and transport are to be made. In this work, neural networks are used to estimate radar facies probabilities from ground penetrating radar (GPR) images, yielding stochastic facies-based models that honour the large-scale architecture of the subsurface. For synthetic GPR images, a neural network was able to correctly identify radar facies with an accuracy of approximately 90%. Manual interpretation of a set of 450 MHz GPR field data from the Borden aquifer resulted in the identification of four radar facies. Of these, a neural network was able to identify two facies with an accuracy of near 80% and one with an accuracy of 44%. The neural network was not able to identify the fourth facies, likely due to the choice of defining facies characteristics. Sequential indicator simulation was used to generate facies realizations conditioned to the radar facies probabilities. Numerical simulations indicate that significant improvements in the prediction of solute transport are possible when GPR is used to constrain the facies model compared to using well data alone, especially when data are sparse.This work was supported by funding to R. Knight under Grant No. DE-FG07–00ER15118-A000, Environmental Management Science Program, Office of Science and Technology, Office of Environment Management, United States Department of Energy (DOE). However, any opinions, findings, conclusions, or recommendations expressed herein are those of the authors and do not necessarily reflect the views of DOE. Further support was provided by a Stanford Graduate Fellowship to S. Moysey. The authors would also like to thank James Irving for his assistance with processing of the radar data.     

井间高分辨率纵横波层析成像研究井间油藏

采用视速度-偏振井间波场分离法分离得到纵波和横波波场,使用跨孔地震走时层析成像级联方法计算井间纵波和横波速度结构以及泊松比空间分布.综合纵横波速度、泊松比和测井参数,分析两井间岩性和储层特性.数值模型试验表明该方法垂直分辨率达到主频波长的1/4,并能准确分辨倾斜断层和垂直断层.处理相同纵波资料结果表明本文方法分辨率远远高于以往的井间走时层析成像方法,比井间声波层析成像分辨率更高.使用该方法处理垦71区井间资料,得到同台同源纵波和横波速度结构和泊松比分布.泊松比的低/高变化与测井参数指示的砂/泥岩层基本吻合.根据泊松比分布,参考岩石地震物理学测试数据,可以区分井间介质中的泥岩层、砂岩层,以及砂岩层中的饱水和饱油区,确定储层的连通性,圈定井间的剩余油藏.     

断控层析速度建模方法研究与应用

在复杂断块区域,速度建模及偏移成像是地震资料处理的重要环节,偏移成像的质量依赖于速度模型的准确度.为更好的提升复杂断块区域的速度建模精度,本文提出断层控制层析技术.该技术通过结构张量准确识别和提取断层位置,并建立断层控制矩阵和目标泛函表达式,基于上述方式实现对层析反演迭代过程的约束,最终获得高精度的速度场.为验证上述方法的有效性,应用常规速度层析方法与本文方法进行对比,用二者方法得到的速度分别制作合成地震记录,并与实际钻井数据合成地震记录进行相关,结果证实本文方法的数据和钻井数据的相关性为0.75,而常规层析建模方法的数据和钻井数据的相关性为0.58.结论表明本文方法的速度和真实速度匹配程度高,具有较好的应用效果.

     

The synthesis of travelling ionospheric disturbance (TID) signatures in HF radar observations using ray tracing

Characteristic signatures are often observed in HF radar range-time-intensity plots when travelling ionospheric disturbances (TIDs) are present. These signatures, in particular the variation of the F-region skip distance, have been synthesised using a ray tracing model. The magnitude of the skip variation is found to be a function of the peak electron density perturbation associated with the TID and radar frequency. Examination of experimental observations leads to an estimate of the peak electron density perturbation amplitude of around 25% for those TIDs observed by the CUTLASS radar system. The advantage of using the skip variation over the radar return amplitude as an indicator of density perturbation is also discussed. An example of a dual radar frequency experiment has been given. The investigation of the effect of radar frequency on the observations will aid the optimisation of future experiments.