Derivation Of Refractive Index And Temperature Gradients From Optical Scintillometry To Correct Atmospherically Induced Errors For Highly Precise Geodetic Measurements

Refraction effects of optical beams are generally caused by an inhomogeneous propagation medium and are a major source of systematic errors in the precise optical determination of angles and distances in the atmospheric surface layer. In this contribution a method for deriving vertical temperature and refractive index gradients from optical scintillation is presented. Knowledge of these gradients is required for the compensation of atmospherically induced errors for highly precise terrestrial geodetic measurements, like direct transfer and levelling. The advantage of the present optical method is, that temperature and refractive index gradients can be derived as line-averaged values over the propagation path, which is not possible by meteorological point measurements. Field observations have been carried out with a displaced-beam scintillometer over flat terrain and under different atmospheric conditions in order to verify this method. The experiments show, that this method allows to derive accurate correction values for precise terrestrial geodetic measurements.     

Effect of a Starting Model on the Solution of a Travel Time Seismic Tomography Problem

In the problems of three-dimensional (3D) travel time seismic tomography where the data are travel times of diving waves and the starting model is a system of plane layers where the velocity is a function of depth alone, the solution turns out to strongly depend on the selection of the starting model. This is due to the fact that in the different starting models, the rays between the same points can intersect different layers, which makes the tomography problem fundamentally nonlinear. This effect is demonstrated by the model example. Based on the same example, it is shown how the starting model should be selected to ensure a solution close to the true velocity distribution. The starting model (the average dependence of the seismic velocity on depth) should be determined by the method of successive iterations at each step of which the horizontal velocity variations in the layers are determined by solving the two-dimensional tomography problem. An example illustrating the application of this technique to the P-wave travel time data in the region of the Black Sea basin is presented.     

一种声波在固体和液体中传播的显示方法

给出了一种声波在固体和液体中传播过程的显示方法.该方法能够显示声波在二维平面内的动态传播过程.与动态光弹法、施利仑(Schlieren)法等声传播显示方法相比,该方法不仅能够同时显示不同时刻固体、液体中的声场,而且还适应于非透明固体;另外,对样品的残余应力没有要求,制作比较容易;但是,测量的声场是间接的,受接收换能器余震和液体直达波影响.将铝板放在液体里,记录到了在铝板中传播的纵波、横波和面波,在液体中传播的直达波以及3种首波,它们分别与铝板中传播的纵波、横波和面波相耦合,并且直观地显示了铝板与液体界面上所传播声波之间的耦合关系;将有机玻璃放在液体中,记录了有机玻璃中传播的声波及其与液体中所传播声波的耦合关系.      

A method for visualizing sound propagation in solids and liquids

Introduction When seismic waves propagate in a medium of complex geometry,such as a dam,it pro-duces stress concentration area at some spots of the dam because of the interface reflection and refraction of the waves,resulting in serious safety risk.Understanding seismic wave propagation in complex mediums becomes very important.In order to reduce the damage from natural earth-quakes,precaution measures should be designed on the basis of the seismic propagation charac-teristics.Experimental geo…     

二维变密度声波波动方程的衍射层析成像

本文将二维声波衍射层析成像原理推广到二维变密度声波波动方程中，导出了在线源激发条件下，不同观测系统的滤波反传播重建算法，取两个以上不同频率值可得到密度和体积模量的成像．数值实例表明，该方法能得到密度和体积模量较满意的成像。     

基于交叉梯度交替结构约束的二维地震走时与全通道直流电阻率联合反演

在利用不同的地球物理勘探方法对地下复杂介质成像时,因观测系统的非完备性及数据本身对某些岩石物性的不敏感性,单独成像的结果存在较大的不确定性和不一致性.对于地震体波走时成像与直流电阻率成像,均面临着成像阴影区问题.对于地震走时成像,地震射线对低速区域覆盖较差形成阴影区,造成低速区域分辨率降低.对于电阻率成像,电场线在高阻区域分布较少,造成高阻区域分辨率较低.为了提高地下介质成像的精度,Gallado和Meju(2003)提出了基于交叉梯度结构约束的联合地球物理成像方法.在要求不同的物性模型拟合各自对应的数据同时,模型之间的结构要求一致,即交叉梯度趋于零.为了更有效地实现基于交叉梯度的结构约束,我们提出了一种新的交替结构约束的联合反演流程,即交替反演不同的数据而且在反演一种数据时要求对应的模型与另一个模型结构一致.新的算法能够更容易地把单独的反演系统耦合在一起,而且也更容易建立结构约束和数据拟合之间的平衡.基于新的联合反演流程,我们测试了基于交叉梯度结构约束的二维跨孔地震走时和直流电阻率联合成像.合成数据测试表明,我们提出的交替结构约束流程能够很好地实现基于交叉梯度结构约束的联合成像.与单独成像结果相比,地震走时和全通道电阻率联合成像更可靠地确定了速度和电阻率异常.     

钱塘江双向声传播测流实验研究

2009年4月和5月在钱塘江杭州段进行了双向声传播测流实验.两个声信号发送与接收站设置在钱塘江两岸,间距为3050 m.在声传播实验期间,沿着声传播路径还进行了12次走航ADP(Acoustic Doppler Profiler)流速观测.根据双向声传播的顺逆流传播时间差计算得到沿声传播路径的平均流速,与ADP走航观测结果基本一致,两者的均方差为0.03 m/s.通过ADP走航期间双向声传播时间差的平均值和ADP走航观测所得断面流量之间的线性关系,得到了整个观测期间断面流量的时间序列.实验较好地捕捉到了因钱塘江涌潮而引起的流量变化,流量的最小值为 -3406 m³/s,最大为3509 m³/s.第一和第二实验阶段钱塘江的净流量分别为1329 m³/s和979 m³/s.以上结果表明用双向声传播测流方法可实现对河道流速和流量连续的监测.     

Rapid least-squares inversion of apparent resistivity pseudosections by a quasi-Newton method1

A fast inversion technique for the interpretation of data from resistivity tomography surveys has been developed for operation on a microcomputer. This technique is based on the smoothness-constrained least-squares method and it produces a two-dimensional subsurface model from the apparent resistivity pseudosection. In the first iteration, a homogeneous earth model is used as the starting model for which the apparent resistivity partial derivative values can be calculated analytically. For subsequent iterations, a quasi-Newton method is used to estimate the partial derivatives which reduces the computer time and memory space required by about eight and twelve times, respectively, compared to the conventional least-squares method. Tests with a variety of computer models and data from field surveys show that this technique is insensitive to random noise and converges rapidly. This technique takes about one minute to invert a single data set on an 80486DX microcomputer.     

电阻率层析成像的原理与初步应用

电阻率层析成像的基本原理是根据研究区域边界处的电压和电流值确定区域内部电阻率的分布,它同声波和电磁波层析成像的原理有着本质的不同。根据电阻率层析成像的数学原理,从积分变换、线性迭代反演和非线性迭代反演３种基本的反演途径系统地分析和归纳了当前几种主要的电阻率图像重建方法,并讨论了雅可比偏导数矩阵的求解。已开展的电阻率层析成像的野外观测试验,给出了实测的井间和地表两种观测系统的成像结果,分析了技术应用中的特点与限制     

Spectral variations of underwater river sounds

Passive acoustic monitoring of the self‐generated noise of particle impacts has been shown to be correlated to bedload flux and bedload size. However, few studies have concentrated on the role of acoustic wave propagation in a river. For the first time, the river environment is modeled as a Pekeris waveguide, where a wave number integration technique is used to predict the transformation of sounds through their propagation paths. Focusing on the distance of a hydrophone from the channel bed and cutting off the low frequencies produced by impacts between gravel particles, we demonstrate that acoustic propagation modifies the spectral content of bedload‐generated sound. Acoustic signals analyzed with the proposed model are interpreted by comparison to Helley–Smith bedload data obtained during flood conditions on the large gravel‐bedded Arc‐en‐Maurienne River, France. This study shows that careful attention to acoustic propagation effects is required when estimating bedload grain size distribution with hydrophones in rivers, especially for rivers with slopes higher than 1%. Bedload monitoring with a hydrophone is particularly appropriate for large gravel‐bed rivers – especially so during large floods, when in situ sampling is difficult or impractical and the impact of acoustic propagation is weaker relative to the self‐generated noise of bedload impacts. Copyright © 2017 John Wiley & Sons, Ltd.     

UK 1-D regional velocity models by analysis of variance of P-wave travel times from local earthquakes

A method is presented for deriving 1-D velocity depth models from earthquake bulletin data. The models can be used as initial models for more advanced modelling techniques such as tomographic inversion. The method is useful when there is little or no refraction and long-range reflection survey data. The bulletin travel times are subjected to an analysis of variance, where they are separated into source, distance, and receiving station terms. The distance terms describe the variation of travel time with distance, and the associated trend lines allow 1-D velocity models for the crustal layers to be determined. The velocity models provide an average crustal model for the region derived from local data. This does not include superficial layers which are necessarily poorly determined. Earthquake bulletin P-wave data from propagation paths across three different regions of the UK are employed to illustrate the use of the technique.     

耗散大气中的声波射线追踪

基于分层大气中声波的局地色散关系方程,建立一种计入真实大气衰减效应的有耗大气声波射线追踪模型.在色散方程的虚部中导出声波在运动大气中的耗散系数和竖直方向上的增长因子,并利用真实大气中的衰减理论对所得到的耗散系数进行修正.利用Hamilton方程组解出大气声波在考虑耗散效应下的射线微分方程组.该有耗射线追踪模型的数值模拟...     

Seismic velocities and thickness of horizontal and dipping subsurface layers

Summary A new approach to the determination of the interval and average seismic velocities and thickness of arbitrarily dipping subsurface rock layers from surface travel time measurement is presented. The method is referred to as thegeometric recursion. Its derivation is based upon the geometry of the seismic ray paths between a shotpoint and a geophone, and is basically an extension of Gardner's equation of dip correction. Travel time measurements at and distance between the shotpoint and the geophone and apparent dip angles of the layers are needed in using the techniques. The values determined with the new method are compared with those from two existing methods that are applicable only for the horizontal layers. The accuracy and shortcomings of these methods are then analyzed and discussed with real examples as illustrations. The main computing technique used in the method is the iterative approximation of the undetermined but interrelated subsurface parameters in common loops. The method is practical to use, particularly when the travel time measurements in an area are not redundant. Its efficiency is dependent primarily upon the accuracy of the travel time measurements as well as the effective control of the shortcomings of the method.     

Analysis of radar and raingauge measurements for a critical meteorological event in Tuscany

In this paper, some considerations are given to the employment of C-band polarimetric weather radars for rainfall estimates. The most common error sources are discussed, such as ground clutter and propagation attenuation effects, together with decorrelation in the sampling at the ground between radar and raingauge measurements, which can be quite significant in radar systems located in hilly regions, as is the case of the Arno basin in Tuscany. Since the main objective from a hydrological point of view is the estimate of rainfall at ground, integrations and comparisons are needed between radar and raingauge data, which are characterized by different time and space sampling. The paper is then focussed mainly on this problem and a technique is presented in order to improve radar based rainfall estimates through the integration with raingauge data, in order to enhance the correlation between the two types of measurements. Such a method is finally applied to a serious meteorological event which affected the Arno basin on October 1992.     

龙门山断裂带中段深地震宽角反射/折射剖面研究

射线追踪方法是研究地震波在横向非均匀地壳介质中传播的重要方法.本文推导了理论走时对网格化节点速度的偏导数公式,提出了针对深地震宽角反射/折射剖面数据反演的联合迭代法,并使用该方法对横跨龙门山断裂带中段的一条深地震宽角反射/折射剖面进行了反演和解释.首先,对每一炮的观测数据进行一维反演,在此基础上插值出一个粗略的二维速度模型;然后,使用射线追踪方法计算理论走时,再根据理论走时与观测走时的拟合程度对二维模型进行调整,以获得更加接近实际的二维速度模型;最后,利用联合迭代法对观测走时进行反演,经反复迭代使所有接收点理论走时与实测走时的残差平方和最小,最终获取该剖面的二维地壳速度结构.反演结果表明:测线东段的沉积盖层明显厚于中段褶皱带和西部高原,中部褶皱带部分地区出现基岩裸露;构造转换带两侧的地层分界面近于水平层状分布,其西侧的中、下地壳内各存在一个层间速度间断面;构造转换带内存在薄厚不等的低速层,自西向东有增厚趋势.此外,龙门山断裂带的3条主断裂向下深切结晶基底,这是由于西部松潘—甘孜地块自西向东运动,受到刚性扬子地块的阻挡,沿铲式断裂向上爬升所致;而在断层上盘距地表约15 km深处出现的最大剪应力极值区,正是发生汶川M_S8.0地震的震源位置.      

井间电磁测量的2．5维层析成像方法

利用正则化最小二乘反演方法实现了井间电磁测量数据的层析成像,对井间地层电阻率进行了重建。在成像算法中,我们假设了井间电磁的激发与接收采用电磁偶极子源,井间介质仅在二维（xoz）平面内变化。在数值模拟中,通过对构造走向（y方向）的Fourier变换,将三维电磁场问题转化为一系列二维问题,用等参有限元方法在波数域求解,使实际地层模型的处理得以实现。对于波数域中每个波数对应的电磁场方程采用等参有限元求解,并用高斯积分将波数域解变换为空间域电磁场。利用源与接收器电磁场的互易原理,实现了电磁场响应对电导率分布灵敏度的快速计算。针对正演模拟中源点的奇异性,我们采用具有一定面积的伪艿函数表达源电流分布,使数值解精度得到提高。用层状介质的解析解与数值计算结果的对比,验证了模拟算法的精度。用介质扰动产生的电磁场变化检验互易性定理计算灵敏度的有效性。对简单块状模型、斜向裂缝带模型及“大”字模型的模拟数据成像结果表明,本文介绍的层析成像方法是正确有效的。     

基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移

时间域常Q黏声波方程,由于含分数阶时间导数项,数值求解需要大量内存,计算效率低,不利于地震偏移的实施.通过一系列近似,可将该方程简化为介质频散效应和衰减效应解耦的分数阶拉普拉斯算子黏声波方程,数值求解内存需求少,计算效率高.本文采用交错网格有限差分逼近时间导数,改进的伪谱法计算空间导数,PML吸收边界去除边界反射,对该方程进行数值离散和地震正演模拟,开展地震数据的黏声介质逆时偏移,实现波场逆时延拓过程中同时完成频散校正和衰减补偿.改善深层构造的成像精度,数值结果表明,基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移可大幅度提高地震模拟计算效率,偏移剖面明显优于常规声波偏移剖面,极大改善深层构造的成像品质.     

S‐wave kinematics in acoustic transversely isotropic media with a vertical symmetry axis

Acoustic transversely isotropic models are widely used in seismic exploration for P‐wave processing and analysis. In isotropic acoustic media only P‐wave can propagate, while in an acoustic transversely isotropic medium both P and S waves propagate. In this paper, we focus on kinematic properties of S‐wave in acoustic transversely isotropic media. We define new parameters better suited for S‐wave kinematics analysis. We also establish the travel time and relative geometrical spreading equations and analyse their properties. To illustrate the behaviour of the S‐wave in multi‐layered acoustic transversely isotropic media, we define the Dix‐type equations that are different from the ones widely used for the P‐wave propagation.     

Acoustic finite-difference modeling beyond conventi-onal Courant-Friedrichs-Lewy stability limit: App-roach based on variable-length temporal and spatial operators

Conventional finite-difference (FD) methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy (CFL) numbers 0.707 and 0.577 for two-dimensional (2D) and three-dimensional (3D) equal spacing cases, respectively, thereby limiting time step selection. Based on the definition of temporal and spatial FD operators, we propose a variable-length temporal and spatial operator strategy to model wave propagation beyond those CFL numbers while preserving accuracy. First, to simulate wave propagation beyond the conventional CFL stability limit, the lengths of the temporal operators are modified to exceed the lengths of the spatial operators for high-velocity zones. Second, to preserve the modeling accuracy, the velocity-dependent lengths of the temporal and spatial operators are adaptively varied. The maximum CFL numbers for the proposed method can reach 1.25 and 1.0 in high velocity contrast 2D and 3D simulation examples, respec-tively. We demonstrate the effectiveness of our method by modeling wave propagation in simple and complex media.